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PREFACE 


This monograph serves multiple purposes. First, the editors feel 
that with increasing publishing costs and increasing demands from 
publishers for introductory books which sell to big audiences, there are 
only a few outlets for technical monographs. Yet since it is through 
technical papers that the field advances, the publication and distribution 
of monographs should not be totally sacrificed to economic necessity. 

Second, the editors feel that it is particularly important to publish 
the results of archaeological research. Archaeology is a destructive 
science. Once artifacts have been removed from the earth, they can 
never be replaced into the same context. Without published reports the 
information is lost to other specialists and to the general public. We 
therefore feel responsibility to report on Stanford's archaeological 
research in Mexico as rapidly as possible. This report presents new 
theoretical models and archaeological data which these models purport 
to explain. In addition, it reports on the successful use of new tech¬ 
niques for procuring archaeological information. 

Third, this monograph partially fulfills the legal obligation which 
the director of the research, the Department of Anthropology, and 
Stanford University contracted with the Federal Government of Mexico. 
Ezra Zubrow, as the director of Stanford's archaeological research in 
Mexico, and as a representative of the Department and University, was 
generously granted a permit to do archaeological research in Guanajuato 
by the Federal Government of Mexico^-. This permit, which was obtained 
from the National Institute of Anthropology and History under Article 31 
of the Federal Law on Archaeological, Artistic and Historical Monuments, 
requires the responsible parties to submit to the National Institute of 
Anthropology and History a variety of monthly, annual and final reports. 
This monograph is one of these reports. 

Fourth, this monograph represents an international, cooperative 
effort of both professional archaeologists and students from the United 
States and Mexico. Both professionals and students were involved in 
all stages of this research--the' development of the research design, 
the writing of the proposals, the field execution of the research, the 
writing and the final editing of the monograph. 

^Support for this field research project was provided under grants from 
the Research Development Fund, Stanford University, and Latin American 
Studies, Stanford University. 







The monograph consists of seven articles. They are: 1) The 
Archaeology of Canada,de Alfaro, Guanajuato, by Ezra Zubrow; 

2) Agricultural Patterning, by Ezra Zubrow, 3) Agricultural Dis¬ 
tribution in Guanajuato, by Thomas Sagehorn; 4) Archaeology and 
the Portable Refraction Seismograph, by Elizabeth Skinner; 5) The 
Relationship between Labor and the Size of Sites in Guanajuato, Mexico, 
by Elizabeth Skinner; 6) Analysis of Surface Survey Data, Canada de 
Alfaro, by Emily Wilkens; and 7) A Growth Model of Subsistence Cul¬ 
tures, by Keith Kintigh. 

The research and this monograph represent a large group effort 
and thus numerous individuals and institutions need to be thanked. 

Three individuals must be particularly noted: Senor Emilio Bejarano, 
the regional archaeologist for Guanajuato and Queretaro and the asso¬ 
ciate director of the research program, was indespensible in every 
research, logistic, and bureaucratic phase of Stanford's archaeological 
program; Senor Chavez Morado, the director of the Museo Ahlondiga 
de Granaditas of Guanajuato, was particularly kind in making the facil¬ 
ities of the museum available and providing laboratory and storage space; 
Mary Lou Johnson-Schwartz was the translator par excellence and was 
not only an important member of the research team, but a remarkable 
facilitator who made the contact between two cultures easily adaptable. 
Other residents of Mexico who provided numerous services included: 
Benjamin Brown, Patsy and Earl Jones, Doug Mellis, A.K. Desmond, 
Joanne Coyne, and Tita Braniff. We were fortunate to be able to assist 
four students from the University of Veracruz with their archaeological 
education and to receive the benefits of their experience. 

The following residents of the United States were members of a 
variety of research design, field research, and data analysis teams: 
Margaret Taylor, James Feurstein, Carl Hagenmeir, Robert Bone, 
Christine Hastorf, James Slater, Barbara Solomon, Amos Zubrow, 

John Juste son, Keith Kintigh, Robert Clemen, Elizabeth Skinner, Andrew 
Willard, and Nancy Botti. Drafting was done by Margaret Taylor, 

A.K. Desmond, Nancy Botti, Andrew Willard, and Kippy Nigh. Artist's 
reconstruction of sites were done by Doug Mellis. Field photography 
was the responsibility of Joanne Coyne and James Slater, while photo¬ 
graphic processing in the United States was the result of the patient 
efforts of Tom Sagehorn. 

Finally, I wish to thank the following institutions: Stanford Univer¬ 
sity's Research Development Fund, Center for Latin American Studies, 
Chiapas Project, Center for Research in International Studies, the 
Office of the Dean of Undergraduate Studies, and the Provost's Funds 
for Computation, all provided funds for various aspects of this research. 

In addition, the Department of Anthropology, the Computer Center and 
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the Linear Accelerator all provided equipment and consultation. The 
National Institute of Anthropology and History and the Museo Ahlondiga 
de la Granaditas generously provided information, staff, and supple¬ 
mental support, which included: field vehicles, office, laboratory, 
and storage space, as well as photographic laboratories and exhibit 
space. The Government of the state of Guanajuato generously allowed 
us to use the Treasury Department’s IBM 360 computer. The staff 
of Centenal helped us find and made copies of numerous aerial photo¬ 
graphs and maps as did the Departmento Hydraulico. The United States 
Geological Survey provided us with ERTS satellite imagery of Guana¬ 
juato; since these space observations were multi-spectral, they also 
provided us with the equipment to overlay the various photographs to 
produce a composite image. International Imaging Systems of Mountain 
View allowed us to use their electronic pseudo color density image 
slicing equipment in order to obtain the maximum information from 
our aerial and satelite photographs. The Ejido of Canada de Alfaro 
allowed us to work on their lands and in their fields as well as gra¬ 
ciously taking us into their homes and into their hearts. 


Ezra Zubrow 
Andrew Willard 
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THE ARCHAEOLOGY OF CANADA DE ALFARO: 
A REGIONAL DESCRIPTION 


Ezra B.W. Zubrow 
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Introdaction 

During the summers of 1972 and 1973 I directed the Guanajuato 
Prehistoric Regionalization Project in Canada de Alfaro, Guanajuato, 
Mexico. It was a joint archaeological expedition of Stanford University 
and the National Institution of Anthropology and History. The latter was 
represented by the Associate Director, Emilio Bejarano, who is the 
Director of the Regional Center for Archaeology for the states of Guana¬ 
juato and Queretaro. This article is a descriptive report of the results 
and an introduction to the volume. 

Objectives 

The project had numerous general objectives. First, we wished 
to do a regional study using the archaeological data as a central theme. 

It was clear from the outset that we could not limit ourselves, however, 
only to archaeological data. Therefore, one will find references in this 
paper and in the other papers in this volume to ethnographic, environ¬ 
mental, economic and other data. The second objective was to test a 
series of alternative, formal models. Models were developed to examine 
modern agricultural patterning and distribution as well as the relation¬ 
ships among prehistoric labor, site size, and site function. In addition, 
there was a continuing examination of models of economic growth and 
development in order to understand the processes of changing prehis¬ 
toric regionalization. 

There is a major gap in our knowledge of the prehistory of Guana¬ 
juato. No complete archaeological survey or generalized chronological 
sequence of the state exists. There has not been a complete excavation 
of a single site in all of Western Guanajuato, an area of approximately 
15, 250 square kilometers. It has been, estimated that there are more 
than 5, 000 archaeological zones within this western area. In short, 
this non-nuclear, marginal region of Mexico has been virtually ignored 
by archaeologists with the exceptions of important work being done by 
a limited number of scholars such as Tita Braniff, Emilio Bejarano, 
and Shirley Gorenstein. The latter two are presently undertaking a 
major project in the southeastern part of the state which will do much 
to correct the situation. Thus, the third objective was to complete a 
systematic site survey of a limited area within western Guanajuato and 
to begin the excavation of a series of sites with the purpose of helping 
to develop an accurate ceramic sequence and site typology. Fourth, 
we attempted to develop new techniques for interpretation of archaeo¬ 
logical materials. We developed a methodology for utilizing the re¬ 
fraction seismograph to determine the interior structure and stratigraphy 
of a variety of structures and applied it in archaeological zones. We 
used a variety of remote sensing imagery to provide additional informa¬ 
tion. This included ERTS satellite photography, aerial photography. 
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and pseudo color density image slicing, as well as analytical techniques 
for information processing. 

These general objectives were partially fulfilled by obtaining precise 
data from the following direct field objectives: 

1. The archaeological field objectives were: 

a. The exploration of the archaeological zone Canada de Alfaro. 

b. A systematic surface survey in order to map the locations of 
all the sites within the zone. 

c. The stratified sampling of the surface material from each 
mound within the zone. 

d. The mapping and photographing of each mound and the collec¬ 
tion of a wide range of information about each mound on site 
survey forms. 

e. The test excavation of a limited number of single and multi- 
component sites to help determine the site dimensions, the 
composition of sites, the number of components, and the nature 
and degree of cultural homogeneity of artifactual material. 

f. Artistic reconstructions of a limited number of sites. 

g. The analysis of aerial photographs to locate other archaeo¬ 
logical zones and sites. 

h. The centralized collection of information from other archaeo¬ 
logists on other known archaeological zones in the region. 

2. The objectives of the environmental field studies were: 

a. The collection of geological data and information through geo¬ 
logical surveys and reports to obtain stratigraphic and topo¬ 
graphic information of the zone and surrounding areas. 

b. The collection of biological materials for ecological informa¬ 
tion about potential carrying capacity, standing crop and 
productivity within the zone. 

c. The analysis of aerial photographs for information on agri¬ 
cultural distribution. 

d. Large scale agricultural field surveys to determine the 
patterning of agricultural fields. 

e. An ethnographic resource survey to determine the distribu¬ 
tion of natural resources used in agricultural manufacturing 
processes in modern subsistence economies in the area. 

f. Palynological sampling from sites to determine the prehistoric 
vegetation component. 

3. The objectives of the field laboratory were: 

a. The recording, cleaning and analysis of the surface and 
excavation artifactual collections from each mound. 
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b. The clarification of the chronological sequence of Canada 
de Alfaro. 

c. The clarification of the structural and architectural remains 
in particular sites of the Canada. 

d. The examination of previously predicted site locations from 
analytically determined patterns. 

e. The analysis of the biological and agricultural materials. 

Location 

one divides the map of M.exico in half, both vertically and 
horizontally, the point of intersection is very close to Guanajuato. 

The state is located on the central Mexican plateau about 3 to 4 
hours north of Mexico City by car. It is bordered by the states of Quere- 
taro, San Luis Potosi, Aguascalientes, Jalisco, and Michoacan. Most 
of the articles in this monograph concern aspects of Guanajuato geography. 
This article will emphasize the area surrounding Canada de Alfaro and 
secondarily the area nearby San Miguel de Allende. Canada de Alfaro is 
a small village, an ejido, and a hacienda located four kilometers east of 
Leon, the largest city in Guanajuato and approximately 35 kilometers 
northwest of the city of Guanajuato, the capital of Guanajuato. More 
precisely, the archaeological zone and the village is located at 101° 36.3' 
west, 21° 8' north (see Figure 1). The archaeological zone is located 
on the lands of the Ejido Alfaro, the Ejido San Pedro de las Hernandez, 
the Predio San Miguel, and upon the private property of Maria Solis 
Segura. The zone belongs administratively to the municipality (county) 
of Leon. It is 340 kilometers from the nearest ocean, the Pacific, and 
240 kilometers northwest of the famous site of Tula, as the bird flies 
(see agricultural patterning). 

The archaeological zones surrounding San Miguel de Allende are a 
hundred kilometers southeast of Leon and Canada de Alfaro. They are 
located between 100° 40' and 100° 52' west, and 20 45' to 21 north 

(see Figure 2). Zones are about 150 kilometers from Tula and about 
33 0 kilometers from the closest points on the Atlantic and Pacific sea- 
coasts. 

General Topography, Hydrology, and Geology 

Guanajuato lies between the two major chains of mountains which 
make up the backbone of Mexico, the Sierra Madre Occidental and the 
Sierra Madre Oriental. The bulk of the state consists of large moun¬ 
tain chains interspersed with large alluvial valleys. A traveler often 
has the impression that he is standing in a large flat bowl surrounded 
by high mountains. Several of these large valleys in the western and 
southern part of the state are interconnected and make up what is called 
the Bajio. These valleys are justly known as the breadbasket of Mexico. 


m 
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Figure 2. Archaeological zones near San Miguel de Allende. 



























Figure 3 is an ERTS satellite photograph of an area 186 kilometers 
by 186 kilometers (115 miles by 115 miles) taken from an altitude of 917 
kilometers (570 miles). It covers the area surrounding the city of Leon 
and a considerable area to the north and west. The cities of Lagos, 

Leon, Silao, and Guanajuato are marked on the photograph. The heavy, 
irregular black lines which have been drawn on the photograph demar¬ 
cate the large valleys, including part of the Bajio mentioned above. Most 
of the modern cities and prehistoric sites that are known are located 
within or upon the margin of these large valleys. The archaeological 
zone of Canada de Alfaro lies where the mountains meet the valley east 
of Leon. These valleys usually show the remnants of two or three river 
terraces as large individual hills or series of hills. A careful examina- 
Figure 3 will show the mountains and the drainages between these 
large alluvial valleys. The small, very irregular white shapes on Figure 
3 are lakes. 

The Bajio Valley near Leon, Silao, and the Canada de Alfaro archaeo¬ 
logical zone ranges in altitude from 1,750 meters (5, 740 feet) near Silao 
to 1800 meters (5,904 feet) near Leon, and at the base of the mountain 
range reaches 1850 meters (6,068 feet). The mountains consist of step¬ 
ped foothill ranges which slowly rise from 2,000 meters (6,560 feet) to 
2,100 meters (6,888 feet) to 2,200 meters (7,216 feet) and 2,300 meters 
(7,544 feet) as one moves east. Peaks reach as high as 2,5 00 meters 
(8,200 feet) at Mesa Gallo and Mesa Moreno, and 2,850 meters (9,348 
feet) at Mesa la Cuatralba. 

The state of Guanajuato has a reasonably complex hydrological 
configuration. The state makes up part of the Federal government's 
hydrological zone 12. The most important drainage is the basin of the 
river Lerma in the southern part of the state. It consists of a series 
of natural (Lago de Yuriria) and artificial lakes, such as the one behind 
Solis Dam. A set of secondary rivers, including the Rio Laja in the 
San Miguel area, the Rio Irapuato, the Rio Turbio, and the Rio Guana¬ 
juato have considerable drainages in their own right and create their 
alluvial valleys. The state has primarily a dichotomic drainage pat¬ 
tern with a major exception (see Figure 4). The area immediately 
surrounding Leon, Irapuato, Salamanaca, and Celaya have clear rec¬ 
tilinear drainage patterns. This clearly demarcates the Bajio from 
other valley areas. As one moves south from Leon to Irapuato and 
east to Celaya, the number of channels per square kilometer increases 
from an average of . 07 to .14 to .25. If one compares the area just 
west of the Canada de Alfaro archaeological zone to the area just east 
of the zone, a change from rectilineal to dichotomic drainage is reflected 
m a slight decrease numerically in the channels. This decrease, . 07 
to . 06 channels per square kilometer, also denotes that one is moving 
into a mountainous configuration. 
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Figure 3. ERTS Satellite photograph including major cities, 

contour lines demarcating valley areas and the Leon 
Silao transect. (NASA, Evos Data Center) 
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Figure 4 
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The area surrounding the archaeological zones in San Miguel de 
Allende also follow a dichotomic drainage pattern. However, in this 
area the dichotomic pattern is in combination with a parallel pattern. 
This area has approximately . 09 channels per square kilometer. All 
of the above figures are based on samples of 400 square kilometers. 

The geology of the large valleys, including the Bajio, is straight¬ 
forward. They consist of an alluvium made up partially from the valley 
beds themselves and partially from the mountains. As one moves closer 
to the mountains, increased amounts of vasicular basalt are found. The 
Bajio consists of middle Cenozoic extrusives (Cmv) derived from lava 
flows, Oligocene andesites, basalt and Miocene rhyolites. The andesites 
are usually found stratigraphically lower than the rhyolite deposits. The 
mountains which border the Bajio in the east near Leon, and which are 
partially incorporated into archaeological zones, consist of the lower 
Cenozoic intrusives (Cii) principally from the lower Tertiary. It also 
includes deposites of Grustaceous materials. As one moves further 
east into the higher mountains and the mining areas, that is closer to 
the city of Guanajuato, the Cenozoic intrusives give way to a complex 
pattern of Mesozoic Triasic deposits (TR). These consist of agglo¬ 
merates and calcareous deposits as well as black and green volcanics. 
These mountains have been and continue to be exploited for silver min¬ 
ing which dates back in the area to the 16th century. The silver ore is 
sulphur and sulphur-antimony and is found in pockets or enclaves. Pro¬ 
found exploitation has taken place over four centuries and the mine 
shafts reach 300-500 meters. 
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As one continues to move east from the city of Guanajuato and 
drops out of the mountains into the large alluvial valley surrounding 
San Miguel de Allende and its archaeological zones, the general valley 
geology returns to middle Cenozoic extrusives (Cmv). Just south of 
the town, however, there exists a component of Cenozoic volcanics 
(Csc) derived principally from basalt lavas, pyroclastic agglomerates 
and the numerous undercones of the recent Pliocene. 


General Demography. Biology, and Agriculture 

Guanajuato is a mixture of small, very productive agricultural 
regions such as the Bajio and large marginal zones located primarily 
in the mountains. It is our contention that prehistoric ally it was a 
marginal area except in localized zones. Although there have been 
some changes in the environment since prehistoric times, a quantita¬ 
tive understanding of the modern biosphere is useful for analytical pur¬ 
poses. For example, considering the state of agricultural technology, 
one may validly use modern demographic and agricultural statistics 
for the extreme upper limits that no prehistoric population could exceed 


Guanajuato, according to the 1970 census, had a population of 
2, 270, 370 people (1,139,123 males and 1,131, 247 females) in an area 
of 30,589 square kilometers for a density of 74.22 people per square 
kilometer. Only four other Mexican states out of the total of 32, in¬ 
cluding the Federal District of Mexico have higher population densities. 
The 1.56 percent of the surface area of Mexico which is Guanajuato 
supports 4.71 percent of the Mexican population. The 1, 800, 345 hec¬ 
tares under cultivation according to the 19 68 agricultural census rep¬ 
resents 1/3 of the entire land area of the state. If one removes the 
urban population (localities of 5,000+), the density drops to 38.17 per¬ 
sons per square kilometer. Adjusting for the amount of land under 
cultivation, the density for agriculturally productive land is approxi¬ 
mately 12.91 people per square kilometer.! These three figures (74.42, 
38.17, 12.91) bracket the range of densities calculated for early Neo¬ 
lithic populations by Allen for Southwestern Anatolia (24. 3 to 1.8), by 
Forge in New Guinea (34.7) and Zubrow in the Southwestern U.S. using 
floral productivity figures (60.5 - 27.7). Although the densities are 
similar, the productivity of the state is much higher than would be 
expected from a Neolithic population. 


How is this agricultural product distributed by crop? Table 1 
compares the agricultural acreage and production of Guanajuato and 
Mexico in 1968. Three generalizations may be drawn from this table 


1 One should note the slight inaccuracy built into this 12.91 figure 
by the combination of the 19 68 agricultural census with, the 1970 demo¬ 
graphic census. 
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Flrst, production is lower (an average 3. 3 tons per hectare) in Guana¬ 
juato than for all of Mexico (an average 4.2 tons per hectare). This 
marginality may be re-stated in other words that Guanajuato has 6.4 
percent of the total land in Mexico under cultivation but only produces 
5 percent of the total crop tonnage. If we could exclude the Bajio from 
these figures, it would show even more marginal situation. Second, a 
large number of Mexican crops are not grown in Guanajuato. Table 1 
reports all crops in the agricultural census for Mexico. 38 percent (26) 
of these crops are not grown in the state. Many of these are of course 
tropical. However, for even the major staples such as corn and beans, 
the state is just able to hold its own. It produces 5 percent and 7 percent 
of the total Mexican tonnage of corn and beans on 6 percent and 7 percent 
respectively of the total land in these crops. Third, one may examine 
the crops for which the environment in Guanajuato is particularly well 
adapted. By the term well adapted we mean either the amount or the 
area under production make up a significant proportion of the total Mexi¬ 
can cultivated tonnage area (over 30 percent). These crops are garlic, 
onions, green tomatoes, and strawberries. They are relatively unim¬ 
portant as a subsistence base, making up less than 2 percent of the 
state's agricultural land and less than 4 percent of Guanajuato's total 
agricultural product. These specialized crops make up less than 1 per¬ 
cent of the total Mexican land under cultivation and less than 1 percent 
of the total agricultural product. 

In order to obtain a more accurate understanding of the biosphere 
in the Canada de Alfaro archaeological region, a three part research 
design was used. The first part analyzed the patterning of agriculture; 
the second, the distribution of agriculture. They are reported in the 
next two articles in this monograph. The third part consisted of the 
collection and analysis of biological materials within the archaeological 
zone. The expedition began a systematic sampling study to analyze the 
floral biomass. The archaeological zone was divided into three sampling 
strata. The first strata consisted of the agricultural fields. The second 
strata was the non-cultivated areas of the valley portions of the archaeo¬ 
logical zone. The third strata consisted of samples taken in the canyons 
of the mountains bordering the archaeological zone. Five nested qua¬ 
drats (each consisting of a 30 meter square, 3 meter square and a 1 
meter square quadrat) were randomly located within each strata for a 
total of 15 nested quadrats or 45 individual quadrats. The thirty meter 
square quadrats were analyzed for plants of tree size, the three meter 
square quadrats for plants of bushes and shrub size, and the 1 meter 
square quadrats for plants of grass size. Christine Hastorf identified 
all the plants in each quadrat by popular name and Linnean name. In 
addition, she counted the number of plants and weighed them for biomass. 

Table 2 presents the number of plants grouped by genera and quadrat 
for each of the 15 samples in the 3 strata. The popular names are given 





Table 1. Agriculture of Guanajuato 
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Weight of plants in 1 meter 

sample (in grams) 33.0 29.5 50.1 40.9 32.1 56.6 125.2 108.9 59.4 89-9 275.1 31.3 66.2 234.7 69-4 
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for tree size and shrub size plants in addition to the generic name. 

For the grass sized plants, only the generic name is given but the total 
biomass per sample is also shown. Several generalizations may also 
be drawn from this table. 1) The diversity of plants clearly increases 
inversely to their size in all the ecological strata. Thus, as one moves 
from tree size plants to grass size plants, the average number of known 
varieties per sample for the agricultural strata increases from 1.0 (tree 
size), 2.2 (shrub size), to 3.4 (grass size). Similarly, for the non- 
agricultural strata and the canyon strata, the figures are 2.4, 3.4, 5.8, 
and 2.6, 3.2, 6.8 (see Figure 5-a). 2) The diversity of the flora in¬ 

creases as one moves from the agricultural to non-agricultural strata 
and from the valley to the canyon strata for plants of all sizes. Thus, 
the average number of known varieties per sample for tree sized plants 
increases 1.0, 2.4, 2.6 (agricultural, non-agricultural, and canyon 
samples respectively). Similarly, for shrub sized plants the figures 
are 2.2, 3.4, 3.2, and grass sized plants 3.4, 5.7, 6.8 (see Figure 
5-b). 3) The amount of biomass in the 1 meter square quadrats in¬ 

creases as one moves from the agricultural to the non-agricultural 
strata and from the valley to the canyon strata. The average amounts 
of biomass in grass sized plants per quadrat, are 37.1 grams, 87.8 
grams, and 135.3 grams respectively. 4) Agricultural products make 
up approximately 65 percent of all known plants in the 3 meter square 
quadrats. The ratio is 9 to 1 corn to beans. On the average, 2. 6 corn 
plants are planted per meter square. In the 1 meter quadrats where 
agricultural products are immature, they make up 5 percent of the 
known plants. There are 7.0 bean plants planted per meter square. 

In short, the general ecological picture for agriculture which may 
be drawn for Guanajuato is a picture primarily of marginality, one in 
which the land is similar to the average agricultural land in Mexico 
for the production of the staple crops and where production is propor¬ 
tional simply to size of land and agriculture. It is an area in which 
agriculture is a significant part of the biosphere (5 to 65 percent), but 
one in which there is more biomass in the non-agricultural and canyon 
zones than in the agricultural zones. Indeed, with larger biomasses, 
greater diversity of plants and thus stability, the non-agricultural and 
canyon zones must be considered an important resource for the prehis¬ 
toric populations even after agriculture was developed. 
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Present Sta te of Knowledge of the Prehistory of Guanajuato 

Archaeological interest has begun to focus on northern and western 
Mexico only during the last 30 years and with the exception of Porter's 
early work, it is only during the last 10 years that significant archaeo¬ 
logical research has begun to take place in the state of Guanajuato. 

The Lerma River, which flows through southern Guanajuato, serves 
archaeologists as a boundary line for the northern limits of "Nuclear 
Mesoamerica" as defined by Kirchoff. During the pre-Columbian era, 
nuclear Mesoamerica included several populations which shared many 
cultural traits. Most of Guanajuato north of the Lerma River is con¬ 
sidered a part of marginal Mesoamerica, a region characterized as 
being sparsely populated by hunting and gathering peoples with occa¬ 
sional agricultural settlements in river valleys. What is interesting 
to note about marginal Mesoamerica is, first, that the boundaries are 
relatively stable but do not always remain at the same location over 
time. Thus, Kelley notes a frontier with the Chichimecs between 100 BC 
and 1000 AD as far north as southwestern San Luis Potosi. This fron¬ 
tier could not be maintained and by the sixteenth century had retreated 
to the south and west. Second, under stress economic and cultural 
change appears to occur rapidly in these marginal areas. For exam¬ 
ple Armillas discussing nearby Michoacan claims that the cultures are 
short lived and isolated up to the advent of the Tarascan Empire. 

As far as I am aware, no early man sites located in Guanajuato 

T™ 611 6ither excavated or published. Litvak-King has described 
a Middle Formative culture in the lower Balsas basin of Michoacan. 

Fluted points have been found in Jalisco and Pleistocene fauna exist at 
Chupicuaro (see below) and along Lake Cuitzeo. 

Within the Preclassic Period (500 BC - 200 AD), a culture known 
as Chupicuaro (400 - 200 BC) flourished in a large area of southern 
Guanajuato and what is today the northern part of the state of Michoacan, 
especially along the Lerma River. From the evidence found at the 

upicuaro site (during Porter’s excavations done in 1946-48) and from 
• ^ Guana i uato Chupicuaro site reported by Matos Moctezuma 

^ Z lr Z 1S apparent that the Chupicuaro culture closely paralleled that 
oi the Valley of Mexico during the Preclassic period. Perhaps the Val- 
ey of Mexico was even surpassed if the quality, technique and variety 

are US6d aS a criteria - During the Classic period 
^UU - 800 AD) nuclear Mesoamerica expanded to the north. By the end 
of Teotihuacan IV, in the late classic, it had been enlarged over the 
Bajio and into the northwest to cover the area bounded by the 1800 meter 
and 2250 meter contours. Guanajuato and the other areas of marginal 
esoamerica incorporated in this expansion seem to have developed 
small ceremonial centers. Ceramics show local developmental trends 
wi certain influences from Teotihuacan. Marginal Mesoamerica does 
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not seem to have ever developed the cultural elements which prevailed 
at this time in the south -- the cult of Tlaloc, Quetzalcoatal, and other 
major dieties, the use of the Mesoamerican calendrical system, urban 
ceremonial centers, and many other traits which Kirchoff lists. 

During the Post Classic period nuclear Mesoamerica again con¬ 
tracts to its Lerma River boundary and nomadic hunting and gathering 
tribes move into areas where the ceremonial centers are deserted. 
During the Classic and Post Classic, it is probable that Guanajuato, 
and particularly the Lerma River Valley was a major route for the 
movement of goods and people. It is intermediate between the west 
coast, the Northern Frontier (i.e. the states of Zacatecas, Durango 
and Jalisco), and the Valley of Mexico, which served as an economic 
and cultural hub for most of Mexico in precolumbian times. 

One theory recently proposed by archaeologists working in the 
area is that not only was Guanajuato a migration and trade route but 
that as a border area it was actually a sort of military frontier with 
a series of fortifications to protect the civilizations of the south from 
the "barbaric, " warlike nomads of northern Mexico who were seen as 
a threat to the agriculturalists. During the summers of 1972 and 1973, 
work was being carried out in Acambaro in southern Guanajuato by 
Columbia University under the direction of Dr. Shirley Gorenstein. 

A structure was explored on a hilltop site which she felt was a possible 
fortification. 

Other recent archaeological research done in eastern Guanajuato 
includes the work of Beatriz Braniff, who has worked almost exclusively 
in the states of Guanajuato and San Luis Potosi, in an attempt to set up 
ceramic chronologies for these areas and correlate them with the valley 
of Mexico. Her work has provided a valuable framework for the study 
of ceramics in this project. Figure 6 is an attempt to look at the chron¬ 
ological relationships between Guanajuato and the surrounding areas of 
Michoacan and Zacatecas -- Durango. It represents a synthesis derived 
from the work of Chadwick, Kelley, Braniff, Gorenstein, Willey and 
others. It should be noted that at best it is a tentative chronology which 
outlines temporal relationships between cultural sequences, ceramic 
phases and archaeological sites. It must be considered tentative be¬ 
cause it oversimplifies an extremely complex geographically specialized 
situation and includes considerable speculation. 

The major works done in the entire marginal Mesoamerican region 
concerning specifically the study of settlement patterns are the ongoing 
studies of J. Charles Kelley. He has done extensive work in the states 
in Zacatecas, Durango and Jalisco on the Northern Frontier. 

The Canada de Alfaro zone is simply one link in a chain of inter¬ 
connected clusters of sites. One of the neighboring clusters was 
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Figure 6. Sequences in Guanajuato and Surrounding Areas 
(All of these sequences are actually far more complicated but have been simplified for ease of presentation) 
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1 S Z 7 S . temati f ally ex P lored and partially excavated by Emilio Bejarano in 
9 1. A circular structure somewhat reminiscent of Toate, Jalisco was 
explored and stratigraphic test pits were dug. The final report on this 
zone, known as Ibarilla (see Figure 1), has yet to be published. It is 
par cu arly interesting, in the context of Mexican archaeology, that 
circular structures are found in this region as they are uncommon in 
most of Mexico after the Preclassic except in the Huasteca, a region 
on the northeastern coast, and along the Northern Frontier. 

Archaeological Research in Canada de Alfaro 

After visiting various known sites in central Guanajuato (Tarimoro 
Agua Espinosa, Tierra Blanca, Ibarilla) in order to familiarize the re-’ 
search group with known sites, materials, and lithic and ceramic taxon¬ 
omies, archaeological research was begun in Canada de Alfaro. This 
zone was chosen for intensive study because it cross-cut two major 
eco ogical zones as well as numerous michrohabitats. Initially there 
was an informal exploration which consisted of attempting to delineate 
where the approximate outer limits of the zone were located. Once these 
lmits were determined, the research group turned to the systematic 
exploration of the zone. . 

Systematic Exploration of the Zone 

The zone was surveyed in the following manner. The total universe 
was defined and divided into two areas bisected by the river. The outer 
oundaries of the first area were then surveyed with the transit. A line 
oi ten people -- spaced approximately twenty meters apart -- was then 
set up, and this line was moved back and forth between the boundaries. 
The line began at the SW corner and searched a 200 meter strip until 
i reac e a oundary, then it turned and surveyed the next strip. This 
process continued until the entire area was covered and was then re- 

JxHfh m (W the SSCOnd area - 111 short ’ it was a systematic field survey 
with 100% coverage. When a prehistoric mound or activity area was 

oca e y t e survey team, it was tagged with surveyor's tape. Each 
si e was t en relocated by a mapping team which surveyed with a transit 
e exact location of each site as shown on the map of the zone (Figure7). 
The total area surveyed was 24.5 square kilometers in which were found 

I -fo* 18 constitutes a to t a l 2.9 sites per square kilometer. Of 
these, 18 sites are m area one and 54 sites are in area two. They 

range in size from 720 square feet to 80,435 square feet with an aver- 

““ ° f l 0 ' 169 s T uar ® feet. Within this group there are 48 mounds, 

lithic and sherd scatters, 6 plazas, 1 ballcourt. There are 24 circular 

an rectangular mounds as well as a variety of miscellaneous shaped 
site s • 

/i iS i arge diversity im PHes either that the archaeological zone was 
probably the result of multiple occupations or a reflection of multiple 
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activity areas. The map clearly shows that the distribution of mounds 
is non-random. Rather there is, in area 1, a generally linear pattern 
between site 3 and site 18. This is complicated by a major cluster of 
sites consisting of numbers 5, 6, 7, 8, 9, 10, 11, which I have labelled 
cluster 1. In area 2 there are 5 more clusters. Cluster 2 is located 
at the top of a hill and is composed of sites 19, 20, and 21. This cluster 
which includes a ballcourt, a site with a plaza, as well as two prehistoric 
roads, ranges in elevation between 1952 meters and I960 meters. Cluster 
3, the central cluster of the zone, consists of sites 22, 23, 24, 30, 31, 

43, 45, 46, 48, 50 to 59, and 60 to 63. This cluster is partially located 
on a large terrace. Sites 22 and 43 have large interior plazas. Cluster 
4 is sites 22 to 26; cluster 5 -- 32 to 37; and cluster 6 -- sites 3 4 

65 to 69, 73, 73b. * ' 

The reality of the clusters are calculated by comparing average 
inter-cluster site distance to average intra-cluster site distance. 

For example, the average distance between sites in cluster 3 is 
approximately 100 meters. The average distance between sites in 
cluster 5 is approximately 135 meters, the average distance,however, 
between the sites in cluster 3 and the sites in cluster 5 is 600 meters 
or approximately a ratio of 1:6. Similar results are obtained from the 
other clusters. If communication is a partial function of distance, it 
is clear that site communication within clusters is greater than site 
communication between clusters. 

If one examines the site distribution from south to north, it follows 
a general pattern of an upside down V. The central cluster ( # 3), makes 
up the vertex and the NE leg. It is also clear that the site density in¬ 
creases as one moves from south to north. This major site cluster is 
not surrounded on all sides by the other clusters. Instead the other 
clusters are on the southeast and southwestern sides and are not on 
the northern side of the central cluster. 

Stratified, Systematic "Quadrant Defined" Surface Sampling of Each Mound 

As each mound was discovered, a standard site form was completed 
and a systematic quadrant defined sample was taken from the surface of 
each site. First, the center of the site and the total area of the site were 
determined. The total area was sub-divided into three sampling strata, 
the top, side and border areas of the mound. A ten percent sample of 
each strata was obtained from five randomly selected surface squares. 

In order to obtain a representative sample, and also to insure that the 
sampling was not clustered, four of the five squares were randomly lo¬ 
cated under the constraint that one fell within each of the four compass 
quadrants (0-89, 90-179, 180-269, 270-359). The fifth square was ran¬ 
domly located anywhere within the strata. We thus sampled the full 
potential variation. Each square was located on the site map and then 
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a complete surface collection of all the cultural material within each 
square was made. All of the collected materials from each square 
were kept in bags labelled by site number, square number, type of 
material, date of collection, and name of the recorder. The total area 
collected from the 72 sites was 77, 337 square feet or an average of 
1,106 square feet per site. There were 41,927 ceramics and 9,088 
lithics collected. Only 5 of the 72 sites (19, 32, 33, 38, 40) have 
more lithics than ceramics. The lithic to ceramic ratio ranges from 
18:1 to 1: 29 with a total ratio of 1:5. 

If one calculates the lithic to ceramic ratios for clusters 1 to 6 
as percentages, the values are 37%, 88%, 18%, 11%, 52%, 18%. There 
are two groups. The first consists of clusters 1, 2, 5. It is more 
peripheral geographically and has higher lithic to ceramic ratios. The 
second more centrally located group, clusters 3 and 4, have lower 
lithic to ceramic ratios. However, the peripheral cluster 6 also has 
a lower lithic to ceramic ratio. It is possible to arrange these to em¬ 
phasize their exponential character. Figure 8 plots these ratios ex¬ 
pressed as percentages, as natural logrithims and as a base 10 logrithim. 

Why should these ratios be exponential? They do not correlate with 
the average size of sites per cluster, with the number of sites per cluster, 
nor with the chronology of the cluster (see chronology section). There¬ 
fore, a functional explanation is reasonable by excluding the other al¬ 
ternative s. 

More information derived from this survey will become relevant 
when we consider the chronology of Canada de Alfaro. A complete anal¬ 
ysis of the survey data occurs in the Wilkens article included in this 
volume. 

Detailed Surveying of Sites 22 and Site 43 

In addition to the field mapping of each site as part of the systematic 
field survey, detailed contour maps were drawn of sites 22 and site 43 
by transit. These sites were chosen on the basis of location, size, 
structural similarities, and excavation plans. Site 22 and site 43 are 
both located near the central intersections of the legs of the V (which 
makes up the distribution of the site). Site 22 is the largest structure 
in the zone; site 43 the eighth largest site (see Skinner n Labor and 
Size 11 in this volume). Both sites have similar structural characteris¬ 
tics. Each is a stepped mound which makes a platform on which a 
smaller pyramid, and a depressed plaza are located. Figure 9 is the 
contour map for site 22. The dimensions are dependent upon exactly 
where one measures the site for it is not exactly symmetrical. It is 
approximately 260 1 by 190* ■ From the base of the platform mound to 
the floor of the sunken plaza is 15 1 , to the surface of the platform 
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Figure 8. Lithic to Cermaic Ratio for Geographical 
Site Clusters Expressed as %, log^, and 
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mound to fee floor o£ the sunken plaza is 15', to fee ^faceof fee^ 
platform mound 20', and to the top of the pyramid is 30 . The sunke 

plaza is approximately 75' by 70'. 

Douglas R. Mellis, an artist and scientific illustrator with con¬ 
siderable experience in Mesoamerican archaeological sties, spen 
several days in the field and many hours with topographic maps an 
field notes. He attempted to recreate the site as accura e y * s p 
sible from a scientific and artistic perspective (Figure 10). He sh 
£e s«e partially uncovered and reconstructed in it. present^mnron-^ 
ment. There is clear archaeological evidence for al 
as the proportions, numbers of steps in the platform moun , 
curved sides). The only exceptions are the staircases an 
man on the platform mound in which he has used ar s c isc 

Fie’ure 11 is the topographic map of site 43 • The site ™' ore 
rectagular than site 22. It is approximately 240' by 280'. A 1 * 0 "® 

it covers a larger area, the vast majority of the sitems. not 
tall as site 22. From the base of platform mound (6038 j to^he sunk 
plaza is 1.7', to the top of the platform mound is 4 , and to fee t p 
the pyramid is 13’. 

From the survey, the south wall of the site showed two small step 

terraces on the platform mound (at 2' and 3'). upper terrace 

corresponds to an actual step platform confirmed by excavation. The 
location of the excavation trench and test pits is shown. 


Excavations; 

In site 43, the excavations must be considered expioratory. Oniy 
a surface area of 288 square feet or 4/10 of one peroent of &e total 
site area were excavated. The excavation may be chvrded into tw 
uarts Part one is the trench. It consists of six squares each 6 by 
^ Since the trench was placed on the edge of the mound, it was neces - 
ta‘rv to survey in a series of relative datum points. Each was of course 
tied into the site datum. Part two are the two test pits, square 9^an 
square 13, which are on the same east-west line as the ren <~ • _ 
trench and square 13 were excavated in a combination of artificia and 
natural levels. Each square within the trench and square 13 were 
excavated in 6" levels. Each 6” level was individually mapped by 
square and sorted for culhiral materials. Tbje.e materials ( were 
placed into separately labelled bags for each level of each squar e. 
phese were then washed and analyzed in the laboratory separately. 
However when there was an actual change in the geological strati- 

trtihcial level ended, and a new level which corresponded 

to the natural stratigraphy was begun. Thus, we have 
natural, cultural, and artificial stratigrap y. 







29 





Figure 10. Artist's reconstruction of Site No. 22, 
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Figure 11 „ Topographic Map of Site 43, Canada de Alfaro 
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Table 3. Summary Data from Excavations in Site 43 


Square Relative 
No. datum 

Approx. 
No. of 
art. levels 
exc avated 

Approx. No. 

of geological 

levels 

excavated 

Absolute 

depth 

excavated 

Number 

of 

ceramic s 

Number 

of 

lithic s 

i 

-5.36 

4 

2 

-24” 

1111 

166 

ii 

-4. 36 

6 

3 

-36" 

2594 

363 

HI 

-2.86 

10 

5 

-60” 

3789 

278 

IV 

-1.86 

13 

12 

-78” 

2313 

176 

V 

-1. 36 

14 

11 

-84” 

3589 

146 

VI 

- .86 

15 

17 

-90” 

3280 

205 

IX 

1. 64 

x (26)* 

12 

-155” 

1285 

269 

XHI 

5. 14 

27 

25 

-159” 

1643 

169 


After each square was finished being excavated, stratigraphic maps 
were drawn of all the walls to scale and mosaic photographs of the 
stratigraphy were made. 

•jU ' 

Square IX, the first test pit to be excavated, was not excavated 
by artificial levels. 


C14 samples were also taken where possible. 

The stratigraphy of site 43 is very complex. However, it is con¬ 
sistent to the extent that over 6' (from the north to south wall) of the 
excavation one is able to trace the same levels usually with small 
changes taking place. Figures 12 and 13 are simplified stratigra¬ 
phies of the south wall of test pit 13 and 9. Figure 14 simplifies the 
stratigraphy of the north waH of squares 3, 4, 5, 6. The platform 
mounds artificiality is most clearly shown in Figure 10 by the deep 
stone level. This level consists of closely fitted stones the entire 
horizontal plane of the square, (6' by 6') by 3 1/2' deep. Figure 12 
partially masks information about the step. What is not apparent from 
the cross section but which apparent from the horizontal sections of 
Square III is the large number of stones and their horizontal locations. 
As one moves down through the layers of the square, there is the fol¬ 
lowing distribution of stone. Two of the top levels (18” - 30” below the 
surface) have a large number (35) of construction size stones (20”+in 
diameter) in the western half of the square only. There are almost no 
stones m the eastern half. Between 30” -42” below the surface, con¬ 
struction size stones cover the entire square. Other indications of the 
step are the faced stones and stones placed upon each other to make up 
ttie wall . The site, however, follows standard Mesoamerican construe- 

usmg a wide variet y of material to fill the "rubble” 
interior oi the mound. 
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Figure 12. Site 43 -- stratigraphy of test pit XIII (south wall). 
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Figure 13. Site 43 -- stratigraphy test pit IX (south wall). 
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Figure 14. Site 43 -- stratigraphy of squares III, IV, V, VI (north wall). 



brown sand with gravel 
brown sand with no gravel 

brown sand with gravel 

light brown sand no gravel 
black clay with gravel 

black sand with clay 

dark brown soil 
brown soil with clay 

compacted brown soil 
sand 

compacted brown soil 





























































































-35- 


Other field analyses in Canada de Alfaro 

Two other types of field work were undertaken in the archaeological 
zone of Canada de Alfaro. First, a series of stratigraphic samples 
were taken. In order to do this we used a porta-hydraulic drill (Model 
DR-2000) made by Soil Test. Inc. The utility of this model I originally 
tested in shell midden sites in Northern California in conjunction with 
excavations being conducted by John Fritz of the University of California 
at Santa Cruz. This drill allows one to excavate down to 25' or more 
below the surface and to take out cores varying between 1" - 3” in dia¬ 
meter. The exact diameter of the core is dependent upon the bit and 
coring sampler used. Complete sequences (vertically) of samples were 
ta en on sites 10, 11, 22, 43, 57, 6l. We thus have a stratigraphic pic¬ 
ture of the sites to depths 44", 48", 82", 40", 57", and 94" respectively 
Since the actual diameter of the holes is so small, the disturbance of 
cultural material was kept to a minimum (less than 1/1000 of 1%). 

Artist and publication costs preclude actual publication of the se¬ 
quences here. However, the drilling showed that the stratigraphy of 
site 43 includes all the structural materials that were used to construct 
the other five mounds . Their order in each sequence varied. In addi - 
tion, the only other site with as complex a stratigraphy (i.e. number 
and variety of strata) is site 22. 

Second, a terra scout refraction seismograph (model R-150, Soil 
lest, Inc.) was used in conjunction with the drills to determine the 
stratigraphy. This machine is basically made up of three components, 

1) the hammer, which strikes the ground creating a sound wave; 2) the 
geophone, which picks up the sound waves; 3) the oscilloscope, which 
records the sound wave on a cathode ray tube and measures the travel 
time of the wave between the hammer and the geophone. It is capable 
of determining the nature of sub-surface materials. A dense material 
such as granite will transmit sound waves much faster than loose mate- 
nal such as suiid ( e .g. 10,000 - 20, 000 feet per second in comparison 
° feet per second). It is also capable of determining depth 

of subsurface materials since sound waves have the property of being 
refracted back to the surface when there is a high velocity material 
beneath a low velocity material. Multiple seismograph courses were 
run on sites 22, 43, 57, 61, 10, and 11 as well as the Rio Guanajuato 

With appropriate sampling, the internal shape of structures may 
also be obtained. A full report on this aspect of the analysis will be 
ound later m this volume entitled "Archaeology and the Portable Re- 
iraction Seismograph. M 
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Laboratory Analysis and Chronology 

While the surface survey and excavation were taking place, labor¬ 
atory analysis of the artifactual materials was begun. Every piece of 
ceramic or lithic material was classified and recorded by substance 
number and weight. The basic ceramic typology had 122 exterior ’ 

Z 3 de . C ° ratlon combinations (e.g. black on buff), 4 interior colors 

/f* U ^ an _ tem P ers (e.g. fine). A short summary of the total depth 
of the ceramic typology is presented in Figure 15. 

Figure 15. Summary of the Ceramic Typology 


Ceramics 



Artifact Sherd 

























-37- 


There is at least one example of every criteria for distinction, 
but not for every possible combination of distinctions. For example, 
there is at least one example of each of the 122 exterior color combina¬ 
tions, at least one example of each of the 4 interior colors, and at least 
one example of each of the three tempers. But there are not examples 
of 1464 combinations of distinctions (i.e. 122 x 4 x 3), instead only 
some 310 or 21% actually show up in all the collections from the sur¬ 
vey (i.e. 41,927 ceramics). A summary of the number of all ceramics 
from the 72 sites surveyed broken down by interior color and temper is 
below. It reveals a major hole in the combination of distinction actually 
occuring. In all the sites for all exterior colors, there are no ceramics 
with medium or course tempers and brown interior. In fact, there are 
virtually no ceramics with brown on the interior side of the vessel. 

Table 4. Summary of Ceramics from 72 Surveyed Sites 
Canada de Alfaro 


Interior Color 


Exterior 

Color 

Temper 

Orange 
or Red 

Buff 

Brown 

Black 

Total 

All 

fine 

6898 

9752 

2 

2364 

19016 


medium 

7895 

8800 

0 

3509 

20204 


course 

1125 

924 

0 

658 

2707 


Total 

15918 

19476 

2 

6531 

41927 


Potentially a piece of ceramics could fall into any one 184,464 cate¬ 
gories (158,112 categories for artifacts and 26, 352 categories for sherds). 

Lithics were classified into tools and debitage and then appropriately 
classified into types of tools and by materials. A short summary of the 
total depth of the lithic typology is presented in Figure 16. The total 
number of lithics by material from the 72 surveyed sites is given below. 

Table 5. Summary of Lithics from 72 Surveyed Sites 
Canada de Alfaro 

Lithic Totals (by material) 


All materials 

9088 

Chert 

5998 

Chalcedony 

2284 

Jasper 

406 

Obsidian 

398 

Basalt Rock 

2 

Quartzite 

0 
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Figure 16. Summary of the Lithic Typology 


Lithic s 



Tool 



Debitage 



Quartzite appears only in the excavated material. 

For either typology, the sites are analyzed by collection squares 
for the systematic stratified survey. For the excavation, each 6' by 6' 
square is classified by level. After the material was categorized, 
counted, and analyzed, the work continued at Stanford. (All the actual 
material remained in Guanajuato except for a 1% study collection which 
was returned after 6 months.) 

Maps of the area were sent up to the Stanford Linear Accelerator 
to be run on the Monroe Optical Scanners. From these Optical Scanners 
were received on card output the X and Y coordinates of all the sites. 

The coded ceramic and lithic data was combined with the X and Y coor¬ 
dinates of the sites and processed. The in-putting of the data into the 
computer was done by a series of Wylbur terminals. Wylbur, a general 
intersystem communication network, allowed for rapid inputting of data, 
as well as r&pid editing and correction. A series of cross tabulation 
matrices were run, by site, by type and by total. Since the result was 
120,000 lines of data, it is impossible to reproduce it all here, as it 
takes up a stack of computer print-out 10 inches tall. This is conceived 
as the first step towards the utilization of a more sophisticated analyses 
such as seriation, hierarchical clustering analysis, discriminant analy¬ 
sis, etc. However, since meaningful generalizations may be made from 
these cross-tabulations, it will be hopefully useful to present some of 
them here. For example, out of the 72 surveyed sites, 

1) 85% of the ceramics are non-distinguished buff, orange, 
orange-red, and dirty orange wares. 











- 39 - 


2) Only 4% of the ceramics are painted wares (e.g. Black on buff), 
but 57 different types occur. 

3) Deborah Schlacter has shown for 46 of the 72 sites that a 
reasonably simple typology of the sites based upon the number 
of ceramics from each site's 10% sample and the density of 
ceramics is identical to a variety of simple typologies based 
upon ceramic types. (See Tables 6, 7). She has shown that 

it is possible to use at least 41 other ceramic criteria and come 
out with exactly the same categories of sites. No other group¬ 
ing of the sites is capable of being created so many ways. In 
other words, this typology has the most redundancy through all 
the ceramics. 

4) As one moves from early to late through the relative chronology 
of the sites, the diversity of types of ceramics decreases and 
then increases. 

Table 6. Typology of 46 Sites Defined by Number of Ceramics 

in 10% Sample and Density of Ceramics per Volume Unit 
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Table 7. Typology of 46 Sites Using One Set of Minimum 
Ceramic Criteria. 


Black 



If one turns to the excavation, it appears that similar generaliza¬ 
tions may be made: 

1) As one moves from the early to the later stratigraphy, the di¬ 
versity of types of ceramics decreases and then increases. For 
example, in one square the number of types decreases 87% and 
then increases by 267% as one moves upward through the strati¬ 
graphic column. 
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2) Most of the ceramics are the same undistinguished wares as in 
the survey. For example, buff, orange, orange-red, and dirty 
orange wares make up 74% while painted ceramics make up 
2% of the ceramics in the square noted above. 

Tita Braniff and her assistants were kind enough to explain and show 
to me examples of pottery classified according to her typology for Grana- 
juato. The following table includes some of the diagnostic pottery types 
by phase. The pluses indicate significant quantities found at Canada de 
Alfaro. 


Table 8. Diagnostic Ceramic Typology by Phase 


Middle Preclassic 


Upper Preclassic 

Early Classic 

Phase a) San Miguel 

Phase b) 


Corporo gray-brown 
Negative 

Hematite redware 
Morales gray-brown 
Morales polychrome 
Morales red on buff 

Blanco Levantado 
San Miguel red on buff 
San Miguel black 
Rio Yerde orange 
Electra polychrome 


+ 

+ 

+ 

-f 

+ 


She also suggested that incized pottery on the interior of the vessels 
occurs primarily in the preclassic while incisions on the exterior occur 
primarily in the classic. Our analysis of a 1% sample (i.e. approxi¬ 
mately 4, 000 ceramics classified diagnostic or artificial) showed that 
for incized ceramics a ratio of 2:1 interior to exterior incisions. This 
suggests that the Canada de Alfaro zone dates to the preclassic-classic 
boundary. 

Turning to the lithics, several generalities appear appropriate. 

1) Most of the lithics in the excavation occur in square 6. No tools 
occur until level -12 and -13 where they are found in considerable 
quantity. Below -13 the number of tools decrease. From level -1 
down, the weight and quantity of each type of stone increases to -6 or 
-7, then it decreases to level -13. 2) The order of lithic quantities 

and quartzite in decreasing order which corresponds closely to what 
is found on all 72 sites surveyed. 

Obsidian hydration was used as a relative dating method. A piece 
of obsidian being exposed to the atmosphere will form a ring around its 
exposed edges by hydration, the longer it has been exposed and the older 
its age. In order to be useful as an archaeological dating technique, it 
must be used on a piece of obsidian which was fractured prehistoric ally, 
a tool or debitage. 
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Forty-five samples were prepared - one from each of 45 sites. 
Each sample showed prehistoric modification. Out of the 45 samples 
only 7 could be dated. Table 5 presents the obsidian hydration data. 

Table 9. Relative Chronology by Obsidian Hydration 


Time 

late 

a 


early 


Sites 

27 

49 

18 

3 

43 

10 

69 


Hydration 

measurement 

1.0 m 
1.1 m 
1.2m 

1.6 m 

1.7 m 
2.5m 
3.0m 


If one attempts to extend this chronology to other sites, an in eres - 
ing phenomena occurs which reinforces the chronology. . 
the sites which were grouped in the same ceramic grouping based on 

Table 6 as each site which had a obsidian hydration measurement. 

calculated the average diversity of ceramic types for each cer ^ 
grouping of sites with an obsidian hydration measuremen an 
ganized these values by obsidian hydration date (see Ta e ). 

Table 10. Tentative Relative Chronology of Sites by 
Ceramic Groupings 

Average Diversity 
of Ceramic Types 

18.0 


12. 3 
10.0 
11.2 
16 


I then took all of the sites in each geographical cluster which had a 
site with an obsidian hydration date and calculated similarly the average 
^r of ceramic typL per cluster. I then organized these averages 
bv obsidian hydration date. In this case also the diversi y 
types begins to decrease as one moves upwards through tec rono ogy 
aYd L/increases. This principle is exac*y what also is operating as 
one moves through the strata in the excavation of site 43. 


Time 

Sites by Obsidian 

Grouping of Sites 



Hydration Date 

Ceramic Categories 

late 

27 

8, 27 

/ 

V 

49 




18 




3 

15, 3, 2, 31, 45 



43 

43, 35 



10 

10, 29, 55, 58, 21 

e; 

irly 

69 

4, 1, 16, 59 
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During the excavation of Site 43 in Canada de Alfaro, several 
samples were taken from various areas and levels of the sample squares* 
These samples consisted of carbonized material. One sample, a small 
burned corn-cob, which was collected from Pit No. 13 was not analyzed 
by Carbon-14 dating. This sample was preserved by the Mexican Govern¬ 
ment as the only example of corn from this time period from the state 
of Guanajuato. Twelve other carbonized samples were taken to the U.S. 
for Carbon-14 dating analysis by Geochron Laboratories in Cambridge, 
Massachusetts. With the aid of Mr. Richard Pardi, the radiocarbon 
manager, I selected two samples of the twelve. These were the only 
samples which could be dated. I received a report on each sample. 

The essential data are as follows: 

Pre-treatment: Selected charcoal fragments were cleaned of 
foreign materials including rootlets, or other contaminating ma¬ 
terial that could be observed. They were then digested in hot 
dilute Hydrochloric Acid (HC1), and in hot dilute Sodium Hydro¬ 
xide (NaOH) to remove chemical contaminants prior to combustion 
and analysis. 

Sample No. 1. From Site 43, Trench 1, Pit 13, Level 1. 

Description: Charcoal fragment, 1/2 to 1/4 inch. 

The age is 1710 + 140 Carbon-14 years,B.P. (240 A.D.). 

Sample No. 2. From Site 43, Trench 1, Pit 13, Level 2. 

Description: Charcoal and bone marrow fragments between 1/16 
and 1/8 inch. 

The age is 1500 + 85 Carbon-14 years, B. P. (450 A. D.). 

It should be noted that at the 1 level, these ages are identical and 
that no particular importance should be given to their relative ages. 

These dates are based on the Libby half life (5570) for C 14 . The error 
stated is +_ 1 6 as judged by the analytical data alone. The modern 
standard is 95% of the activity of N.B.S. Oxalic Acid. The age is re¬ 
ferenced to the year A. D. 1950. These dates put site 43 at levels 1-2 
very close to the Preclassic Classic border as did Tita Braniff's ceramic 
chronology. 

In summary the laboratory analyses have not only shown a variety 
of generalizations about ceramics and lithics, but indicate that the 
Canada de Alfaro zone dates to the Preclassic Classic boundary approxi¬ 
mately A. D. 350. Furthermore, there is some evidence that not all the 
sites are contemporaneous based upon relative chronologies. 
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Other Zones 

Icist dccsfic "the focus in circhcicology 1ms chciugcd from, 
emphasis on locating cultures in time and space to questions of cultural 
process, resulting in more explicit theoretical formulations and testing 
procedures. One area in which considerable interest has developed is 
regionalization and geographical patterning. Using archaeological and 
ethnographic evidence dealing with site distributions, we attempted to 
construct a communication net with sites being the nodes, branches 
being direct communication links, and weighing such attributes as cost, 
desirability, and capacity. Then we used a "shortest tree" solution in 
order to determine the median of the graph and predicted the existence 
of other zones at the medians. For example, Canada de ALfaro lies 
near the median of the network between Ibarilla and the Cerrito de 
Rayas areas, as exemplified below: 



We attempted to predict a median site between the Rio Guanajuato 
sites (see below) and Ibarilla clusters. We underestimated the amount 
of weighing, and it was not until considerably later that we located the 
L-S zone. This zone which I have labelled the L-S site is located 5.8 
miles from Leon toward Silao along the Leon-Silao highway. (Measured 
from the Leon-Silao highway. (Measured from the Leon city limits sign). 
A track turns off the highway towards the east 2/10 of a mile from which 
the zone is easily visible. Someone has done considerable excavation 
previously. There is a 30 m by 3m trench running North-South and 
two large 4m by 4m excavations. All of these are nowin the process 
of gradually caving in. As of the time I am writing this I have not been 
able to find out who did the excavation and if it has ever been reported 
in print. 

Two seismograph courses were laid out and ran (see Skinner). 
Otherwise, this zone was left undisturbed. I am particularly interested 
in this zone since it is analomous being over a kilometer away from 
any primary drainage, which is a secondary pattern. This pattern 
postulates the correlation of zones with primary drainage. This hypo¬ 
thetical pattern holds valid for Ibarilla (partially excavated by Emilio 
Bejarano) and Canada de Alfaro. The Rio Guanajuato site was located 
on the edge of the first terrace of the Rio Guanajuato. It thus continues 
the pattern into the central Bajio. On this site a 10% surface collection 
was made for comparative purposes to the ceramics of Canada de Alfaro 
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and Ibarilla. In addition, three multiple seismograph courses were 
run (see Skinner). 

It is appropriate to close this descriptive regional study on what 
other archaeological zones are known in Guanajuato. A wide variety 
of zones have been examined by expedition personnel, Emilio Bejarano 
as the regional archaeologist, as well as by other interested archaeo¬ 
logists and individuals. I have attempted to centralize and standardize 
information. Table 11 presents a summary compiled by Mary Lou 
Johnson-Schwartz, Benjamin Brown, A1 Desmond, and Doug Mellis. 

It should be noted that this information comes from a wide variety of 
sources and people whose interest in the particular zones were varia¬ 
ble. Therefore, this information is tentative and is clearly unequal in 
completeness, value, or veracity. It does not include the intensive 
work being undertaken by Shirley Gorenstein whose results will soon 
be available in a far more complete, form. However, I present it for 
its potential value in stimulating more research in Guanajuato. Inten¬ 
sive research will surely modify this summary. 
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Agriculture as an Invasion 

It is possible to look at agriculture as the artificial introduction, 
an invasion, of a series of plants into an already developed natural 
ecological community. Considerable work has been done by various 
scholars dealing with invasions. Ecologists include such well-known 
names as Elton, Watt, Pielou and Mac Arthur. More closely, however, 
to the subject at hand, Binford and Flannery have used the concepts of 
artificial radiation, diversity, and the stability of environment for 
analyses of the origins of agriculture. Geertz and Rappaport have 
used these concepts to explore the relationship between cultural and 
environmental variables. By culture, I mean man's extra-somatic 
adaptation to his environment. By environment, I mean that part of 
man's natural world which he does not manipulate. From this per¬ 
spective the patterning of agriculture is a cultural mechanism for 
adapting to and controlling diversity, stability and productivity. 

Definitions 

In order to understand this theoretical stance a series of definitions 
will be needed. By invasion we mean the process by which an organism 
spreads from one, "the original, " community to a new community and 
successfully adapts to its new habitat. Invasions must be seen in the 
context of barriers and refuges. A barrier is any physical or ecolo¬ 
gical factor which inhibits dispersal or invasion and thus effectively 
acts as a mechanism of isolation. A refuge is an isolated habitat often 
protected by barriers which has preserved ecological (such as lack of 
predators) and/or environmental conditions (such as climate) which 
were more common in the past. Agriculture is the artificial introduc¬ 
tion of plants into a habitat of plants for human use, i. e. , an artificially 
stimulated invasion. By agricultural patterning is meant the repetitive 
contiguous succession of agricultural fields. 1 

By food chain we mean the transfer of energy through the consump¬ 
tion of one type of organism by another. For example, herbivores con¬ 
sume plants and carnivores herbivores. Although food chains vary from 
one to five trophic levels, there are rarely more than five because 
the amount of available energy for additional organisms is greatly 
reduced by heat loss at each successive consumption. However, ex¬ 
ceptions do occur since the efficiency of energy transferred varies 
enough for this to become possible. Thus, food chains may vary in 

’''For example, a transect showing the following sequence of fields, 
[corn, cotton, beans, corn, cotton, beans, cotton, cotton, sorghum, 
squash, corn] contains two patterns. One is [corn, cotton, beans], the 
other is [cotton]. Each is repeated twice. 
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complexity and in the number of different types of organisms which are 
dependent upon each other. 

In order to measure the complex series of related food chains 
which form an ecological community, the concepts of stability and 
diversity are necessary. By stability is meant the persistence through 
space and time of an ecological system. Diversity is a non-temporal 
measure of the distribution of the biomass. Most measures relate the 
number of species to the size of the population of each species in various 
ways. Diversity is related in two ways to community stability. First, 
by increasing the diversity of species of prey one adds alternate food 
sources for the predators. This enhances stability because there are 
more chains or routes by which the predators may be supported if there 
is a sudden decrease in one species of prey. Second, diversity among 
predatory species feeding on the prey will result in greater stability 
and control of the prey population. Changes in the size of one predatory 
population will be compensated by the size of the changes in the other 
predatory species, thus minimizing the change in the population of the 
prey. Later, it will be our contention that although development of 
agriculture is an invasion which takes place in communities with low 
diversity and thus short food chains, the patterning of developed agri¬ 
culture allows man to partially control the diversity of the community, 
the size of the food chain, and most importantly the productivity. The 
rationale for this viewpoint is explained in the following sections. 

Characteristics of the Invaders and the Invaded Environments 

The invaders, whether they are supposedly beneficial plants to 
man, such as agricultural strains,or harmful pests such as crabgrass, 
have certain common biological characteristics. First, they have the 
ability to adapt to a variety of new environments and are thus reason¬ 
ably generalized. Second, in comparison to most species they have a 
relatively high reproduction rate. Third, they tend to have competitive 
advantages which allow them to obtain exclusive occupation of their 
niches. Fourth, they tend to be relatively more transient than most 
species. 

Man's cultivation of plants is an attempt to create a monocul¬ 
ture, an artificially simplified community, by weeding out competitors 
and thus increasing production. This reduces the diversity index of 
the community. This initial competitor reduction is reinforced two 
ways. First, man's conscious attempts at herbivore control even at 
such a simple level as keeping animals out of the fields, reduces the 
range available to herbivores and decreases the number of species 
actually subsisting off the fields. Second, this change in herbivores 
is reflected in unintentional changes in the carnivore distribution 
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, Q c o x Th artificial instability of these communities are appar- 

(Elton observed the relatively rapid changes which take 

ent to anyone who has = fi alds a L abandoned. The history 

place in agncu ura selection has been a continual attempt to 

°£ agricultural breeding and 5electI °“ “ S - eld . This selection pro- 

i-nrrease the reproduction rate and thus, tne yi , 

“t often made species more specialised but there have also been 

periods of selection for more generalized species thus increasing the 
yield by allowing more land to be put under cultivation. 

What type of communities are invaders most likely to succeed 
and have a population explosion? The answer is natural commun^ 

with low diversity indices (l. e. , species poor comm uni 

Ire both fewer competitors and predators . 2 In addition, there are 

f must be met by the new environment. Not only 
r^rrirr:lentbTs S ^“ntly y benevolent as to allow the invader 
m survL. but it must also be possible for the invader to reproduce. 

Principles for Invaders 

It is now possible to state a series of principles involved in the 
successful invasion of foreign agricultural plants into a new community. 

They are: 

The stability of a community is directly related to the diversity 
of species to populations, the number of trophic links (i.e. , 
depth of food chain), the productivity of the species^ an 

predictability of the environment. S*|DI, TL, F, F j. 


1 . 


2 There are three types of species poor situations: new commun- 
the propertie^^f Ae 

ries are not mutoaUy gg ^ ^ predictabi iity areas will be subject to 

-npries noor communities at an early stage. . , 

P6 It P ill thus behoove us to briefly look at the c a 

teristics of new environments“ b ‘ a domin ant species 

new environment are usua y . , ,, , , as t b e community 

which is a general list and ie ^ o£ more specialized 

becomes more mature, it is repia y 

species which are "smaller niched. " 
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2 . 


The probability of a successful invasion is inversely related 
to the stability, diversity, productivity, predictability, and 

barriers of the recipient community, Io<.fi.» _L» J_* J_> J_> 1.1, 

y LS DI TL P Pr BJ 

and directly related to the stability, etc. , of the donor com¬ 
munity. 

The probability of a population explosion of the invading 
population is inversely related to the diversity of the reci¬ 
pient community and the number of refuges within the reci¬ 
pient community which allows competitors to survive 

P°T_L’ J_f. 

LDI R. 


Catch 22 and/or Hedging Man's Bets 

Previously we stated that agriculture is an artificially simplified 
community, ideally a monoculture with increased production because 
of minimal competition. However, this viewpoint is based on a sim¬ 
plified view of reality. First, even if there were a single ideal agri¬ 
cultural crop which would fulfill all of man's nutritional needs (i.e. , 
the ideal monoculture), man's penchant for novelty would probably 
lead to an increased number of crops with differing requirements. 
Reality, however, shows that even a large number of different crops 
are insufficient to meet man's nutritional needs. In short, diverse 
nutritional requirements cause diverse agricultural crops. Second, 
for man it is not enough that there is a successful artificial invasion. 

The size of the payoff in the crops must be sufficiently large to offset 
the varied economic costs, such as labor invested into producing the 
crop, etc. Third, for agriculture to be successful it cannot be a 
haphazard proposition. Over long time spans and in a variety of 
environments there must be a high probability of success in terms 
of large payoffs for a variety of crops. 

A farmer is faced by Catch 22. If he attempts to increase 
production by decreasing competition, the size of the payoff goes 
up, but the long term probability of success goes down. His overall 
nutrition may also go down because he is planting fewer types of 
crops (less diversity). If he attempts to increase the probability 
of success (as well as his interest in his diet) his payoff goes down. 
Conservative peasant communities are well known for opting for 
small, traditional plots, planted in a diverse set of crops. Agri- 
corporations with their larger savings base opt for land consolida¬ 
tion and monoculture. 

One may view these land use strategies as two ends of a continuum. 
Monoculture is a uniform distribution of one crop. Multiculture may be 


IfcS'lfnCf1 v I.CTRA3 
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defined as the most non-uniform (random) distribution of a large number 
of crops. Patterned agriculture is, of course, the middle part of the 
continuum — an intermediate number of crops partially patterned (see 
Figures la. and, lb. below). I use the terms mono and multiculture 
rather than simply mono and multicrops to distinguish that these strate¬ 
gies are the results of both the number of crops and the patterning of the 
land use. 


Figure la. 
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Figure lb. 

Agricultural Land Use Strategies 
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Recipient Communities and Land Use Patterning: Hypotheses 

Let us examine what kind of land use patterning one would expect 
under differing types of natural communities. As previously noted, 
diverse communities are poor candidates for population explosions 
of invaders, i.e., large production of agricultural plants (Principle 3). 
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Thus, in a diverse community the cultural problem is to increase 
production. One would thus expect a large amount of land use pat¬ 
terning which would culturally shift the system closer to monoculture 
and the desired increased production (Figure 2, Case A). In addition, 
since the recipient community is diverse, the total system is reason¬ 
ably stable and can tolerate greater variability in the amount of land 
use patterning. 

Under conditions of a non-diverse recipient community, the 
cultural problem is to increase stability of the recipient community. 
(One needs only note Principle 1. Stability is directly related to di¬ 
versity. ) Since the non-diverse recipient community is the ideal can¬ 
didate for the population explosion of agricultural plants, the potential 
for high productivity may already be assumed. The cultural solution 
is to decrease the amount of patterned land use (Figure 2, Case B). 

In addition, since the recipient community is not too stable, it cannot 
tolerate a large degree of variability in the amount of land use pattern¬ 
ing. 

Figure 2. 

Cultural Solusions to Diversity Problems 
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A Test Case 

The minimum requirements for a test case are two communities 
preferably in two environmental zones and the measured amount of 
patterned land for agriculture across these zones and communities. 
Although the two zones would guarantee differences in diversity, the 
transition between the zones theoretically would be more diverse than 
either of the homogeneous zones. During the summers of 1972 and 
1973 Stanford University, the University of Guanajuato, the Mexican 
National Institute of Anthropology and History have been conducting an^ 
archaeological research project in Canada Alfaro, Guanajuato, Mexico. 


3 This work financed by the Stanford University Research Develop¬ 
ment Fund, Latin American Studies, and a variety of administrative 
agencies of Stanford University. 
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As an adjunct to the archaeological data from 73 preclassic-classic 
sites, environmental and agricultural patterning data was collected 
on the ground along the entire length of a 180 mile transect. This 
transect began about 21 miles north of Mexico City and continued north¬ 
west to the city of Leon. 

Definition of Environmental Zones and Transitions: 

A cursory examination of the general environment showed that 
the transect crossed two zones. One zone, closer to Mexico City, was 
a semi-mountainous area with small dessicated valleys, the other con¬ 
sists of large alluvial valleys ringed by mountains. The boundary be¬ 
tween the zones is not a clear cut knife's edge, but a transition whose 
location is partially dependent on what environmental parameters; one 
is measuring. Table 1 presents the location on the transect of the en¬ 
vironmental transition points for a variety of environmental parameters. 

These are graphed spatially on the transect m Figure 3b. The zone 
and the transition are identified and located also m Figure 3b. 
transition is determined by the relatively larger number of transition 
points in the relatively smaller area. (Figure 3a. is included to orient 
the reader to the relationship between the graph's field-number axi 
and the actual transect. ) 


S. E. end of 
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Figure 3a. 

Explanation of Sequential Field Numbers 
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(Each unit on the axis equals 20 fields. Each 20-field group is located 
geographically (i.e. sequentially on the transect). 


Figure 3b. 

Location of Ecological Transition Points on Transect 
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Table 1 

Location of Environmental Transition Points on 



Land Use 

Transect 



Transition 

1. Koppens Climate 
Classification 

Transition 

From To 

Cwb Cwa 

Mile No. 

47 

Land Use 
Field No. 

213 

2. Mean annual 
precipitation 

75 0 mm 

600 mm 

41 

150 

3. Relative varia¬ 
bility of mean 
annual precipita¬ 
tion 

20% 

25% 

4 

16 

4. Relative varia¬ 
bility of annual 
precipitation 
anomalies 

-5. 0 

-2.5 

39 

13 9 

5. Surface geology 

quatenary 

alluvium 

quatenary 

limestone 

3 

12 


quatenary 

limestone 

tertiary 

extrusive 

35 

111 


tertiary 

extrusive 

quatenary 

limestone 

54 

274 


quatenary 

limestone 

tertiary 

extrusive 

109 

648 

6. Natural vegetation 
classification 

me squite - 
grasslands 

pine- 

oak 

41 

15 0 


pine- 

oak 

mesquite- 

grasslands 

66 

360 
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Measuring Agricultural Land Use : 

Along the transect land use was measured sequentially. Every 
time there was a change in crop or land use new measurements were 
made in the field. Thus, a segment of the transect might read: 


Field No. 457 Corn 100 m 

45 8 ^ Fallow 25 m 

459 Corn 50 m 

460 Non-agricul- 

tural Land 200 m 

461 Fallow 25 m 

462 Corn 25 m 


The data will be classified here into four categories: Corn (C), Fal¬ 
low (F), Non-agricultural (NA), and other crops (OC). 

It is possible to show part of the actual land use pattern along 
part of the transect in a clear visual manner. There exists a very 
powerful technique for analyzing remote sensing data called pseudo¬ 
color density slicing. It makes it possible to distinguish slight grada¬ 
tions in photographs. Figure 4 shows a black and white aerial photograph 
of the northwest end of the transect. Figure 5 shows the same photo¬ 
graph using a delay to emphasize the land use patterning. The pseudo¬ 
color density slicing technique makes it possible to demarcate different 
types of land use through the electronic mapping of different colors to 
show the differences in fields; e.g. fallow, non-agricultural, plowed, 
newly planted and planted agricultural fields. The figures above and 
the multiple color photographs were made possible by the exceptionally 
fine equipment and courteous help given me by International Imaging 
Systems of Mountain View, California. Cost, however, precludes 
their publication here. 


) 
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Figure 4. Aerial photograph of the archaeological zone of Canada de 
Alfaro and surrounding area including the northwest end 
of the transect. 



v.» 
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Figure 5. The area of the inset in Figure 4 using a 250 nanosecond 
delay to emphasize field patterns (negative). 
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The land use data was searched for patterns. All possible repetitive 
patterns up to an eight character segment were compared to the transect 
data (see Figure 6 below). 


Figure 6 

Comparison of Repetitive Pattern to Transect Data 
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The number of actual repetitive patterns per 20 fields has been plotted 
sequentially in Figure 7a. In other words, the number of repetitive 
patterns in the first 20 fields of the transect, sequential field numbers 
1- 20, were counted and plotted as a point; then, the number of patterns 
for fields 21 - 40, etc. Thus, as one moves along the graph from field 
numbers 1 - 20 to field numbers 981 - 1000, one is moving sequentially 
from one end of the transect to the other recording the amount of pattern¬ 
ing. A three point running average of the amount of patterning has also 
been plotted sequentially in Figure 7b. 


Figure 7 a. 

The Amount of Land Use Patterning per 20 Fields 
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No. of 
repetitive 
field pat- 
terns per 
20 fields 


Figure 7b. 

The Amount of Land Use Patterning per 20 Fields 
(3 point running average) 
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Results: 

If one examines Figure 7b., two facts are immediately apparent. 
First, the amount of land use patterning as one moves through space 
follows a cyclical pattern similar to a sine wave. Second, the size of 
the oscillation increases for the first two waves or 240 fields and then 
decreases for the next four waves. Comparing Figure 7b. to Figure 2, 
one sees that the ecological transition not only demarcates the separa¬ 
tion between the increasing and decreasing oscillation, but it corresponds 
to the greatest amount and variability of the amount of land use patterning. 

In short, it is clear that in the area of greatest ecological diversity, 
the transition, the amount and variability of land use patterning is great¬ 
est. As one moves further into the more homogeneous (less diverse) 
zones, the amount of land use patterning and the amount of variability 
in patterning decreases. This is exactly what is suggested by the hypo¬ 
theses and is what is expected if land use patterning is a cultural mech¬ 
anism to help maximize production and stability in a variety of ecological 
situations. 

Postscript: Hedging tlie Author’s Bets: 

The author wishes to hedge his bets two ways. First, I wish to 
point out that these distributions might be caused by other variables. 

One of the most obvious would be suggesting that land patterning is a 
function of the settlement pattern. Figure 10 shows the density of 
housing along the transect. There is no obvious relationship between 
land use patterning (Figure 7a. and 7b.) and housing density (Figure 8). 
Second, this theory and test case was based upon generally accepted 
ecological theory developed in the last 20 years by Robert Mac Arthur 
and others. However, it should be noted that a new book by Robert 
May entitled Stability and Complexity in Model Ecosystems (1973) has 
stood the more traditional theory on end. One of his critical argu¬ 
ments follows. 









-65- 


Figure 10 

Settlement Pattern Data 
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n The central point remains that, if we contrast simple few- 
species mathematical models with the analogously simple multi¬ 
species models, the latter are in general less stable than the 
former . . . 

"This straightforward mathematical fact contradicts the 
intuitive verbal arguments often invoked . . . , to the effect that 
the greater the number of links and alternative pathways in the 
web, the greater the chance of absorbing environmental shocks, 
thus damping down incipient oscillations. The fallacy in this 
intuitive argument is that, the greater the size and connectance 
of a web, the larger the number of characteristic modes of 
oscillation it possesses: since in general each mode is as likely 
to be unstable as to be stable (unless the increased complexity 
is of a highly special kind), the addition of more and more modes 
simply increases the chance for the total web to be unstable . . . 

"On the other hand, the balance of evidence would seem to 
suggest that, in the real world, increased complexity is usually 
associated with greater stability . . . 

"There is no paradox here. 

"One facet of the explanation can be that, although increased 
complexity makes for a more unstable system, it is advantageous 
on other grounds ... if there is a generalization, it could be 
that stability permits complexity. " 

Our model is only partially dependent upon "the stability diversity 
axiom. " It is not critical for this paper that complexity causes stability 
(the traditional viewpoint) or that stability permits complexity (May's 
viewpoint). In either case, the cultures will attempt to use land pat¬ 
terning as a mechanism for maximizing both invader production and 
probability of long term success (specialized types of decreased di¬ 
versity and increased stability). 
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Archaeological sites are subject to continuous modification by man. 
Therefore it is necessary to study man's effect on the present environ¬ 
ments well as the prehistoric changes in the environment m order m 

place archaeological data in the proper cultural contact, 
primary modification in the environment surrounding the Guanajuato 
sites is agriculture, three models were developed to better define the 
contemporary agricultural systems operating in the site-clus er are . 

The Bajio region near Leon is generally characterized as an extremely 
fertile and productive region of Mexico. However, as one travels thro g 
this region it becomes apparent that certain areas, most notably alo g 
the hills that flank the valley, exhibit a subsistence or marginal type of 
agriculture. Next to the cluster of excavated sites, an m par icu , 
adjacent to site-19, an ejido (subsistence) farming community called 
Canada de Alfaro exists with a population of 450. . A P pr03 f™ at ® h ^ ^ ; or 
meters north of this village the regional market is located. This J 
market area, Leon, has a population of 364, 990. The village of Canada 
de Alfaro is located at the flanks of hills. Leon, while also located nea 
the Mils, is nearer the valley-center such that these population centers 
“em an excellent test case to evaluate some effects of 
and social variables on field location, producfivtty, and 

that influence the 'cash-crop versus subsistence-crop' dectstons a farmer 
makes in this area* 

The aerial photographs used in this study cover an area of. 128 
kilometers (see the 'Agricultural Land Map' on page' 81 

cation of photographs). A systematic analysis of totalled ‘ 

cultivation 1 number of fields per grid-unit, transportation ne ’ 

and Variability between fields according to 'total land(under 
(techniques explained later in the article) was carried out with respect 
to the specific requirements of each model. 

Agricultural Models for Canada de Alfaro and Leon 
Model One : 

This model hypothesizes that there should be highly marginal land 
in the site-cluster area due to the ejido farm's low productivity. T 
model further contends that low productivity results directly from 
tually reinforcing political, historical and economic factors. 

Codigo Agrario y Ley de Colonization (a government publication that 
dafin® am! legaiVequLements for agriculturai activity h » 

known that an ejido farm collective cannot extend beyond a 7 kil 

radius of the ejido center. Canada de Alfaro, lo £ at ^ ^iTs^b sistence 

site-cluster, is a community of farm-workers who work at 

levels. The ejido farm is characterized by individually tilled, sm , 
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subsistence plots* In these farm units it is also assumed that mecha¬ 
nization is not widespread and agricultural technology is low. 

When expropriation of hacienda land occurred in Mexico, the legis¬ 
lative codes (e. g. Codes of 1934 and 1942) assured that the former 
landowner could retain the more fertile land to cultivate. In the Canada 
Alfa.ro area the ejido land is located near the hills, away from the 
fertile valley center (see Figure I below). The haciendado (the land- 
owner), in short, offered the poorest quality land for expropriation and 
retained the best irrigated land surrounding the casco, the center of the 
estate where the more fertile soil is found. 

Figure 1 

Site-cluster 
Canada de Alfaro 

© 

Marginal Land _ 


Hacienda 

© 


Valley-center 


Distance 


Highly Productive 
Fertile Land 


This model also contends that economic factors inhibit the growth 
of ejido farm technology and the ejidatarios (the ejido farmer) ability to 
acquire more land (possibly near the central valley). The ejido farmer 
is incapable of countering the use of marginal land and low productivity, 
because credit is difficult to obtain. Carr ( 1969 ) points out that the pri¬ 
vate banks never favor small-scale agricultural credit. The Banco 
Nacional de Credito Ejidal only offers the ejidatario credit on the se¬ 
curity of his crop rather than on his land. However, the land is collec¬ 
tively owned so that only ejidatarios coming together as collectives can 
get credit. Belshaw (1967) contends that the lack of social cohesion 
inhibits the cooperation necessary to obtain credit. In the village case 
of Huecorio, a village Belshaw studied, located 110 miles due south of 
Leon, even small groups of 10-20 farmers would not pool their resources 
and apply for credit since none were willing to accept the collective re¬ 
sponsibilities of an immediate debt and undefined returns. 

As a result of (1) political and socio-economic restraints placed on 
the ejido farmer, and (2) the highly marginal farm land in the site- 
cluster Canada de Alfaro area, two other subsistence farming charac¬ 
teristics can be seen to inhibit agricultural production: 

1. It is possible that an innovation will enable the farmer to produce 
more corn in place of supplementing his crop with purchases, with the 
effect of releasing sufficient cash from other activities to cover the cost 
of innovation. Belshaw maintains that this type of logic ignores the fact 
that the cash income arrives largely in the summer from vegetable crops 





- 70 - 


(authors note: if cash income arrives at all). At this point in the season 
all accumulation of corn has been consumed. Consequently, the capital 
is needed to support the family until the following corn harvest in Decem¬ 
ber. In short, economic conditions define the action that the farmer 
takes. In the Canada de Alfaro area, characteristics of subsistence 
farming impose a heavy penalty for failure whenever an attempted inno¬ 
vation does not succeed. A failure which places the farmer below the 
margin of subsistence or in debt with a local creditor could prove counter¬ 
productive in the long-run. 

2. Mellor (1968) contends that even if ejidatarios are motivated to 
change and innovation is profitable, it is necessary that the individual 
himself participate in the net increase in wealth which accompanies suc¬ 
cessful innovation. Whether the source of capital is a local creditor or 
the comparatively wealthy inhabitant of the near-by hacienda, in the ejido 
or share-cropping economy (both exist in the site-cluster area), the 
farmer has little control over the benefits of capital investment. It is 
assumed that capital is scarce in the ejido system and that capital is 
highly complementary with other inputs. 

This model maintains essentially that highly marginal land within 
at least a 7 kilometer radius of the ejido center is due to mutually re¬ 
inforcing political, historical, and economic factors specific to the area. 
These restraints also create a relatively static environment that pre¬ 
vents the ejidatario from either moving to a more fertile land or gaining 
the capital resources to develop a more technologically advanced agri¬ 
cultural environment (i. e. to create an economic surplus to invest in, 
and benefit from, technical innovations, and/or establish a strong eco¬ 
nomic base to influence and possibly restructure the political and eco¬ 
nomic institutions that restrain growth). 

Test Implications: 

1. According to the basic assumption that ejido farms have typically 
small and fragmented plots, the model predicts that within an approximate 
7 kilometer radius of the village , the plots should be highly concentrated 
and have smaller dimensions as opposed to fields near the valley-center. 
There should be a gradient effect when plot size is graphed moving from 
Canada de Alfaro toward the valley-center and Leon. 

2. If the plots are fragmented and located in farm-land of marginal 
utility, then the total land under cultivation should be highly variable 
going across fields radiating from the site-center. 

3. There should be physical characteristics in the environment 
indicating that the farmer in the village lacks an economic surplus to 
counter the economic and political restraints on farm productivity. In 
this case, transportation networks will be used as an indicator of the 
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economic status of the farmer. In the ejido area, one would expect the 
transportation routes to be either trails and rural roads as opposed to 
the more modern transportation networks, such as paved roads, major 
highways, etc. , near the valley-center and Leon. 


Model Two: 

This model hypothesizes that a farmer makes a decision at any one 
location whether to grow cash-crops or subsistence-crops. Acknowledg¬ 
ing Model One's assumption that there should be a high market demand 
for agricultural products in Leon, this model predicts that a farmer will 
want to grow cash-crops as close as possible to the market area. Since 
the village is classified as an ejido, it can be assumed that the village 
has subsistence crop needs. Assuming that there are economic and phys 
ical restraints on the subsistence farmer (i.e. lack of capital to acquire 
trucks or hire labor and a distance factor that inhibits growing crops 
beyond a point attainable by walking from the village), this model pre¬ 
dicts that the farmer will want to grow subsistence crops as close as 
possible to the place of consumption (in this case, the Canada de Alfaro 
area). Because of these restraints placed on the farmer as to the loca¬ 
tion of his plots, the model further establishes that a farmer will grow 
subsistence or cash-crops according to a finite number of hours desig¬ 
nated for transportation to the market. If one assumes that the subsis¬ 
tence farmer grows crops up to 3 1/3 kilometers outside the village 
area, then the distance to solution/market is 2 to 1. At a point 3 1/3 
kilometers outside of Canada de Alfaro, going through space toward 
Leon, it is predictable that subsistence-crops be grown in this area and 
that beyond this point, cash-crops will be grown. For instance, as one 
nears the regional market, Leon, a fewer number of man-hours are 
required to transport cash-crops to the market. Although the cash- 
crop is grown nearer to Leon (the crop destination), the farmer who 
is based in Canada de Alfaro must typically travel the distance by foot 
or burro thus increasing the number of man-hours invested. Graph 1 
below illustrates this relationship between distance to each respective 
market and time alloted for transportation: 

Graph 1 



de Alfaro km 
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From this graph it is apparent that transportation cost (in terms 
of hours) are maximized approximately 3 1/3 kilometers outside of 

Canada de Alfaro and 12/3 kilometers outside of Leon. At this point, 
transportation costs, from the subsistence farmer's point-of-view, are 
prohibitive if the farmer wishes to engage in subsistence-crop agriculture 
since he can only designate a certain ration for transportation (in terms 
of hours). 

Test Implications: 

1. If plot size is used as an indicator of crop-type (subsistence or 
cash) and one assumes that subsistence plots are smaller and fragmented, 
then the size of plot should increase the nearer one comes to the regional 
market, Leon. Specific elements of the model require that plot size in¬ 
creases dramatically approximately 3 1/3 kilometers outside the village 
of Canada de Alfaro. 

2. If cash-crops have larger plots and are closer to the regional 
market, there should also be less variability in the total square miles 
under cultivation in this area. Travelling the reverse course toward 
Canada de Alfaro, there should be a high degree of variability in total 
land under cultivation. 

3. As in Model One, there should be a pattern to the type and spa¬ 
tial orientation of transportation networks. In an area where subsistence 
crops are grown one would expect an abundance of trails and rural roads 
(dirt or gravel). Again, one would expect a more advanced (not implying 
greater utility) system of transportation in regions where cash-crops are 
grown. 

Model Three: 

This model hypothesizes that any variation in the total square miles 
under cultivation in the area should correspond to the proximity of the 
regional market, Leon. From socio-economic and housing characteris¬ 
tics of both population centers it can be inferred that there is a multiple- 
goods market demand for agricultural products grown in the local area 
(see Tables 2 and 3 for comparative socio-economic and housing com¬ 
paratives between the two locations). Since the village of Canada de Alfaro 
is located adjacent to the site-cluster area, this model bases its predic¬ 
tions on an assumption that individual farmers will grow crops as close 
as possible to the market center. In short, from the perspective of a 
farmer, the location of farm plots, and hence, the area of highest agri¬ 
cultural production, is a function of distance from the market area. 

From this assumption, it follows that the total square miles under cul¬ 
tivation should also increase as one nears the perimeter of Leon. An 
important condition to this model requires that agricultural production 
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be concentrated between the two geographical points, Canada de Alfaro 
and Leon, as opposed to an equal distribution of total square miles under 
cultivation in the area covered by the aerial photographs. Figure 2 
below illustrates this relationship between proximity of the regional 
market and the total square miles under cultivation: 

Figure 2 


Canada Increasing Sq. Miles 

_ _ f. . Leon 

de Alfaro —— under Cultivation —— 

Test Implications: 

1. The model predicts that there should be a gradient when total 
land under cultivation is plotted starting in Canada de Alfaro and termi¬ 
nating (at the highest level) near the perimeter of Leon. 

2. The model also implies that the agricultural land near the regional 
market is highly desirable. If this is true, there should be a low degree 
of variability of total land under cultivation around the perimeter of Leon. 

3. Going through space in any direction other than north (or toward 
Leon), the model requires that the total land under cultivation decrease. 
This behavior would support the model's condition that production be 
concentrated between the two geographical points, Canada de Alfaro and 
Leon. 

4. When one travels through space in any direction other than toward 
Leon there should be a greater variability in total land under cultivation. 

It follows that there should be lower figures for land under cultivation as 
opposed to higher figures between Canada de Alfaro and Leon. 


Test Results 

Square Kilometers of Land under Cultivation per Photograph: 

Graph 2 was plotted to show the total square kilometers under cul¬ 
tivation as one progresses north from the base of the series of photo¬ 
graphs.^ - Each line represents one series of photographs. Line 1 shows 
that totals per photograph for series 37, photographs 57 through 46 
(following the south to north progression of photographs on the 'Agricul¬ 
tural Land Map'). Line 2 represents series 36, photographs 6 through 

*^See Table I for the total square kilometers under cultivation per 
photographic unit. 




















- 74 - 


18. Line 3 represents series 34, photographs 9 through 17. 


Graph 2 


Total Sq. 
miles of 
land under 
cultivation 



Two characteristics of the graphs are particularly important to the 
models proposed. The first pattern that develops is that the totals for 
land under cultivation generally increase as one nears the northern 
area of photographs. Second, there is erratic behavior of total land 
under cultivation near the Canada de Alfaro area. Line 2, which repre¬ 
sents the land nearest the ejido and the valley-edge, increases drama¬ 
tically as a result of a total of 5.7 square kilometers under cultivation 
in photograph 36-13. At approximately the same geographical location in 
photograph 34-12 , there is a low total of land under cultivation. In both 
lines, however, an upward trend develops after this decrease and levels 
off at approximately the same total as Line 1. 


Delta Fractions in Land under Cultivation per Photograph: 

To test for degrees of variation in total land under cultivation per 
photograph, a series of delta fractions were calculated. Table 1 shows 
the total square miles under cultivation in each photograph. Using 
each series of photographs as arbitrarily defined units of space, the 
delta fractions for every two consecutive photographs can be calculated. 2 


x- 


where y = the southern-most photograph 


The formula used was 

of the two photographs and' x = the northern photograph. The delta fractions 
can be used as an indicator of variability of land under cultivation in the 
area covered by the photographs. For example, in column 1, photograph 
series 37 starting at the base, delta fractions progress as such: 


Each unit is an equal spatial area for comparison with other units 
in the same series. The delta fraction measures the percentage change 
in the amount of cultivable land between each pair of contiguous photo¬ 
graphs . 
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Photographs 

North 37-54 

37-55 
37-56 

South 37-57 


Land under Cultivation 

2 . 6 square km 
2 . 6 square km 
2.1 square km 
2 . 8 square km 


Delta Fraction 

. 00 

. 24 fraction increase 
-.25 fraction increase 


The complete delta fraction results are as follows: 


Photos 

37-46 

-47 

-48 

-49 

-50 

-51 

-52 

-53 

-54 

-55 

-56 

-57 


Delta Fraction Photos 


-.08 
. 00 
.09 
-.10 
.18 
. 00 
-. 11 
.42 
. 00 
.24 
-.25 


36-18 

-17 

-16 

-15 

-14 

-13 

-12 

-11 

-10 

-9 


Delta Fraction 


. 00 
. 04 
. 04 
-.07 
-.49 
1 . 11 
.50 
. 00 
.64 


Photos Delta Fraction 


34-17 

-16 

-15 

-14 

-13 

-12 

-11 

-10 

-9 


-.11 
-. 03 
. 03 
.17 
. 80 
-.33 
. 00 
.14 


Graphs 3a, 3b, and 3c (found on pp. 79, 80) show these figures 
plotted on a graph. The points plotted are adjusted taking into account 
the lack of spatial continuity between the photographs travelling east to 
west across the base of the map. 


From the results of the ’delta fraction variation, ' it is apparent 
that a pattern develops for total land under cultivation travelling north¬ 
ward from the base of the photographs. In all three series of photographs, 
the fractions appear to oscillate with amplitudes between highs of .42 
(graph 3a), 1.11 (graph 3b), and .80 (graph 3c) and lows of .08, -.09 
(graph 3a), . 04 (graph 3b) and .03, -.03 (graph 3c). The behavior of 
the oscillation also follows a pattern that corresponds to major geograph¬ 
ical locations discussed in the models. Following the spatial course of 
the photographs, particularly high delta fractions can be identified south 
of any point west of the site-cluster, Canada de Alfaro area. However, 
as one nears a point parallel to Canada de Alfaro and progresses north¬ 
ward to Leon, the amplitude of the oscillation damps and continues to 
exhibit low delta fractions until the series of photographs end. 


Plot Size 

To find an approximation of plot size as one nears the valley-center 
(west of Canada de Alfaro) and Leon (north of Canada de Alfaro) a grid 
pattern was super-imposed on the aerial photographs. The photographs 
used for this test were 34-13 through 34-17 and 36-14 through 36-18. 
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Each grid-unit was „ 63 kilometers by . 63 kilometers and contained 
.40 square kilometers of land. Two transects were drawn from the 
Canada de Alfaro area with grid-units adjacent with the transect and 
contiguous with each other. There were a total of 28 grid-units with 
each transect equalling a distance of 5 kilometers. The grid-pattern 
looked as such: 


Transect la lb 



After this grid-pattern was super-imposed on the photographs the 
number of fields per grid-unit were counted. With the number of fields 
known, these figures served as an indicator of field size, recognizing 
that some interference is to be expected from dwellings and roads. 

Graphs 4a, 4b, 4c, and 4d (located on p. 80 ) show the number 
of fields per grid-unit for the series of fields. Geographical points of 
reference have been marked on the x coordinate of each graph for com¬ 
paring the number of fields per grid-unit. Two characteristics of the 
graphs are relevant to the models being tested. First, the decreasing 
number of plots going either west or north of Canada de Alfaro indicates 
that plot size increases steadily as one nears the valley-center or the 
city of Leon. Second, there appears to be a dramatic increase in num¬ 
ber of plots per grid-unit as each transect passes through Canada de 
Alfaro. Although Graph 4d is less dramatic in its steady decline in 
fields per grid-unit, this particular series of units is located approxi¬ 
mately 1. 2 kilometers south of the site-cluster, Canada de Alfaro area 
and initially travels through hills that flank the valley. A third charac¬ 
teristic of field size is apparent when one observes the pattern of field 
dimension and location. Within the first two groups of contiguous fields 
in the grid-pattern (transect 2a and 2b) fields have highly irregular 
shapes (near-circular, triangular and oval). It is worth noting that 
these fields are located near the hills that surround Canada de Alfaro 
on 3 sides. 

Transportation Network Map 

Along the flanks of the hills, as mentioned throughout the article, 
subsistence agriculture predominates. From the mountainous areas 


» 
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near the site-cluster Canada de Alfaro area to the valley floor one sees 
a transition in the pattern of transportation networks. A transportation 
Map 1 which includes geographical formations and population centers, 
illustrates this transition: 



= trails 

----- = rural roads 
i inn . w—rm - modern highway 

ISBfererei a t - airport 

The community of Canada de Alfaro has limited access to Leon, 
the regional market, through the use of narrow rural roads. In the 
specific area of Canada de Alfaro there are many trails concentrated in 
the village area with two trails leading back into the mountain area. 
Toward the center of the valley one sees that at least 4 major highways 
converge at the city of Leon. The highway that runs parallel to the rural 
road eventually turns into a 4-lane highway about 7 kilometers outside 
of Leon. Adjacent to the highway is an airport. 

Conclusion 

Model One attempts to analyse the political, historical, and economic 
restraints placed on the ejido farm collective. The model 1 s prediction 
that plot size and fragmentation are aberrant in regard to other agricul¬ 
tural areas near Leon and the valley-center is supported by test results. 

It is true that field size increases going north or west of Canada de Alfaro 
and that there is a high degree of variability of total land under cultiva¬ 
tion near Canada de Alfaro. However, the field numbers per grid-unit 
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decline rapidly at a point approximately 1 kilometer outside of the village 
in both directions (west and north) of the transects. This information 
suggests that the ejido does not have geographical boundaries that extend 
to the legal maximum of 7 kilometers. Rather, the ejido boundaries 
probably end 1 to 2 kilometers outside of the village. This condition 
is also supported by the results of the transportation network study. 
Although this particular test is not comprehensive in quantitative terms, 
the test supports the notion that there is a form of subsistence farming 
existent in the site-cluster area. The concentration of trails and rural 
roads within the small geographical area of the ejido makes this test 
result highly complementary with the 'delta fraction 1 and 'plot size 1 
results. 

This understanding of the geographical boundaries of the ejido vil¬ 
lage and land also effects the accuracy of Model Two. This model pre¬ 
dicted that transportation costs would prohibit the subsistence farmer 
from growing cash-crops beyond a point of 3 1/3 kilometers outside the 
village. Although test results clearly contradict this prediction, it is also 
apparent that the model is inaccurate for a more political reason. The 
assumption that the farmer has a choice as to where he will grow his 
crops is erroneous. In fact, the farmer cannot legally make his decision 
beyond the geographical boundaries of ejido land as implied by law. Model 
predictions in regard to (1) size of subsistence farm plot, (2) increasing 
plot size as one nears the regional market, (3) variability between defined 
units of cultivated land, and (4) a transition in transportation networks, 
have a great deal of value in understanding the environmental factors 
influencing cash-crop versus subsistence-crop "decision" that a farmer 
must make. However, results do not support the model's specific predic¬ 
tion concerning the subsistence-crop versus cash-crop transportation 
costs (in terms of hours) ratio. 

Model Three is a deceptively complicated hypothesis since it seems 
intuitive that agricultural production and total land under cultivation should 
increase the closer one comes to a large market with diverse needs. This 
model assumes that production is concentrated in a relatively local area 
and that any land outside this area (e. g. in the 2a and 2b transect area) 
must have unusually low totals of cultivated land. In fact, the model's 
assumptions are generally supported by test results as total cultivated 
land variability between units/photographs decreases according to the 
proximity of the market, Leon. When total square miles under cultiva¬ 
tion figures are considered, the model's prediction is not supported. It 
should be mentioned that the decline might be considered negligible due 
to minor interferences from dwellings and roads that intersperse the 
periphery of the city-limits. In any case, generally high totals for land 
under cultivation prevail in the perimeter area compared to points father 
south in the valley, particularly along the flanks of the hills. The model's 
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major deficiency (and a deficiency associated with Model Two) is that 
it assumes a clear land-use relationship between the regional market, 
Leon and Canada de Alfaro. The model does not incorporate into its 
assumptions a condition that other factors might influence the location 
and size of fields. However, the model’s theoretical orientation does 
have gravity. According to the proximity of the valley-center from any 
point west of Canada de Alfaro, plot size increases, there are higher 
figures for land under cultivation, and low delta fractions between units 
of land. 

Delta Fractions Graphs 









- 80 - 


Delta Fractions Graph 3 c 



Numbers of Plots per Grid 


Graph 4a Transect la 



Graph 4b Transect lb 



Graph 4c Transect 2a 


Graph 4d Transect 2b 
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Each unit within a series 
is an equal spatial area 
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The aerial photographs were supplied by 
Cetenal, a federal agency of Mexico. 
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APPENDIX 


!• Table 1 shows the square kilometers of land under cultivation 
on each photograph when they are arranged according to the order on 
the Agricultural Land Map. 1 


Photo 

Cult. Sq.km 

Table 

Photo 

1 

Cult. Sq.km. 

Photo 

Cult. Sq 

37-46 

3.6 

36-18 

2. 9 

34-17 

3.1 

37-47 

3.9 

36-17 

2. 9 

34-16 

3.5 

37-48 

3. 9 

36-16 

2. 8 

34-15 

3. 6 

37-49 

4. 3 

36-15 

2. 7 

34-14 

3.5 

37-50 

3. 9 

36-14 

2. 9 

34-13 

3. 0 

37-51 

3. 3 

36-13 

5.7 

34-12 

1.6 

37-52 

3.3 

36-12 

2. 7 

34-11 

2.4 

37-53 

3. 7 

36-11 

1. 8 

34-10 

2.4 

37-54 

2. 6 

36-10 

1. 8 

34-9 

2. 1 

37-55 

2. 6 

36-9 

1. 1 



37-56 

2. 1 





37-57 

2. 8 





2. The 

margin of error 

on any figure in Table 

1 is no 

greater the 


2 1/2 per cent, or in some cases less. The two factors that contribute 
to the margin of error are: 

a. the difficulty associated with fitting photographs together and 
measuring them for land under cultivation, and 

b. the minor scale differences in some photographs, related to 
a variable camera angle and airplane altitude. 

3. Rogers (1968) points out that the degree to which farmers are 
efficiency-minded, market oriented and economically ,r rational M (quotes 
added), depends in a large part on the degree of modernization of the 
village. Tables 2 and 3 below show the socio-economic and housing 
characteristics of the city of Leon and the village of Canada de Alfaro, 
located approximately 5. 7 km southeast of Leon and next to the site- 
center. This information is relevant to the models presented in this 
article as the model's assumptions are contingent upon significant popu¬ 
lation, economic, and agricultural production differentials between 
Canada de Alfaro and Leon. 


•j* 

-r* 

The study to measure the totals of land under cultivation was com¬ 
pleted by James Pendergrass in May of 1973. The figures was arrived 
at by using a compensating polar planimeter, a device for accurately 
measuring irregular surface areas. 
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Table 2 


Comparative Socio-economic Characteristics of the 
Population (expressed in percentages): 


# of People 

Econ Active 

Agric 

Indust 

Commerce 

Literacy 

Alfaro 450 

28. 0 

63.5 

23.0 

5. 6 

22.3 

Leon 364,990 

26. 9 

4. 9 

53.9 

34.4 

79. 0 


Table 3 

Comparative Housing Characteristics of the Population 
(expressed in percentages): 

Outside Inside Non-Dirt 


# of Houses 

Private 

water 

water 

Pipes 

floor s 

Elect'y 

Radio 

T. V. 

Alfaro 69 

18. 8 

85. 5 

5.8 

2. 9 

5.8 

. . . 

72. 5 

• • • 

Leon 54,625 

55. 7 

66.2 

9.1 

68 . 1 

85.9 

85. 8 

85.5 

45.4 


Note: The information in Tables 2 and 3 was found in: 
IX Censo General De Poblacion 1970, volumen 1; 
A Guascalientes A Guerrero, 1973. 
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I. Introduction 

A portable refraction seismograph was tested during an archaeo¬ 
logical field expedition in Guanajuato, Mexico in the summer of 1973. 
This paper attempts to analyze the seismographic data and to determine 
the seismograph's archaeological utility. 

The paper has three sections. The first section deals with the 
theory and use of seismic instruments. The next section explains the 
methodology of the fieldwork and the method of analyzing the data. 

The seismograph was used on several sites. On each site linear 
courses were run, in which data was collected relevant to the sub¬ 
surface near that course. The raw data consists of a list of seismo¬ 
graph readings (time and distance) for each course. A series of 
hypothetical underground structures were developed, accompanied by 
hypothetical time-distance graphs of what the seismographic readings 
of these cases would look like. 

From the hypothetical cases and the real time-distance graphs, 
interpretations of the actual profile of each site were made. Conclu¬ 
sions were then drawn as to the validity of the use of the portable re¬ 
fraction seismograph for archaeological purposes. 


II• The Theory of Seismic Explorati on 

Basic seismic principles 

Seismic methods are based on the fact that "the velocity of propa¬ 
gation of a wave or impulse in an elastic body is a function of the modulus 
of elasticity, the Poisson ratio, and the density of the material, and that 
very great differences exist between wave velocities in solid rock and 
in loose, sedimentary deposits. An impulse travelling through the 
earth will produce longitudinal or compression waves, transverse or 
shear waves, and various types of surface waves. In the seismic 
methods of exploration, it is usually the velocity of the longitudinal 
or compression waves which is important. 

Rocks and many accumulations of unconsolidated sand, gravel, and 
clay are elastic, and resist an applied stress by a force that is approxi¬ 
mately proportional to displacement. Therefore, when disturbed, they 
oscillate with simple harmonic motion. When one portion of such a 
medium is displaced, adjacent particles are displaced. As time passes, 
more and more particles are included in the displaced portion. The 
result is the advance of an elastic wave through the medium. 


M.J. Hvorslev, Subsurface Exploration and Sampling of Soils for 
Civil Engineering Purposes . p. 29. 
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The longitudinal or compression wave 

The longitudinal wave is a body wave, that is, it travels through 
the interior of an elastic medium. It advances by alternating compres¬ 
sion and rarefactions, and thus is also known as a compressional wave. 

It is called longitudinal because the sense of motion of a particle is 
along the direction of the wave’s advance. The symbol P denotes this 
wave type. 

Wave propagation 

If a material in which P waves travel with a certain velocity, 
adjoins another material in which they travel’with a different velocity, 
a wave encountering the boundary between the materials will produce 
displacements at that boundary which result in the setting up of (1) a 
reflected P-wave, and (2) a refracted P-wave advancing into the new 
material. The distribution of the initial energy among the derived waves 
is determined by the angle at which the wave approaches the boundary 
and the contrast in elastic properties of the two materials. ^ 

The refraction principle 

The refraction method is based on the fact that a wave in travelling 
across a boundary is refracted toward the plane of the boundary when it 
enters a material which transmits the wave with a higher velocity and 
toward the perpendicular to the boundary when it enters a material which 
transmits the wave with a lower velocity. ^ 

The example below will illustrate how the seismic method operates: 
There is a wave front in a 2 layer structure consisting of about 80 feet 
of gravel on bedrock. The wave travels about 5, 000 feet per second, 
in the gravel and about 20, 000 ft/sec. in the bedrock. In such a condi¬ 
tion a wave front will be an approximate sphere of 5 0 feet radius . 01 
seconds after it leaves the shot point. The wave is created by a blow or 
a ’’shot” at a particular place and time upon the surface of the earth. 
About „ 016 seconds after the shot, it reaches the rock surface. From 
this point of contact, a wavelet starts into the rock with a velocity of 
20, 000 ft/sec. At . 02 seconds after the shot the wavelet has reached 
the distance noted, not only vertically but radiating in all directions 
within the rock, so that even along the gravel-rock contact the wave 
front has reached the point B (see Figure 1). However, along the 
boundary the wave in the rock starts a new series of wavelets that 
emanate into the gravel and are directed to the gravel surface. 

- 

D.L. Leet, Earth Waves, p. 59. 

3 

M.J. Horslev, Subsurface Exploration and Sampling of Soils for 
Civil Engineering Purposes, p.30. 
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Figure 1 



The wave front travels through the gravel on to each instrument 
at 5, 000 ft/sec, but there is a point along the line, or profile, where 
the wave front refracted from the rock, owing to its higher velocity 
gained from the rock, will reach an instrument before the direct wave. 
In the diagram this is a point a little more than 200 ft. from the shot 
point. Thus the velocity of the materials may be determined; those 
of the earlier arrivals at the geophone may be computed to be 5, 000 
ft/sec, the velocity gravel in the problem, and the later arrivals 
may be computed to be 20, 000 ft/sec, or the velocity of the bedrock. 
The distance from the shot point to the point of intersection of these 
two lines is controlled ky the thickness of the gravel, and the greater 
the thickness, the greater will be the distance of the intersection from 
the shot point. ^ 


P- 


4 


78 


Daniel Linehan, 1 Seismology applied to shallow zinc 
(Figures 1 and a ^ ter Lenehan’s diagram). 


research. 11 









Diagram of wave fronts, and the relevant travel-time graph 
for a 400-foot refraction profile 
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III. Basic Soil Mechanics and the Seismic Principle 

The use of the seismic method to explore subsurface structures is 
based on the fact that the velocity of wave propagation differs a great 
deal in different materials. 

Velocity characteristics of different materials 

The age of the material through which a given wave propagates is 
the critical variable determining the wave's velocity. In general, it can 
be stated that the velocity in unconsolidated materials is low and in con¬ 
solidated materials it is high. The velocity in new or re-deposited soils 
is usually much lower than in undisturbed soils. Often, unusually low 
velocities are recorded in deep fills, especially if they contain refuse, 
ashes, cinders, and other similar materials. Also very low velocities 
have been recorded in a mass of loose rock and boulders which contained 
little or no earth, but in which there was no continuity of structure. 

A table has been constructed showing the speed of wave propagation 
in different materials, and shows in a general way the approximate 
range of velocities which may be expected (see Table 1). (For a more 
complete listing, see Shepard, 1939: Table 1.) 


Table 1 


Velocity Table 


600 - 1, 200 ft/sec. 

1, 200 - 1, 900 ft/sec. 

1, 600 - 2, 600 ft/sec. 

1, 800 - 7, 000 ft/sec. 
1,300 - 2,500 ft/sec. 
4,600 - 6, 000 ft/sec. 

2, 600 - 11, 000 ft/sec. 

3, 000 - 9, 000 ft/sec. 
10, 000 - 20, 000 ft/sec. 

6, 000 - 20, 000 ft/sec. 


Sand and Topsoil 
Sandy Clay 
Gravel 
Glacial Till 
Rock Talus 

Water in Loose Materials 

Shale 

Sandstone 

Granite 

Limestone 


When a layer of homogeneous soil overlaps one of higher elastic 
properties, there will be a critical shooting distance for which the time 
of travel through the upper layer is just equal to the time of travel over 
the longer path which penetrates the deeper and faster layer. At this 
critical distance there will be a break in the time-distance graph to a 
different slope, which represents the velocity of wave propagation in 
the second layer. Similarly, when the critical shooting distance for 
the third layer is reached there will be another break in the graph to 
a new slope, which represents the velocity of the third stratum. 
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The above assumes that the path of the seismic wave is vertical 
from the shot point (or hammer point) to the rock, along the rock to a 
point directly under the detector, and thus vertical to the detector. 

This assumption is only useful when dealing with shallow depth deter¬ 
minations o 

When the surface of the rock is parallel to the ground surface, 
the slope of the graph (see Figure 2) will give true velocities• 
However, when the interface is not parallel to the ground surface, the 
profile line must be shot from opposite directions to determine the 
true velocity. 

If the rock surface slopes downward from the shot point, the re¬ 
fracted wave has to travel farther to the surface at each successive 
reading, and the apparent velocity will be slower than the true velocity. 
If the rock slopes upward from the shot point, there is less distance 
for the wave to travel and the apparent velocity will be greater than the 
true velocity. The apparent velocities do indicate the direction and the 
approximate slope of a rock surface. 


IV. The Refraction Seismograph 

Principle of the refraction seismograph 

The following is a description of the basic one-channel seismograph. 
This instrument consists of a timing unit, an amplifier, and a display 
unit, all built into one box. It has a separate power unit, geophone, 
and a 10 pound sledge hammer. The hammer is connected to a switch, 
and when the hammer hits the ground, the switch closes and a pulse is 
fed into the timing unit. After a certain time delay, which is controlled by 
the operator, a marker appears on a cathode-ray tube screen. The 
seismic wave is picked up by the geophone, and appears on the same 
screen. The timing dial is then rotated so that the beginning of the 
first impulse coincides with the marker. The time delay, which is 
equal to the travel time of the fastest wave, is then read on the timing 
dial. 


In the field, the geophone is placed at the location where the depth 
to bedrock is required, and the hammer blows are along a profile at 
spaced intervals. 

5 

The terra-Scout model R-15Q 

The Terra-Scout essentially follows this basic design. It is made 
up of three parts: the hammer, geophone, and the instrument which 


The portable refraction seismograph used in this fieldwork was 
a Terra-Scout model R-150, manufactured by Soiltest Inc. of Evanston, 
Illinois. 
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records the wave from the geophone on a cathode-ray tube screen, and 
measures the travel time of the wave between the hammer and the geo¬ 
phone. 

The instrument itself consists of various controls: the off-gain, 
delay, range, test button, intensity, focus, and horizontal and vertical 
positioning. The gain control increases or decreases the size of the 
wave on the screen. The delay read-out varies the interval of the time 
between the hammer blow and the start of the sweep across the screen. 
The range provides 2 ranges by acting as a multiplier of the delay dial 
readings. Therefore, this timing read-out can be direct or by the multi¬ 
plier in milliseconds. It is accurate to -I- 0.25 milliseconds from 0 to 
25 0 milliseconds. The test button acts as the hammer switch to produce 
a sweep across the screen, and is used as a check. The wave display 
is on a 5-inch cathode-ray tube controlled by horizontal and vertical 
centering, intensity and focus controls. The tamper, or hammer, has 
a triggering accuracy of jf5 microseconds and is damped with nitrogen 
and helium gases. This instrument commonly gives results to depths 
of 75 feet to 100 feet, and in some areas can be used up to 200 feet 
in depth. 

Problems with and limitations of the refraction seismograph 

When using the seismic method of investigation, there are numerous 
problems and limitations of the instrument itself that must be taken into 
account in the interpretation and analysis. 

The seismic refraction method is based on the assumption that the 
velocity of the seismic waves will increase with increasing depth. A 
layer which transmits the waves at a lower velocity than the overlying 
layer cannot be detected or its thickness determined. (In fact, a layer 
must have a certain minimum thickness in order to be detected at all. ) 

If the masked layer is underlain by a layer of higher velocity than the 
overlying layer, the resulting graph will be a 2-layer curve with the 
middle layer absent, but the calculated depth to the lowest layer will 
be so much greater than the actual depth that any familiarity with the 
geology of the area will lead to doubts as to the validity of the results. 
When working in a given area, some of the graphs might show a second 
layer while on other graphs the layer will be absent or only slightly 
indicated. If the presence of a second stratum is indicated by some 
graphs, it is probable that its extent v/ould be general, and it would 
be logical to include the second phase wherever there is any indication 
of it. This problem has been dealt with in more detail by Kaila and 
Narain (1970). 

Advantages and disadvantages of the refraction method 


The refraction method of investigation is most powerful for re¬ 
connaissance in a new area. Records of the first pulse to arrive at 
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the detector can be made for a wide range of source-detector distances. 
Figure 3 shows the various paths a refracted wave may take: 

Figure 3 



In addition to the direct wave AD, a refracted wave can travel by the 
minimum time path ABCD. It is necessary for the lower layer BC to 
transmit compressional waves faster than the upper layer. At long 
distances the refracted wave overtakes the direct wave, and can be 
clearly read because it is the first wave recorded. 

The greatest strength of the refraction method is that it allows a 
picture of the strata to be drawn and provides a velocity for each layer. 

One of the limitations placed on this method is that the length of 
the refraction spread must be many times longer than the depth to the 
deepest refractor, usually 5 or 6 times. 

In order for the refracted wave to return from a layer to the sur¬ 
face, the velocities must increase in each subjacent layer. This in¬ 
capacitates this method in areas where the upper layer of the ground 
freezes and masks arrivals to a distance many times the thickness of 
the frozen layer. In other areas, layers of rock are usually not laid 
down in order of increasing velocity downwards, giving velocity re¬ 
versals. A hard thin layer can be sandwiched between low velocity 
layers, such as clay, and the refracted wave from the thin layer will 
mask that from the lower velocity layer underneath. 

The author is aware of other relevant problems such as "surface 
interference, " "wave dispersion, " "groundroll, " "reflection problems, " 
"surface regularity, " and "weathering. " However, space precludes 
their discussion here. For a complete discussion see Skinner, 1974. 


V. Methodology of Seismograph Fieldwork 

The decision of which site to seismograph was based on various 
factors. Seismographing was done on sites previously excavated and 
on sites where cores were taken by the hydraulic drill. This method 
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was used in order to supply a data base for correlation. Where excava¬ 
tion and/or drilling was not done, seismographing was done on sites 
where a fair amount of archaeological information on subsurface struc¬ 
tures was known. A transect for a course on a particular site was 
somewhat arbitrarily chosen. 

Using the instruction manual for the Terra-Scout model R-15 0 as a 
guide, the spacing of the hammer blows was set at intervals of 10 feet. 
The course was usually run for 100 feet, but where it seemed necessary 
to obtain continuous information over a long distance, the course was 
run to 180 feet — the length of the hammer wire. 

In dealing with homogeneous layers, the course length is usually 
determined by the depth required. The depth tested is about 1/4 times 
the course length. A high velocity material near the surface would 
normally decrease the required length because one would only need 
enough points to show the slope of this high velocity material. However, 
when the intent is to follow a subsurface structure, the length of the 
course will necessarily be arbitrary, because the length of the struc¬ 
ture is unknown. 

In normal homogeneous layers, a reversed course is run to check 
for sloping subsurface layers, and to get an accurate depth if there are 
sloping layers. A reverse course is run by moving the geophone to the 
last hammering point of the forward course and the course is run back 
toward the point where the first course geophone was placed. Each 
course by itself gives only an apparent velocity. By averaging the re¬ 
versed course with the forward course, the true velocity is determined. 

However, this was not clearly understood by the operators, and a 
reversed course was run where the hammering was done in a reversed 
direction, but without moving the geophone. This simply checks the 
consistency of the operator without adding new information. 

Once the reading on the oscilloscope begins, the size of the wave 
must be adjusted. This is done with the GAIN control and must be 
adjusted when the hammering begins at each new distance from the 
geophone. The first wave to reach the geophone is the important one 
for timing purposes. Often near the critical distance there is a change 
in the wave shape and a decrease in wave amplitude. The GAIN control 
must be increased enough to recognize the smaller wave and to use it 
fof- timing. 

A marker is used as a reference point in the determination of the 
time delay between the hammer blow and the wave arrival at the geo¬ 
phone. This marker is placed on the bottom of the first wave at each 
blow by turning the DELAY dial clockwise while hammering until it is 
in the bottom of the first wave. From this, the time is read off of 
the DELAY dial and recorded. 
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VI. Method of Analysis 

All readings were converted into milliseconds. A computer program 
was then developed to plot a time-distance graph, with distance on the 
horizontal axis and time on the vertical axis. The scale of this graph 
varied, depending on the length of the course and the range of milli¬ 
second readings. 

All possible archaeological cases for that area were hypothesized, 
and graphs of those cases were drawn. 

The points on the time-distance graphs were connected, keeping in 
mind at all times the slope of the ground surface under the transect. 
Interpretations were made, based in a large part on knowledge of the 
sites, such as excavation data, core samples from drill borings, re¬ 
construction drawings, and general knowledge of the most likely struc- 
tures to occur. 

The critical distances - where the layers change - were determined, 
and the velocities of each layer calculated. Equations for depth calcu¬ 
lations, height of walls, and width of depressions were developed, and 
actual values calculated. The depth equation only applies as long as 
the velocity of each successive layer increases. Correlations were 
attempted where possible with excavation data and drawings. Last of 
all, a profile was drawn of the interpretation of each case. 
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VII. Data 

Lists of Seismograph Readings 

Site - Rio Guanajuato, Course 1, Range 0-500 Reverse Course 


Seismograph Seismograph 



readings 

Milliseconds 


readings 

Milliseconds 

10 ft. 

36 

18. 0 

100 ft. 

263 

131.5 

20 ft. 

102 

51. 0 

90 ft. 

244 

122. 0 

30 ft. 

no 

55. 0 

80 ft. 

222 

111. 0 

40 ft. 

142 

71.0 

70 ft. 

200 

100. 0 

50 ft. 

170 

85. 0 

60 ft. 

181 

90.5 

60 ft. 

183 

91.5 

5 0 ft. 

170 

85, 0 

70 ft. 

203 

101.5 

40 ft. 

144 

72. 0 

80 ft. 

225 

112.5 

3 0 ft. 

120 

60. 0 

90 ft. 

243 

121.5 

20 ft. 

101 

50.5 

100 ft, 

268 

134. 0 

10 ft. 

32 

16. 0 

Site - 

Rio Guanajuato, 

Course 2, Range 

0-500 

Reverse Course 


Seismographs 


Seismograph 



readings 

Milliseconds 


readings 

Milliseconds 

10 ft. 

125 

62.5 

100 ft. 

387 

193.5 

20 ft. 

170 

85. 0 

90 ft. 

380 

190. 0 

3 0 ft. 

272 

136.0 

80 ft. 

357 

178.5 

40 ft. 

308 

154. 0 

70 ft. 

333 

166.5 

50 ft. 

333 

166.5 

60 ft. 

340 

170. 0 

60 ft. 

342 

171. 0 

50 ft. 

325 

162.5 

70 ft. 

333 

166.5 

40 ft. 

310 

155.0 

80 ft. 

351 

175.5 

30 ft. 

268 

134. 0 

90 ft. 

378 

189. 0 

20 ft. 

15 0 

75.0 

100 ft 

3 91 

195.5 

10 ft. 

128 

64. 0 

Site - 

Rio Guanajuato, 

Course 3, Range 

5 0-500 

Rever: 

se Course 


Seismogrphs 



Seismograph 



readings 

Milliseconds 


readings 

Milliseconds 

10 ft. 

40 

20. 0 

100 ft. 

500 

250. 0 

20 ft. 

89 

44. 5 

90 ft. 

456 

228. 0 

3 0 ft. 

109 

54. 5 

80 ft. 

375 

187.5 

40 ft. 

148 

74. 0 

70 ft. 

348 

174. 0 

5 0 ft. 

202 

101.0 

60 ft. 

308 

154. 0 

60 ft. 

276 

138. 0 

50 ft. 

201 

100.5 

70 ft. 

345 

172.5 

40 ft. 

144 

72. 0 

80 ft. 

374 

187.0 

3 0 ft. 

103 

51. 5 

90 ft. 

452 

226. 0 

20 ft. 

90 

45.0 

100 ft 

502 

251. 0 

10 ft. 

53 

26.5 
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Site - Rio Guanajuato, Course 4, Range 0-500 Reverse Course 



Seismograph 

readings 

Milliseconds 


Seismograph 

readings Milliseconds 

10 ft. 

37 

18.5 

170 ft. 

241 

120.5 

20 ft. 

61 

30. 5 

160 ft. 

235 

117 .5 

3 0 ft. 

80 

40. 0 

15 0 ft. 

233 

116 .5 

40 ft. 

96 

48. 0 

140 ft. 

215 

107.5 

5 0 ft. 

114 

57. 0 

13 0 ft. 

209 

104.5 

60 ft. 

133 

66.5 

120 ft. 

207 

103.5 

70 ft. 

147 

73.5 

110 ft. 

199 

99.5 

80 ft. 

154 

77. 0 

100 ft. 

188 

94. 0 

90 ft. 

167 

83.5 

90 ft. 

169 

84.5 

100 ft. 

192 

96. 0 

80 ft. 

154 

77. 0 

110 ft. 

196 

98. 0 

70 ft. 

146 

73. 0 

120 ft. 

210 

105.0 

60 ft. 

134 

67. 0 

13 0 ft. 

211 

105.5 

5 0 ft. 

113 

56.5 

140 ft. 

216 

108. 0 

40 ft. 

93 

41.5 

15 0 ft. 

233 

116.5 

30 ft. 

78 

39. 0 

160 ft. 

244 

122. 0 

20 ft. 

63 

3 1.5 

170 ft. 

247 

123.5 

10 ft. 

36 

18. 0 


Site - 

Canada de Alfaro, 

Course 1 


Reverse 

Course 


Site 43, Range 0 

-100 





Seismograph 



Seismograph 



readings 

Milliseconds 


readings 

Milliseconds 

10 ft. 

167 

16.7 

100 ft. 

423 

42.3 

20 ft. 

308 

30. 8 

90 ft. 

419 

41. 9 

3 0 ft. 

340 

34. 0 

80 ft. 

412 

41.2 

40 ft. 

359 

35.9 




50 ft. 

391 

39.1 

60 ft. 

405.5 

40. 6 

60 ft. 

409 

40.9 

50 ft. 

389 

38.9 

70 ft. 

412 

41. 2 

40 ft. 

364 

36.4 

80 ft. 

413.5 

41.4 

30 ft. 

343 

34.3 

90 ft. 

418 

41. 8 

20 ft. 

307 

30. 7 

100 ft 

418.5 

41.9 

10 ft. 

164 

16.4 
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Canada de Alfaro, Site 43, Course 2 , Range 0-100 Reverse Course 



Seismograph 

readings 

Milliseconds 


Seismograph 

readings 

Milliseconds 




0-500 Range 


10 ft. 

170 

17. 0 

100 ft. 

185 

92.5 

20 ft. 

225.5 

22.6 

90 ft. 

181 

90.5 

3 0 ft. 

340.5 

34.1 

80 ft. 

179.5 

89.5 

40 ft. 

386.5 

38.6 

70 ft. 

181 

90.5 

50 ft. 

400 

40. 0 

60 ft. 

165.5 

82.5 

60 ft. 

421 

42.1 

0-100 Range 


70 ft. 

43 0 

43.0 

50 ft. 

389 

38. 9 

0-500 Range 


40 ft. 

382 

38.2 

80 ft. 

180 

90. 0 

3 0 ft. 

337 

33.7 

90 ft. 

187 

93.5 

20 ft. 

330 

33.0 

100 ft. 

185.5 

92.5 

10 ft. 

319 

31. 9 


Canada de Alfaro, Site 43, Course 3, Range 0-100 Reverse Course 


Seismograph 



readings 

Milliseconds 

10 ft. 

89 

8. 9 

20 ft. 

181 

18. 1 

3 0 ft. 

267 

26. 7 

40 ft. 

294 

29.4 

5 0 ft. 

306 

30. 6 

60 ft. 

314 

31.4 

70 ft. 

325 

32. 5 

80 ft. 

331 

33. 1 

90 ft. 

342 

34. 2 



Seismograph 

readings 

Milliseconds 

100 ft. 

386 

38.6 

90 ft. 

343 

34. 3 

80 ft. 

332 

33.2 

70 ft. 

326 

32. 6 

60 ft. 

315 

31.5 

5 0 ft. 

308 

30. 8 

40 ft. 

293 

29.3 

3 0 ft. 

267 

26. 7 

20 ft. 

180 

18.0 
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Canada de Alfaro, Site 22, Course 1, 


Seismograph 

readings Milliseconds 


10ft. 

18 

9. 0 

20 ft. 

36 

18. 0 

3 0 ft. 

56 

28. 0 

40 ft. 

67 

33.5 

5 0 ft. 

77 

38.5 

60 ft. 

80 

40. 0 

70 ft. 

82 

41. 0 

80 ft. 

86 

43. 0 

90 ft. 

91 

45.5 

100 ft. 

91 

45.5 

110 ft. 

102 

51. 0 

120 ft. 

104 

52. 0 

130 ft. 

107 

53.5 

140 ft. 

111 

55.5 

15 0 ft. 

118 

59. 0 

160 ft. 

122 

61. 0 

170 ft. 

13 0 

65.0 

180 ft. 

128 

64. 0 


Range 0-500 Reverse Course 



Seismograph 

Millis ec 


readings 

180 ft. 

128 

64. 0 

170 ft. 

128 

64. 0 

160 ft. 

122 

61.0 

150 ft. 

117 

58.5 

140 ft. 

no 

55. 0 

13 0 ft. 

106 

53. 0 

120 ft. 

104 

52. 0 

110 ft. 

102 

51. 0 

100 ft. 

90 

45.0 

90 ft. 

90 

45.0 

80 ft. 

84 

42. 0 

70 ft. 

82 

41. 0 

60 ft. 

79 

39.5 

5 0 ft. 

77 

33.5 

40 ft. 

68 

34. 0 

30 ft. 

54 

27. 0 

20 ft. 

37 

18.5 

10 ft. 

29 

14.5 


Canada 

de Alfaro, 

Site 22, Course 


Seismograph 



readings 

Milliseconds 

10 ft. 

18 

90. 0 

20 ft. 

37 

18.5 

3 0 ft. 

53 

26.5 

40 ft. 

69 

34.5 

50 ft. 

94 

47.0 

60 ft. 

no 

55. 0 

70 ft. 

132 

66. 0 

80 ft. 

145 

72.5 

90 ft. 

152 

76.0 

100 ft. 

177 

88.5 

110 ft. 

216 

108. 0 

120 ft. 

229 

114.5 

130 ft. 

244 

122. 0 


Range 0-500 Reverse Course 

Seismograph 



readings 

Milliseconds 

13 0 ft. 

244 

122. 0 

120 ft. 

229 

114.5 

110 ft. 

218 

109. 0 

100 ft. 

180 

90. 0 

90 ft. 

168 

84. 0 

80 ft. 

143 

71.5 

70 ft. 

132 

66.0 

60 ft. 

111 

55.5 

50 ft. 

94 

47. 0 

40 ft. 

67 

33.5 

30 ft. 

65 

32.5 

20 ft. 

37 

18.5 

10 ft. 

18 

9.0 
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Canada de Alfaro, Site 22, Course 3 


Seismograph 

readings Milliseconds 


10 ft. 

15 

. 75 

20 ft. 

36 

18. 0 

3 0 ft. 

53 

26.5 

40 ft. 

69 

34.5 

50 ft. 

73 

36.5 

60 ft. 

81 

40.5 

70 ft. 

88 

44. 0 

80 ft. 

95 

47.5 

90 ft. 

101 

50. 5 

100 ft. 



110 ft. 

127 

63.5 

120 ft. 

135 

67.5 

13 0 ft. 

157 

78.5 

140 ft. 

162 

81. 0 

15 0 ft. 

176 

83.0 

160 ft. 

186 

93. 0 

170 ft. 

191 

95.5 

180 ft. 

214 

107. 0 


Canada de Alfaro, Site 22, Course 

Seismograph 

readings Milliseconds 

10 ft. 

16 

8. 0 

20 ft. 

48 

24. 0 

3 0 ft. 

58 

29. 0 

40 ft. 

65 

32. 5 

50 ft. 

75 

37.5 

60 ft. 

83 

41. 5 

70 ft. 

89 

44.5 

80 ft. 

96 

48. 0 

90 ft. 

104 

52. 0 

100 ft. 

106 

53. 0 


Range 0-500 Reverse Course 


Seismograph 

readings Milliseconds 


180 ft. 

214 

107..0 

170 ft. 

192 

96. 0 

160 ft. 

185 

92.5 

15 0 ft. 

177 

88.5 

140 ft. 

15 9 

79.5 

13 0 ft. 

155 

77.5 

120 ft. 

133 

66.5 

110 ft. 

124 

62. 0 

100 ft. 

104 

52. 0 

90 ft. 

99 

49.5 

80 ft. 

92 

46. 0 

70 ft. 

83 

41.5 

60 ft. 

77 

38. 5 

50 ft. 

69 

34.5 

40 ft. 

64 

32. 0 

30 ft. 

60 

30. 0 

20 ft. 

31 

15.5 

10 ft. 

15 

7.5 


Range 

0-500 Reve: 

rse Course 


Seismograph 



readings 

Milliseconds 

100 ft. 

106 

53. 0 

90 ft. 

104 

52. 0 

80 ft. 

98 

49. 0 

70 ft. 

89 

44.5 

60 ft. 

85 

42.5 

50 ft. 

76 

38. 0 

40 ft. 

67 

33.5 

3 0 ft. 

59 

28.5 

20 ft. 

43 

21. 5 

10 ft. 

16 

8. 0 
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Canada de Alfaro, Site 11, Course 1, 

Seismograph 



readings 

Milliseconds 

10 ft. 

14 

7. 0 

20 ft. 

27 

13.5 

30 ft. 

36 

18. 0 

40 ft. 

45 

22.5 

5 0 ft. 

50 

25. 0 

60 ft. 

56 

28. 0 

70 ft. 

87 

43.5 

80 ft. 

94 

47. 0 

90 ft. 

92 

46. 0 

100 ft. 

98 

49. 0 

Canada de Alfaro, Site 11, Course 2, 


Seismograph 



readings 

Milliseconds 

10 ft. 

14 

7. 0 

20 ft. 

30 

15 . 0 

3 0 ft. 

40 

20. 0 

40 ft. 

45 

22. 5 

5 0 ft. 

53 

26.5 

60 ft. 

57 

28.5 

70 ft. 

61 

30. 5 

80 ft. 

65 

32.5 

90 ft. 

67 

33.5 

100 ft. 

71 

35.5 

Canada de Alfaro, Site 10, Course 1 : 


Seismograph 



readings 

Milliseconds 

10 ft. 

7 

3.5 

20 ft. 

24 

12. 0 

3 0 ft. 

47 

23.5 

40 ft. 

60 

30. 0 

50 ft. 

68 

34. 0 

60 ft. 

69 

34.5 

70 ft. 

70 

35.0 

80 ft. 

71 

35.5 

90 ft. 

72 

36. 0 

100 ft 

74 

37. 0 


Range 0-500 Reverse Course 

Seismograph 


readings Milliseconds 


100 ft. 

97 

48.5 

90 ft. 

91 

45.5 

80 ft. 

93 

46.5 

70 ft. 

88 

44. 0 

60 ft. 

55 

27.5 

50 ft. 

50 

25. 0 

40 ft. 

45 

22. 5 

3 0 ft. 

37 

18.5 

20 ft. 

28 

14. 0 

10 ft. 

14 

7. 0 

Range 0-500 

Reverse Course 

Seismograph 


readings 

Milliseconds 

100 ft. 

70 

35. 0 

90 ft. 

66 

33.0 

80 ft. 

63 

31.5 

70 ft. 

60 

30. 0 

60 ft. 

55 

27.5 

5 0 ft. 

52 

26. 0 

40 ft. 

45 

22.5 

3 0 ft. 

39 

19.5 

20 ft. 

31 

15.5 

10 ft. 

18 

9. 0 

Range 0-500 

Reverse Course 

Seismograph 


readings 

Milliseconds 

100 ft. 

75 

37.5 

90 ft. 

73 

36.5 

80 ft. 

71 

35.5 

70 ft. 

70 

35. 0 

60 ft. 

69 

34.5 

5 0 ft. 

68 

34. 0 

40 ft. 

60 

30. 0 

3 0 ft. 

47 

23.5 

20 ft. 

24 

12. 0 

10 ft. 

4 

2. 0 
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Canada de Alfaro, Site 10, Course 2 



Seismograph 

readings 

Milliseconds 

10 ft. 

2 

1 . 0 

20 ft. 

26 

13. 0 

30 ft. 

51 

25.5 

40 ft. 

60 

30. 0 

50 ft. 

63 

31.5 

60 ft. 

66 

33. 0 

70 ft. 

67 

33.5 

80 ft. 

71 

35.5 

90 ft. 

71 

35.5 

100 ft. 

72 

36. 0 

Canada de Alfaro, Site 57, Course 
Seismograph 


readings 

Milliseconds 

10 ft. 

8 

4. 0 

20 ft. 

21 

10.5 

30 ft. 

42 

21.0 

40 ft. 

56 

28.5 

5 0 ft. 

61 

30.5 

60 ft. 

63 

31.5 

70 ft. 

64 

32. 0 

80 ft. 

70 

35. 0 

90 ft. 

73 

36.5 

100 ft. 

73 

36.5 

Canada de A.lfa.i*o 3 Site 57, Course 

Seismograph 


readings 

Milliseconds 

10 ft. 

7 

3.5 

20 ft. 

24 

12. 0 

30 ft. 

43 

21.5 

40 ft. 

59 

29.5 

5 0 ft. 

62 

31. 0 

60 ft. 

68 

34. 0 

70 ft. 

74 

37. 0 

80 ft. 

72 

36. 0 

90 ft. 

71 

35.5 

100 ft. 

70 

35. 0 


Range 0-500 Reverse Course 

Seismograph 


readings Milliseconds 


100 ft 

72 

36. 0 

90 ft. 

71 

35.5 

80 ft. 

70 

35.0 

70 ft. 

66 

33. 0 

60 ft. 

66 

33. 0 

50 ft. 

63 

31.5 

40 ft. 

60 

30. 0 

30 ft. 

51 

25.5 

20 ft. 

27 

13.5 

10 ft. 

2 

1 . 0 

Range 0-500 Reverse Course 


Seismograph 



readings 

Milliseconds 

100 ft. 

73 

36.5 

90 ft. 

73 

36.5 

80 ft. 

69 

34.5 

70 ft. 

64 

32. 0 

60 ft. 

63 

31. 5 

5 0 ft. 

61 

30.5 

40 ft. 

55 

27.5 

3 0 ft. 

43 

21.5 

20 ft. 

26 

13.0 

10 ft. 

6 

3. 0 

Range 

0-500 Reverse Course 


Seismograph 



readings 

Milliseconds 

100 ft. 

70 

35. 0 

90 ft. 

70 

35.0 

80 ft. 

71 

35.5 

70 ft. 

74 

37. 0 

60 ft. 

69 

34. 5 

5 0 ft. 

63 

31.5 

40 ft. 

59 

29.5 

3 0 ft. 

43 

21.5 

20 ft. 

21 

10.5 

10 ft. 

7 

3.5 
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Canada de Alfaro, Site 61, Course 1, Range 0-5 00 



Seismograph 



Seismograph 

Milliseconds 


readings 

Milliseconds 


readings 

10 ft. 

8 

4. 0 

100 ft. 

80 

40. 0 

20 ft. 

25 

12.5 

90 ft. 

75 

37.5 

3 0 ft. 

43 

21.5 

80 ft. 

69 

34.5 

40 ft. 

54 

27. 0 

70 ft. 

64 

32. 0 

5 0 ft. 

58 

29. 0 

60 ft. 

60 

30. 0 

6 0 ft. 

61 

30.5 

50 ft. 

56 

28. 0 

70 ft. 

64 

32. 0 

40 ft. 

52 

26. 0 

80 ft. 

70 

35. 0 

30 ft. 

43 

21.5 

90 ft. 

76 

38. 0 

20 ft. 

20 

10. 0 

100 ft. 

80 

40. 0 

10 ft. 

4 

2. 0 

Canada de Alfaro, Site 61, Course 2, 

Range 

0-500 Reverse Course 


Seismograph 


Seismograph 

Milliseconds 


readings 

Milliseconds 


readings 

10 ft. 

2 

1 . 0 

100 ft. 

73 

36.5 

20 ft. 

29 

14. 5 

90 ft. 

68 

34. 0 

3 0 ft. 

41 

20.5 

80 ft. 

62 

31. 0 

40 ft. 

48 

24. 0 

70 ft. 

58 

29. 0 

5 0 ft. 

50 

25. 0 

60 ft. 

53 

26.5 

60 ft. 

54 

27. 0 

5 0 ft. 

49 

24.5 

70 ft. 

58 

29. 0 

40 ft. 

47 

23.5 

80 ft. 

63 

31.5 

30 ft. 

40 

20. 0 

90 ft. 

68 

34. 0 

20 ft. 

30 

15.0 

100 ft. 

73 

36.5 

10 ft. 

2 

1 . 0 

Site - 

Leon - Silao, 

Course 1, Range 0-500 

Reverse Course 


Seismograph 



Seismograph 

Milliseconds 


readings 

Milliseconds 


readings 

10 ft. 

5 

2. 5 

100 ft. 

13 8 

69. 0 

20 ft. 

23 

11. 5 

90 ft. 

136 

68. 0 

30 ft. 

37 

18. 5 

80 ft. 

121 

60.5 

40 ft. 

47 

23.5 

70 ft. 

115 

57.5 

50 ft. 

66 

33. 0 

60 ft. 

97 

48.5 

60 ft. 

95 

47.5 

5 0 ft. 

66 

33.0 

70 ft. 

No reading 


40 ft. 

46 

23. 0 

80 ft. 

121 

60.5 

3 0 ft. 

40 

20. 0 

90 ft. 

135 

67. 5 

20 ft. 

23 

11.5 

100 ft 

140 

70. 0 

10 ft. 

5 

2. 5 
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Site - 

Leon - Silao, 

Course 2, 

Range 0 - 500 

Reverse Course 


Seismograph 



Seismograph 



readings 

Milliseconds 

readings 

Milliseconds 

10 ft. 

5 

2.5 

100 ft. 

70 

35. 0 

20 ft. 

25 

12.5 

90 ft. 

66 

33. 0 

3 0 ft. 

38 

19. 0 

80 ft. 

59 

29.5 

40 ft. 

42 

21. 0 

70 ft. 

54 

27. 0 

5 0 ft. 

45 

22.5 

60 ft. 

50 

25. 0 

60 ft. 

50 

25. 0 

5 0 ft. 

46 

23.0 

70 ft. 

54 

27. 0 

40 ft. 

41 

20.5 

80 ft. 

59 

29.5 

3 0 ft. 

38 

19. 0 

90 ft. 

65 

32.5 

20 ft. 

25 

12. 5 

100 ft. 

71 

35.5 

10 ft. 

5 

2.5 


F ormulae 

lo Depth formulae, for constant-velocity layers are given below where: 

K = constant 

n = number of a given layer or strata 
V = velocity 
d = depth 

C = critical distance, or the distance from the geophone to the point 
at which the slope of a given velocity line changes 
t = time 

g = width of a depression 
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2. Depth formula for sloped ground surface with pyramids: 

Possible only if the angle of the walls are known. 

3. Flat ground surface with a platform: 

Again, the solution is possible only if the angle is known or if a 
true reversed course was run. 

4. Height of a wall, as in Case 2. 

h = height of wall 

_ dj (v 2 2 - 2 V] 2 ) 
h= 

This formula only applies if the velocity of the second stratum is at 
least \f2T times the velocity of the first stratum. 

5. Width of a depression, as in Cases 4 and 5 in the following section: 

g = (v?H) + _ C 1 V ] 

c' is any point after which the floor is picked up beyond the 
depression, at t Q . 


General Operating Principles 

1. A straight line indicates that the hammer's path is parallel to the 
stratum of the quickest path. 

2. A curved graph will veer up (away from the straight) if the quickest 
path is veering away from the hammer path. 

3. A curved graph veers down if the quickest path is getting closer 
progressively to the hammer path. 
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Hypothetical Archeological Cases and Their Time-Distance Graphs 

Case 1. A flat ground surface and homogeneous strata parallel to the 
surface. 


Hypothetical Stratigraphy 


XX x xxxxx 

X X * X 

X X X XX XXX 

oooooooo 

O 0 O O 0 ^ O o ft 

o o o o o o o ~ 


/ / ✓ / // /^ /// 
/ / / // // / ' /< ' 


Hypothetical Time-Distance Graph 



Case 2. A flat ground surface with the strata parallel to the surface, 
and a wall and floor of stone. 


Hypothetical Stratigraphy 


Hypothetical Time -Distance Graph 



O O G & * O O O O O 

o o O o a o O G o c 



a. For exact numbers, the wall and floor must be about the same density. 

b. Otherwise the above graph applies to all cases and the wall occurs 
between the bounds of the dip. 

c. The only way to differentiate different materials used for the wall 
and floor will be in the sharpness of the curve. If the floor is of 
packed earth and the wall of stone, the curve will approach a 
straight line. 
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Case 3. Similar to case 2 above, except that the transect is following 
a wall which ends, and the path jumps down to pick up a floor which 
is lower than the wall 


Hypothetical Stratigraphy 


"y —J - 5 T 7 X . X X 

* * x X_ X 


X X 

X x XX 
X X X X 


f r.i i vLin x 

y X X X HJ r l I n-rnmi 
X < X x x x x y X x. 


ooooooo o 

o o o o o o o o 


Hypothetical Time-Distance Graph 



Case 4. A flat ground surface with the strata parallel to the surface. 
A depression filled with approximately the same material as the first 

stratum. 


Hypothetical Stratigraphy 


y X X X X X 

X X X X X X 
X > X X X X X 

1 11 ctii x x pcxxxn 
o o °\x * y o o o 

O O O o ~ 0 " O o o o 


Hypothetical Time -Distance Graph 



A graph of this type will apply to 
any interruption in the stratum 
not necessarily a man-made _ 
depres sion. 


a. 
































> o 
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ase 5. A flat ground surface with the strata parallel to the surface, 
depression filled with a layer of rocks at the bottom. 


Hypothetical Stratigraphy Hypothetical Time-Distance Graph 



a. This is the same general graph as 
Case 4 above, but the third slope 
will be shorter. This case will 
rarely occur. 

Case 6. A flat ground surface, one platform with one pyramid on top. 


Hypothetical Stratigraphy Hypothetical Time-Distance Graph 



a. A single platform only, or any 
number of pyramids will give a 
variation of this graph, with 
more or less curves. 
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Case 7. A sloped ground surface, and a floor which is not parallel 
to the surface. The course runs downhill. 

Hypothetical Stratigraphy Hypothetical Time-Distance Graph 




Case 8. A sloped ground surface, and a floor which is not parallel to 
the surface. The course runs uphill. 

Hypothetical Stratigraphy Hypothetical Time-Distance Graph 




Distance 
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Case 9. A sloped ground surface, 1 single platform. The course 
begins at the top and runs downhill. 


Hypothetical Stratigraphy 



Hypothetical Time-Distance Graph 



Case 10. A sloped ground surfac 
The course runs uphill. 

Hypothetical Stratigraphy 



, and any number of pyramid walls. 


Hypothetical Time-Distance Graph 


















- 111 - 

VIII. Results 


List of Critical Distances, Velocities and Depth Values 

Layer Critical 


Site Site # 

Course 

i 

of 

Strata 

Distance 
in feet* 

Velocity 

in ft/sec Depthinft 

Rio Guanajuato 


1 

2 

20. 0 

303. 0 
1,013.0 

7. 

6. 

3 

2 

Rio Guanajuato 


3 

1 

2 

20. 0 

541.0 
387. 0 



Rio Guanajuato 


4 

1 

20. 0 

741. 0 

1,685. 0 



Canada de Alfaro 

43 

1 

1 

20. 0 

709. 0 

7,207.0 

9. 

1 

Canada de Alfaro 

43 

3 

1 

30. 0 

1,124.0 

12. 

0 



2 

50. 0 

5,128.0 

25. 

5 




3 

70. 0 

10,526.0 

31. 

8 




4 

90. 0 

11, 764.0 






5 


2,273.0 



Canada de Alfaro 

22 

1 

1 

30. 0 

1, 053.0 

8 • 

1 



2 

50. 0 

1, 905.0 

23. 

4 




3 

90. 0 

5, 714.0 








4, 865.0 



Canada de Alfaro 

22 

Reverse 

1 

40. 0 

1, 229.0 

16. 

7 





7, 000.® 



Canada de Alfaro 

22 

2 

1 

50. 0 

1, 052.0 

2. 

1 





1, 067.0 



Canada de Alfaro 

22 

Reverse 

1 

30. 0 

852. 0 






2 

40. 0 

100.0 






3 

60. 0 

909. 0 






4 

80. 0 

1, 250. 0 






5 


436. 0 



Canada de Alfaro 

22 

3 

1 

40. 0 

1,111. 0 






2 

50. 0 

5, 000.0 






3 

60. 0 

2, 500.0 






4 

100. 0 

2, 911. 0 






5 

130. 0 

1, 224.0 






6 


1, 755.0 



Canada de Alfaro 

22 

Reverse 

1 

2 

3 

4 

20. 0 

30. 0 

80. 0 

1, 250.0 

690. 0 
3,125.0 

1, 818. 0 



Canada de Alfaro 

22 

4 

1 

2 

3 

20. 0 

50. 0 

625. 0 

2,222.0 
3,290.0 

7 

, 17 

.5 

. 7 

Canada de Alfaro 

22 

Reverse 

1 

2 

3 

4 

20. 0 
40. 0 
60. 0 

741.0 

1,666. 0 
2,222.0 

3, 810.0 

6. 2 

12. 8 
24. 7 


* The critical distance is the distance from the geophone to the point at 

which the slope of a given velocity line changes. 
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Site 

Site # 

Course 

Layer 

of 

Strata 

Critical 
Distance 
in feet 

Velocity 
in ft/sec 

Depth in ft« 

Canada de Alfaro 

11 

1 

1 

2 

3 

4 

5 

20. 0 

40. 0 

60. 0 

70. 0 

1,538. 0 
2,222.0 

3,636.0 
645.0 
5,454.0 

4.3 

12. 7 

Canada de Alfaro 

11 

2 

1 

20. 0 

1, 250.0 

5. 9 




2 

50. 0 

2, 610.0 

18.6 




3 

80. 0 

5, 000.0 

30. 8 




4 


6,666.0 


Canada de Alfaro 

11 

Reverse 

1 

20. 0 

1, 538.0 

5. 8 




2 

50. 0 

2, 858.0 

18.1 




3 

80. 0 

5,455.0 

22. 8 




4 


5, 714. 0 


Canada de Alfaro 

10 

i 

1 

20. 0 

1, 176. 0 

5.5 




2 

70. 0 

2, 174.0 

32.6 




3 


10, 500. 0 

. 

Canada de Alfaro 

10 

2 

1 

30. 0 

817.0 

12. 7 




2 

80. 0 

5, 000.0 

42.7 




3 


20, 000. 0 


Canada de Alfaro 

10 

Reverse 

1 

30. 0 

817. 0 

13.3 




2 


6,666.0 


Canada de Alfaro 

57 

i 

1 

40. 0 

1, 224. 0 

17. 0 




2 


7, 500. 0 


Canada de Alfaro 

57 

Reverse 

1 

30. 0 

1, 080. 0 

6. 9 




2 

40. 0 

1, 667. 0 

19.9 




3 


6, 666. 0 


Canada de Alfaro 

57 

2 

1 

40. 0 

1,154. 0 

18. 0 




2 


10, 909.0 


Canada de Alfaro 

61 

1 

1 

30. 0 

1,143.0 

7. 2 




2 

40. 0 

1, 818. 0 

18. 0 






4, 615.0 


Canada de Alfaro 

61 

Reverse 

1 

20. 0 

1, 250. 0 





2 

30. 0 

870. 0 





3 

40. 0 

2, 222. 0 





4 

70. 0 

5, 000.0 





5 


3, 750. 0 



Canada de Alfaro 

61 

2 

1 

2 

3 

4 

5 

6 

20. 0 

30. 0 

40. 0 

50. 0 

90. 0 

741.0 

1, 667. 0 

2, 857.0 

10, 000.0 

3,478. 0 

4,000.0 

6. 2 
12. 7 
24.8 

Canada de Alfaro 

61 

Reverse 

1 

20. 0 

714. 0 

6.2 




2 

40. 0 

2,353.0 

22. 0 




3 

50. 0 

10, 000.0 





4 

80. 0 

4, 615.0 





5 


3,636.0 


Leon-Silao 


i 

1 

30. 0 

1, 142. 0 

5. 0 




2 


1,429. 0 


Leon-Silao 


2 

1 

20. 0 

1, 000. 0 

6.4 




2 

40. 0 

2, 352. 0 

17.1 




3 

70. 0 

5, 000.0 





4 


3,535.0 
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Site - Rio Guanajuato Course 3 Range 0-500 
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Site - Canada de. Alfaro Course 1 Site 43 Range 0 - 100 
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Canada de Alfaro Site 22 Course 2 Range 0 - 500 
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Canada de Alfaro Site 22 Course 3 Range 0 - 500 
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Canada de Alfaro Site 22 Course 4 Range 0 - 500 
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Canada de Alfaro Site 11 Course 1 Range 0 - 500 
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IX. Structural Conclusions 

The structural interpretations that follow were based on the hypo¬ 
thetical cases, the real time-distance graphs as plotted by the computer, 
and general knowledge of the sites. 

The assumption supporting these interpretations is that the wave 
path will follow the structure and will reflect the surface of that struc¬ 
ture. For example, inCase 2 of the hypothetical cases, the path fol¬ 
lows the wall (assuming that it is wide enough), and the curve in the 
time-distance graph reflects that wall. 


A * ^ 

—~r. 

A 

x X x 

X >C y 


X x x x 

\J v/ 

^ a Y 

n — l — m — r — r 

.J-J 

A A X 

h—mi t~i 


Q 0 C 0 V 0 C & 0 O & 

0 © Q O Q Q O & & & 



The path follows the floor in a straight line of a certain velocity* 

As it follows the wall up towards the surface, the velocity speeds up, 
thus the veering down from the straight line. As it moves down the wall 
further away from the surface, the velocity begins to slow down again, 
and as it reaches the level floor, the velocity again resumes the same 
speed as before. The operating principles (page 91) apply to any case; 
the shape of the graph reflects the shape of the structure. 

It must be made clear that the following interpretations are prob¬ 
ably not the only ones possible, but these seemed the M best fit ,r to real 
data, and seemed most probable with the interpreter's knowledge of the 
sites themselves. 

Structural Interpretations 

Site Rio Guanajuato - Course 1, Reverse Course: 



The course ran down hill, and showed a layer of disturbed topsoil 
down to a depth of 7.345 feet directly under the geophone. At this point 
it appears to be picking up a single platform of shale or limestone. 
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Forward Course 



A layer of disturbed topsoil down to a depth of 7.345 feet directly 
under the geophone, then picking up what appears to be a smaller plat¬ 
form on top of a larger one. The structure continues after the present 
day surface levels off. 

Site Rio Guanajuato - Course 3 



The forward and reverse courses show basically the same structure. 
The depth of the topsoil is unknown. A structure is detected which appears 
to be a larger platform on top of a smaller one. Another platform con¬ 
tinues where the present-day level of the surface begins to flatten out. 

Site Rio Guanajuato - Course 4, Forward Course 


Geophone 

^ To ps61 



The forward course shows a layer of topsoil to a depth of 6. 238 feet 
directly under the geophone. Then a platform is picked up, and the 
mound on the very top could either be another small platform or a pile 
of loose boulders. What appears to be two platforms are picked up on 
the other side, but this could be only a reflection of the damage to the 
structure where a hole has been dug and several stones removed. 



The reverse course also shows a layer of topsoil to a depth of 
6 23 8 feet. In this case, it seems to detect a platform on top of at 
least one other platform. The smaller mound on top appears to con¬ 
sist of loose dirt from the hole dug on the other side. Again, it is 
probable that the damage near the end of the course is being picke up, 
affecting the readings. 
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Canada de Alfaro : 


Site 43 - Course 1 




This course shows a layer of topsoil to a depth of 9. 06 feet, then 
appears to pick up three platforms with very badly weathered sides. 



The alternative interpretation also shows a layer of topsoil to a 
depth of 9. 06 feet, but with a single platform on top of a platform after 
the surface levels off, somewhere around 60 feet. 

Site 43 - ^^ r* ^^- 



The forward course shows a layer of topsoil to a depth of 12 feet. 
Then it appears to pick up a platform leading into another platform, all 
under a relatively flat ground surface. Where the surface begins to slope 
■upward, another platform shows. Its nature remained unknown because 

the COUrSe on^ v m f ?.n r\-n a m nnliin cl A-no 



Geophone 


Reverse Course 


rock 


The reverse course shows no structure at all in the plaza, only 
normal homogeneous layers. However, a portion of a platform is picked 
up where the course begins to go uphill at a distance of 70 to 90 feet. 
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Site 22 - Course 1 


Geophone 



The forward course in this case seems to be the most probable. 

It appears that directly under the geophone, there is no structure. Here, 
there is a layer of topsoil to a depth of 8. 05 feet, then a layer of sandy 
clay or gravel to a depth of 23.4 feet. Somewhere around 5 0 feet, a 
platform is picked up, which goes down steps onto another platform in 
a sunken plaza. An alternative interpretation would be that there is no 
structure at all in this area, but that which is detected is downslope 
layering over a prehistoric sunken plaza dug into the ground. 



The forward course again seems to be the most likely case. A 
layer of topsoil to a depth of 2.10 feet directly under the geophone over- 
lies a platform in the sunken plaza. The platform becomes an unknown 
number of pyramids as the slope of the ground surface goes uphill. 

Site 22 - Course 3 



The most likely interpretation for the forward course is that a 
layer of topsoil is picked up to an unknown depth. Then the pyramid 
wall is detected to the side of the transect rather than underneath it. 
Damage to the wall due to pothunting is also detected at the beginning 
of the course. Somewhere around 100 feet, the pyramid ends and the 
path jlamps down to pick up the platform under the pyramid. The plat' 
form continues approximately 13 0 feet. The course then begins to 
follow steps down the side of the platform. 
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Site 22 - Course 4, Forward course 


Geophone 


-¥r 

\ 


7.49 ft. topsoil 

\ 

~ / q r. gravel or 

D - 17.69 ft. , . . 

r Z glacial fill 


shale, sandstone, or limestone 


Reverse course 


Geophone 


-i-p 


_ / , „ - topsoil 

D = 6.19 ft. * 

> 

r 

sandy clay 

1Z. 77 ft. or gravel 


gravel, glacial fill 
24. 73 ft. or rock talus 


shale, sandstone, or limestone 


Both the forward and the reverse course show a natural homogeneous 
layering, and no structures. The only difference is that the reverse course 
shows one more layer than the forward course, and differing depths. 


Site 11 - Course 1 


Geophone 


JL 


Dj= 4. 26 ft. sandy clay 


D 0 = 12.66 ft. f* ave1 ’ S 1 ^ 1 

2 fill or rock talus 



shale or sandstone 


Both the forward and reverse courses show no structure, but a 
depression somewhere between 50 and 80 feet* The nature of the de¬ 
pression is unknown. 
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Site 11 - Course 2 


Geophone 

T 

Dj= 5.47 ft. sandy clay 

_> 

gravel or 

D “ 18. 06 ft. glacial fill 

r * 


shale or sandstone 


Again, both courses show no structures, but homogeneous layers. 
The curve near the end of the course probably indicates a downslope in 
the layering. 

Site 10 - Course 1 



The most likely possibility here is that a pyramid to the side of the 
transect was detected as the transect moved closer to the pyramid. 

The structure ended and the natural layer of shale or granite was picked 
up at an unknown depth. 

Site 10 - Course 2 

The transect is most likely picking up a structure to the side of it, 
rather than underneath it, but this interpretation is highly speculative. 

Both Course 1 and Course 2 are probably invalid, as there was a 
great deal of noise at the site, making accurate readings difficult. 

Site 57 - Course 1 



The reverse course shows a layer of topsoil to a depth of 6. 93 feet. 
Then what appears is the wall of a pyramid to the side of the transect as 
the transect moves towards that wall. A small portion of a platform in 
the plaza is indicated, then another platform begins as the slope of the 
surface goes uphill. 
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Site 57 - Course 2 


\ 

\ 

\ 

i 

* Geophone 

Both courses show a layer of topsoil to an unknown depth. Again, it 
seems that the path is picking up the wall of a pyramid to the side of the 
transect, as it moves away from the wall. Then it jumps down to pick 
up a platform in the plaza, which continues for a short way after the 
ground surface begins to slope downward. 



Site 61 - Course 1, Forward course 
Geophone^r 


T 


D^= 7 ft. topsoil 


D = 17. 9 ft. sandy clay 

Lt 


shale or sandstone 

The forward course shows homogeneous layers and no structures. 


Reverse course 



The reverse course would seem to indicate that the wall of a pyramid 
is picked up to the side of the transect, which is moving towards the wall, 
then the side of a platform as the transect is going down while the ground 
level remains flat. 


Site 61 - C 
Geophone 


ourse 2, Forward course 


> 

r 

Dj=6.19ft. topsoil 

> F 

D^=12. 69 ft. sandy clay 


D = 24. 80 ft. glacial fill 


shale or limestone 


The forward course indicates homogeneous layers and no structures. 
The curve near the end of the graph could be caused when the transect 
moves onto the cornfield where the layer of topsoil thickens, i.e. the 
wave slows down. 
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Site 61 - Course 2, Reverse course 



D^=7.31 ft. topsoil 


gravel or 
D =21. 94 ft. rock talus 

\f 2 _—— 


shale or limestone 


The reverse course also shows homogeneous layers, but with one 
less layer than the forward course. The same curve, only less pronounced, 
is shown at the end of the graph. 

Site Leon-Silao - Course 1 



A layer of top soil is shown to a depth of 5. 01. This is probably 
due to the fact that the geophone was placed on what appeared to be a 
mound of old backdirt. Then a platform seemed to be picked up which 
continued after the ground surface began to slope uphill. After 7 0 feet, 
nothing can be said about the structure because of damage from an old 
excavation. 

Site Leon-Silao - Course 2 



gravel or 
D2=17. 09 ft. rock talus 


shale or sandstone 


This interpretation would indicate homogeneous layers and no 
structures. The curve caused by the slowing velocity beginning at 
about 5 0 feet would be due to the fact that the course goes into a corn¬ 
field where the layer of topsoil would be thicker, causing the wave to 
slow down as it passed through this layer. 

Site Leon - Silao - Course 2, Alternative interpretation 



This interpretation would show the platform continuing under the 
pile of old backdirt, the platform ends, and the path picks up the natural 
bed of shale or sandstone at a lower level under the cornfield. 
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X. Final Conclusion 

Effective use of the seismograph has previously been based on the 
assumption that subsurface strata are homogeneous, each successive 
layer being of a higher velocity. However, the investigators of this 
project were not only interested in the general stratigraphy of the area, 
but in whether or not subsurface anomalies, i. e. , structures, could be 
detected without having to spend a great deal of time, money, and man¬ 
power for excavations. This involved dealing with layers which were not 
homogeneous, and which were not laid down with uniformly decreasing 
densities from the lower levels up. Previous investigations have dealt 
with irregularities or anomalies as 'problems' that must be recognized 
and accounted for in the interpretations, but which were not interpreted 
themselves. Therefore, except for the Rainy expedition, ^ this investi¬ 
gation is the first of its kind. 

Is it possible then to detect subsurface structures with a portable 
refraction seismograph? The answer appears to be 'yes, ' but with a 
great many limitations and reservations. 

The basic assumption when using the seismograph for archaeological 
purposes is that a wave will pick up a structure and will follow that struc¬ 
ture for the length of the course because the material of a structure will 
be of a much higher velocity than any layer near that same depth, and 
thus will be the quickest path for the wave to take. Therefore, the path 
will never pick up layers below a structure, unless the structure ends. 
The most common problem in most seismographic work, that of the 
hidden layer, becomes inrrelevant for our purposes. 

Methodology of field procedures 

In areas where the subsurface material is unknown, a systematic 
method of placing the courses, and the number of courses run should 
be developed in order to reduce the degree of arbitrariness. The best 
approach seems to be to run four connecting courses, running north, 
south, east, and west. If more courses are run, for example, systema¬ 
tically covering the entire site, the interpreter could never be sure from 
where he was picking up structures — from the side, or below the tran¬ 
sect, and the plethora of detail becomes too confusing. 

Spacing of 10 foot intervals is small enough when one simply wishes 
to determine the slope of a homogeneous layer. However, when the goal 
is to follow a specific structure, smaller spacings of 5 feet or less are 
desirable. This allows the investigators greater accuracy in drawing 

-r- 

Rainy, "Archaeology and Its New Technology, " Science , 153: 1490. 

In the search for the lost city of Sybaris, Rainy reported that seismic 
instruments were "not suitable for detecting small archaeological fea¬ 
tures at relatively shallow depths. " 
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the graphs by helping to determine if a specific portion was a straight 
line or a curve. 

I 

Running a reversed course for this type of problem is essential, 
not for determining average velocities in order to draw a continuous 
layer, but to validate the shape and slope of the subsurface structure. 

Assumptions and limitations for interpretations 

In drawing the graphs, it was often difficult to know when to draw 
a straight line between points — assuming that those points varying from 
the straight line were due to either variations in the velocity, or varia¬ 
tions in the depth of the material through which the wave passes — or 
whether to draw a curve by picking up every point in the curve. Error 
of the machine and error of the computer in drawing the graphs by 
rounding off the readings should also be taken into account. 

An intelligent interpretation of velocities on a time-distance graph 
always requires some knowledge of the area tested. Information about 
what may be expected will reduce the number of courses required 
for checks. It is difficult to pin-point a specific velocity for a certain 
material since the material itself may vary with depth or geographic 
location. For example, weathered, broken, or seamed rock has a slower 
velocity than solid rock; clean gravel has a slower velocity than clayey 
gravel. 

In interpreting the graphs, the knowledge of the mound's appearance 
and the most likely subsurface structures were always kept in mind. 

The general rule followed was that wherever the course was run over a 
sloping ground surface, the resulting line would probably be a curve — 
indicating a structure - unless the points approached a nearly straight 
line, or unless the surface of the mound was parallel to the general 
ground surface. Other data was used to help interpret the graphs, such 
as excavation data, reconstruction drawings, site forms, and general 
knowledge of the area. 

In cases where the time readings varied at the same distance, a 
line through both points was drawn, and the interpretation allowed for 
both cases. 

The course ended at an arbitrary point, not necessarily at the end 
of a structure, making it very difficult to interpret the end of the graph 
correctly. Also, it is often difficult to decide if that which is detected 
by the seismograph is a structure directly underneath the transect, to 
the side of it, or running ninety degrees to the course. The only way 
to decide would be to match up the profiles of cross-connecting courses. 

It must be remembered that depth calculations give only the depth 
of layers directly under the geophone. A different formula must be 
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used to calculate the depth to the structure, and this can only be done 
if the angles of the structure are known. The general depth formula 
was used on the graphs until a point was reached where a layer showed 
a lower velocity than the preceding layer. Otherwise, the velocities of 
the curves were averaged and used to give the depth of the previous 
layer. 

The exact location of the seismograph must be known and recorded 
to help create an accurate interpretation. After interpretations are 
made, they should be tested by excavating in those areas where struc- 
tures were indicated to validate the interpretations. 

It appears that this method for testing for archaeological material 
may prove to be useful, but further experimentation and refinement of 
the technique needs to be done. Only in cases where one is working on 
a specific mound and has other data with which to correlate the seismo- 
graphically generated data can intelligent interpretations be made as to 
the shape and location of subsurface structures. 

Finally, I would like to express my deep appreciation to all of those 
people who have helped me on this project, especially Dr. Zubrow, who 
served as my advisor, Keith Kintigh, who wrote the computer programs, 
John Justeson, who devised the formulas and explained many of the prin¬ 
ciples to me, and Andrew Willard, who edited this manuscript. Need¬ 
less to state, any errors are the responsibility of the author. 
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I ntroduction 

During the summer of 1972, an archaeological expedition in 
Guana iuato Mexico was directed by Dr. Ezra Zubrow of Stanford 
Snive sity This work consisted of extensive surveying, with very 
UmUed excavation. Seventy-three separate mounds, prospective ball- 
courts some plazas, and a road were surveyed. Random system 
sampling was applied to each mound to collect pottery sherds lithic 
material, and other artifacts. Areal measurements were t ~ 
suring the top, sides, and border of each mound. A height measur 
ment was also taken, and the shape of each mound was approximate . 


Method 

This essay attempts to determine if there are any relationships 
between the amount of labor invested in each site an the sizeof th 
sites, and to see if there are any general groupings of the site by 
size This will be done by calculating the volume of each site a 
by Z determining how many man-days of labor were invested in each. 

Using the areal measurements, a volume was calculated for each 
site Formulas for circular and oval cones, and pyrami s yi 
Volume in terms of cubic feet. This number was changed to cubic 
mirrs because the following experiments in labor investment were 
conducted in the metric system. 

Charles Erasmus (1965) conducted field experiments totest theamoint 
of labor invested in the construction of Maya ceremonial J 

S-Ei; -sc sxs; =■,:==* 

of 0 fill k would f thefrequirltbout three man-days of labor for 

content and he allows one man-day for the dirt content of 500 kilos. 

He arrives at a figure of one man-day for the excavation of earth and 
rock, per cubic meter. This gives a total of five man-days needed 
for the excavation and transportation of one cubic meter of i , 
taining both rock and earth. 

From what is known, the mounds at Guanajuato are rubble filled, 
some contain walls of masonry, and therefore the total figures would 

have to include the amount of time invested in construction of masonry. 

Erasmus had a Maya stone mason construct a test wall I 
that 1400 kilos of rock, 300 kilos of sascab-the u = ed “ 

mortar, and 100 kilos of lime were used per cubic meter o rnasonry. 
He estimated that one-half a man-day per cubic mete Y 
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the sascab content would be maximal, and he uses the previous igure 
of four man-days for the excavation and carrying of the rock. He a - 
lows three and one-half man-days for the manufacture of the lime used 
in the mortar. This gives a total of eight man-days per cubic meter. 

Next, he had two Maya building contractors calculate the time 
invested in the stone walls of houses they were building. He arrived 
at a figure of four man-days per cubic meter of construction. This 
gives a total of twelve man-days needed per cubic meter of excavation, 
transportation, and construction. 

By taking the average between the twelve man-days which includes 
construction, and the five man-days for excavation and transportation, 
I decided a figure of eight and one-half man-days invested per cu ic 
meter of fill would be appropriate. 


Since the volume was calculated in terms of cubic meters, it was 
multiplied by 8.5 man-days per cubic meter. This index is the num er 
of man-days invested per site. This figure was then multiplied by 5 
man-hours per man-day to give the total number of man-hours in¬ 
vested per site. This number was multiplied by 50 percent to allow 
for possible room space. 


The sites were rank-ordered by thousands of man-days needed 
for construction. Next, a graph was constructed showing the number 
of sites per each thousand man-days. It showed that there were four 
distinct groupings of sites - a large number between 0 -4000 man-days, 
fewer between 4500 - 7500, one between 12, 000 - 14, 000 man-days, an 
two between 21, 000 - 23, 000 man-days. Using these four groups, the 
chi-square test was applied. It showed that the probability of these 
groupings occurring by chance is so minimal that it may be legitimately 

ignored. 


Data is given in tiie following tables. 


Data 
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Table I 

Calculations of Volumes, Man-days, and Man-hours 



Volume 
in feet^ 

Volume 

in meters 

Number of Number of 
man-days man-hours 

Number of 

man-hrs. t 2 

Site #1 

2379. 2 

67.37 

572.65 

2863.25 

1431.61 

Site #2 

5265.1 

149. 09 

1267.29 

6336.48 

3168. 24 

Site #3 

863 9.1 

244.62 

2079.35 

‘10396.76 

5198.38 

Site #4 

4569. 25 

129.00 

1096.5 

5482.5 

2741.25 

Site #5 

9448.53 

267.55 

2274.17 

11370. 87 

5685.43 

Site #6 

540.13 

15.29 

130. 00 

650. 02 

325.01 

Site #7 

1155 . 2 

32.71 

278.04 

1390. 23 

695. 11 

Site #8 

12800.00 

362.45 

3080. 84 

15404. 21 

7702.10 

Site #9 

209.38 

5.92 

50.39 

251.97 

125.98 

Site #10 

5655.56 

160.14 

1361. 24 

6806.20 

3403.10 

Site #11 

366.00 

10.36 

88. 09 

440.46 

220.23 

Site #12 

137.30 

3. 88 

33. 04 

165.23 

82.61 

Site #13 

6900. 83 

195.40 

1660. 96 

8304.83 

4152.41 

Site #15 

14389.69 

407.46 

3463.46 

17317.33 

8658.66 

Site #16 

4599.84. 

130.25 

1107. 13 

5535.69 

2767.84 

Site #19 

8267.46 

234.10 

1989.90 

9949.51 

4974.75 

Site #20 

5927.24 

167. 83 

1426. 62 

7133 . 00 

3566.57 

Site #21 

1908. 39 

54. 03 

459.33 

2296.66 

1148.33 

Site #22 

92799.08 

2627.75 

22335.89 

111679.48 

55839.74 

Site #23 

25025.82 

708.64 

6023.48 

3 0117.43 

15058. 71 

Site #26 

89427.66 

2532.28 

21524.42 

107622.13 

53811.06 

Site #27 

15339. 83 

434.37 

3693.15 

18460. 78 

9230.39 

Site #28 

23028. 50 

652.08 

5542.75 

27713.75 

13856. 87 

Site #29 

7479.97 

211.80 

1800.36 

9001.80 

4500. 90 

Site #30 

5808.00 

164.46 

1397. 93 

6989.66 

3494.83 

Site #31 

9014.18 

255.25 

2169. 63 

10848.15 

5424.07 

Site #32 

22543.90 

638.36 

5426.11 

27130.56 

13565.28 
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Volume 

in feet 

Volume 

3 

m meters 

Number of 
man-days 

Number of 

man-hours 

Number of 
ma'n-hrs.T 2 

Site #35 

27235.08 

771. 20 

6555.23 

32776.18 

16388.09 

Site #36 

3266.66 

92.50 

786.25 

3931.27 

1965.63 

Site #37 

4100.00 

116. 09 

986.83 

4934.16 

2467.08 

Site #41 

16255.41 

460. 29 

3912.52 

19562.64 

9781.32 

Site #42 

1494. 29 

42.31 

359.66 

1798.31 

899.15 

Site #43 

21850. 00 

618. 71 

5259. 09 

26295.48 

13147.74 

Site #45 

10598.61 

300.11 

2550. 98 

12754.94 

6377.47 

Site #47 

6893.65 

195.20 

1659.23 

8296.19 

4148. 09 

Site #48 

14903.37 

422. 01 

3587. 10 

17935.52 

8967.76 

Site #5 0 

1223.18 

34. 63 

294.40 

1472. 04 

736.02 

Site #51 

50061. 71 

1417.57 

12049.39 

60246.99 

3 0123.49 

Site #53 

7314. 91 

207.13 

1760. 63 

8803.16 

4401.5 8 

Site #54 

2733.20 

77.39 

657.85 

3289.28 

1644.64 

Site #55 

3889.64 

110.14 

936.20 

4681.00 

2340. 50 

Site #57 

4062.88 

115. 04 

977.89 

4889.49 

2444.74 

Site #58 

2015.42 

57. 07 

485.09 

2425.46 

1212. 73 

Site #5 9 

4069.17 

115.22 

979.41 

4897.06 

2448.53 

Site #60 

75.41 

2.13 

18.15 

90. 75 

45.37 

Site #63 

1323.09 

37.46 

318.45 

1592. 27 

796.13 

Site #69 

6403.99 

181.33 

1541.38 

7706.91 

3853.45 
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Table II 

Rank Ordered Sites by Number of Man-days 



-23000 

-22000 

-21000 
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Table III 

The Number of Sites per Each Thousand Man-days 



X 


X 


r 

1—» I—» !— 1 I—* 

-J 00 NO 

o o o o 
o o o o 
o o o o 



Thousands of man-days 


1-24000 
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Table IV 

Results of the Chi-square Tests 

' 

Test #1 - x 2 = 83.24 
p = 1 % 

2 

Test #2 - x = 320653.82 
P =1% 

I 

Test #3 - x^ = 40.51 
P = 1% 

•JU 

•nr* 

Three different chi-square tests were applied to the same 
data, all showing the groupings to be significant. 


Table V 

Totals of Volumes, Man-days, and Man-hours 

1. Total volume in cubic feet - 573, 363* 98 ft. 3 
2„ Total volume in cubic meters - 16, 235. 70 m. 3 

3, Total man-days 

Group I - 0 -4000 man-days = 53,292. 09 man-days 
Group II - 5 000-7000 man-days — 28, 806.66 man-days 
Group III - 8000 - 13,000 man-days = 12, 049.49 man-days 
Group IV - 14,000 - 23,000 man-days = 43, 860. 31 man-days 
Total - 138,008.45 man-days 

4. Total man-hours - 690, 038. 04 man-hours 
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Cone lusions 

Table III shows that these sites fall into distinct groupings, and 
the chi-square tests indicate that these groupings are not random or 
by chance. The graphs indicate the inverse relationship between num¬ 
ber of sites and labor invested. What is important to note is that the 
amount of labor invested determines site size, the groupings of sites 
by similar amounts of labor invested, and the size of these groupings. 

One interpretation of these conclusions is that these large sites 
served some special purpose, such as ceremonial centers. The smal¬ 
ler sites may have been individual house mounds. Therefore, only a 
few of these large sites would have been constructed, requiring large 
numbers of man-days. The entire labor force may have been brought 
together for these community projects, while only small groups would 
have been needed for the individual projects. 

It may be found that the amount of labor invested in a site is de¬ 
termined by the function of the site, i.e. , whether it is a common 
structure or a religious or ceremonial structure. However, this paper 
deals only with large sites versus small sites, and their relationships 
to the amount of labor. It says nothing about the actual function of the 
sites. Further research is needed to test this conclusion. A study of 
the spatial distribution may give some clues as to whether the group¬ 
ings of the sites had any special significance or patterns. More ex¬ 
cavation would probably show the types of structures at each site. A 
study of the types of artifacts found at each kind of site would also give 
clues as to the function of each site. 
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Introduc tion 

A survey of archaeological sites employing systematic survey 
design and recording of data was performed during two summers in 
Guanajuato, Mexico. Seventy three sites were found. A preliminary 
analysis of 52 variables was recorded across the surface of each site. 
These data were tabulated (see Appendix A) and then a series of cal¬ 
culations were applied (see Appendix B). The tabulations and calcula¬ 
tions are the basis of the following proposals. 

Although it is not the purpose of this paper, the data contained in 
Appendix A may be used for comparisons with other aspects of the total 
project. These include ecology samples, geographic locations, seismo¬ 
graph data and others.^ 

From the surface survey data, I suggest four proposals about pre¬ 
historic activities which took place in Canada de Alfaro. They are 
based upon the relationship between ethnographic analogy and the actual 
distribution of the analytical results of the survey. 

Proposal 1 

Proposal 1 states that the domestic food preparation area for all 
73 sites exists on the East side of the valley. The sites are 15, 

16 and some surrounding sites. The most revealing data 

supporting this proposal is that these sites contain the only hearths 
in the entire region. A hearth is one of the most essential items for 
food preparation. A second factor is the relatively scarce amount of 
broken sherds. This suggests that since utensils are critical for the 
function of the area, people are more likely to be careful with them. 

This results in fewer broken ceramic food preparation utensils. The 
only hand tools found among the entire array of sites were found in this 
area. Hand tools are another essential item for food preparation. Water 
is permanently available because the river runs close, but not so close as 
to infest the area with mosquitos. Another intermittent stream runs 
close to the sites on the North side. Two more suggestions are 1) these 
sites are centrally located and easily accessible from all other sites; 

2) the hills, which border these sites on the North and East, protect 
the sites from adverse winds that could easily extinguish food prepara¬ 
tion fires. 


National Science Foundation Grant Proposal for Guanajuato 
Project, Ezra B. W. Zubrow, Stanford University, 1972. 

2 

For a map of the site area see Zubrow's introductory article. 
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Proposal 2 

Proposal 2 suggests that plaza sites are the trade centers for 
the site community. Two representative sites are 52 and 57 (see map). 
The primary reason supporting this proposal is that these sites lie close 
to the river, close to each other, and are located on trails. Rivers 
have a prominent function in trade and communication. In this case, 
the river is far too small to carry boats. However, its banks create 
the easiest areas in which to walk. 

The area which contains the plaza sites is separate from, but still 
within easy access to the area prescribed above as the food preparation 
area. The separate clustering of the two site types (plaza and food 
preparation) derives from the necessary separation of two different 
surpluses. One area deals with a surplus of merchandise for trade, 
while the other handles the surplus food for the community's immediate 
consumption. A mixing of the two could result in confusion. Outsiders 
could intrude on the community and the community could distract pro¬ 
spective traders from their objective. Buying and selling rarely occurs 
in the domestic area; even if a business is in a home, the two are 
separate. 


Proposal 3 

Proposal 3 suggests that the hilltop sites are the lithic workshops. 
The data show that these sites contain the highest concentration of stone 
points and chips. Since it is suggested that the trade and food prepara¬ 
tion areas are located in the protected valley, it seems probable that 
the industry would be separate from the two, due to lack of space and 
undesirability for industry to be close to the domestic area. 

There are several advantages to use hilltop sites as workshops. 

It is easy to defend possible secret methods of production or possible 
stockpiles of supply. It is possible to see (the data show that these 
sites have the best visibility) whether or not trade below in the valley 
is progressing well and thus continue to produce or halt production 
when necessary. 

There is one plaza site located in the hills which could indicate 
1) a trade center for prestigious customers or 2) an alternate trade 
center for when conditions are adverse in the valley (e. g. high mos¬ 
quito density). Assuming the development of this adversity, regard¬ 
less of its nature, it is more probable that an alternate trade center 
would be closer to its industry than to its domestic center. 

These sites are the closest spaced of all sites (high density) sug¬ 
gesting that distancing between sites, which may be hard to accomplish 











- 162 - 


on hill sites, was not necessary. Indeed, the most crowded areas 
today are the industrial centers and business districts. 


Proposal 4 

Proposal 4 suggests that valley circle and oval sites are housing 
sites. These sites account for all the remaining sites in the valley 
which are not already designated food preparation or trade center. 
They are spaced, within a given region, the farthest from each other 
(low density) of all site areas examined. This is significant when con¬ 
sidered against the contemporary trend that sought after housing has 
a space buffer zone but is still in contact with the community. The 
area,^ or square footage of these sites, revealed an average closely 
approximating the area of contemporary living units. Other sites are 
much larger or smaller. These sites are within easy access to the 
river, thus water is always available for the inhabitants. The most 
important factors supporting this proposal are 1) these sites lie in a 
protective valley, 2) are close to a permanent water supply, and 3) 
are spaced so as to maintain a buffer zone of privacy while remaining 
within easy access to food preparation, trade center, and industry 
sites. Finally, these proposed living spaces are on flat ground (0-10 
percent slope). 


Conclusions 

The primary criticism of the data in Appendix A is that it is 
incomplete for certain variables. This was a result of the exigencies 
of the field. 

In order to formally test proposals 1 through 4, it would be neces¬ 
sary to obtain an independent set of data. 


-^-— 

The area was not tabulated or calculated in this paper because 
this study has already been completed for the Guanajuato Project. 
See the paper in this volume, *•' The Relationship between Labor and 
the Size of Sites in Guanajuato, Mexico. 
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Appendix A 


Appendix A includes the list and tabulations of relevant variables. 
The tabulation charts exclude variables 9, 26, and 29-50 because 
they were not recorded enough to warrant inclusion in the charts. 4 
Variables 21-27 are tabulated in patterns corresponding to the form 
of the questions in the list of variables. 

Variables 21-24 consist of 1) five blade types, 2) four columns 
indicating the quantity of each blade type. On the tabulation chart, 
five spaces should be occupied within each box by the letters A, B, 

C, or D. Whenever a dash occurs in a space within a given box in the 
tabulation, it is because the variable representing the space was not 
recorded. The vertical arrangement of the letters in the box corre¬ 
sponds to the vertical arrangement of the five variables in the variable 
list. This arrangement holds for variables 25 and 27, (see tabulation 
chart), though the number of letters in the box changes to 3 and 4 
respectively. 

f Variables 6, 28 and 37 were types on separate sheets titled, 
Additional Data Sheets 11 and should be considered as part of the tabu¬ 
lation chart. 

Variable 51, in the tabulation chart, shows 5 columns for each 
site number. These columns correspond to the order of the obstacles 
as listed in the variable list; river, hill,mountain, rugged, and arroyo. 
The distances recorded in each box are given in yards unless otherwise 
indicated (all fractions are in terms of miles). Variable 52 in the 
tabulation chart has two columns per site. One shows the number of 
sites visible from the given site; the other shows a letter indicating 
the shape of the given site (R = round, O = oval, rec = rectangular, 
s q = square, and Plaza = plaza). 


This information is. included in the calculations (Appendix B) 
by site group. The information may be found in the original site forms 
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Variahles Relevant to the Tabulations and/or Calculations 


1. Nearest trail 


2 . 


3. 


4. 


5. 


Nearby buildings 


Nearest neighbor sites ( # and distance) 

Type of Holding a) Ejido Plot, b) Pequena Propiedad, c) Hacienda, 
d) House Plot, e) Other 

Possibility of Destruction a) Agriculture, e) Erosion, p) Pot hunting, 

o) Other 


7. 


Type of Site 


a) Rubble mound, b) Ramada, c) Ceremonial mounds, 

d) Sherd scatters, e) Activity Areas, f) Houses mounds, 
g) Chip stone scatters, h) Pit houses, i) Burials, j) Other 


8. Exposure a) Exposed, b) Partly exposed, c) Burried, d) Other 

10. Condition a) Preserved, b) Part-preserved, c) Rubble 


11. Dwellings (#) 

12. Construction a) Daub Ik Wattle, b) Rubble, c) Worked stone 

13. Sherds a) very abundant (20/sq. meter), b) moderate (10/sq. meter) 

c) sparse (5/sq. meter), d) very sparse (1/sq. meter), 

e) absent 

a) relatively uniform distribution, b) localized, c) variable, 
a) painted, b) textured, c) plain, d) figurines 

14. Stone a) very abundant (20/sq. meter), b) moderate (10/sq. meter), 

c) sparse (5 / sq. meter), d) very sparse ( 1 /sq. meter), 

e) absent 

a) relatively uniform distribution, b) localized, c) variable 

15. Perishables a) Textiles, b) Cord, c) Other, d) None 

16. Other evidence of construction a) cut stone, b) lajas, c) tepetate 

d) adobe, e) Other 

17. Hearths a) abundant (10+), b) moderate (5), c) few (1), d) none, 

e) exact # 

18. Trash a) % Rock, b) % Bone (animal), c) % Human, d) % Tools, 

e)% perishables, f) Other 


19. Petroglyphs # 

20. Tools a) Mano #, b) Metate #, c) Mortar #, 

f) Scrapers #, g) Other 


d) Pestle #, 


e) Awls #, 
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f 


a) Abundant b) Moderate c) Sparse d) None 
21. Obsidian Blades (50/sq. meter) (25/sq. meter) (10/sq. meter) 

" Scrapers " u ,, ,, 

" Cores " ii H ,, 

" Chips 11 ii n ,, 

" Flakes & " n ,, „ 

other debitage 


22. Chert Blades " 

" Scrapers " 

" Cores " 

" Chips 

" Flakes & " 

other debitage 


i r 
T I 
it 
n 
11 


11 

f 1 

If 

11 

T1 

r i 

II 

M 

M 

ii 


23. Slate Blades 

" Scrapers 
" Cores 
" Chips 
" Flakes & 
other debitage 

24. Basalt Blades 

" Scrapers 
" Cores 
" Chips 
" Flakes & 
other debitage 


25. Ceramics: 

Preclassic 

Classic 

Postclassic 

26. Trails 


a) 0-10% b) 10-25% c) 25-50% d) 50-75% 


if 

it 

r i 

it 

it 

n 


it 

ii 

11 

ii 

11 

it 


e) 75%+ 

I! 

ft 

If 


27. Agriculture 1) Evidence of a) absent, b) present 

2) Irrigation a) River, b) Runoff, c) Reservoir, d) None, 

e) Other 

3) Patterning a) Terracing, b) Waffling, c) Strips, 

d) None, e) Other 

4) % arable land a) 0-10%, b) 10-25%, c) 25-50%, 

d) 50-75%, e) 75%+ 

29. Soil Type a) Sand, b) Silt, c) Clay, d) Aeolian, e) Alluvium, f) Gravels, 
g) Screes 


31. Erosion a) Extreme, b) Moderate, c) Slight, d) None 


A 
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32. Water: Drainage a) Primary, b) Secondary, c) Tertiary 

Distance from site a) 100m, b) < 1 km, c) < 2 km, d) > 2 km 
Type a) Spring, b) Reservoir, c) River, d) Lake, e) Other 
Type a) Permanent, b) Intermittant 

33. Topography a) Valley, b) Valley side, c) Hill, d) Bluff, e) Mesa 

f) Mountain, g) Other 

Exposure a) North, b) South, c) East, d) West 

Slope a) 0-10°, b) 10-20°, c) 20-3 0°, d) 3 0-40°, e) >40° 

34. Geology a) Predominant rock type by % 

sandstone, shale,limestone, conglomerate, other 
sedimentary, slate, schistose, marble, quartzite, coal, 
other metamorphic, volcanic lava, basalt, granite, 
quartz, other igneous, no predominate outcropping 

36. Ecotype a) Desert 1) Cool, 2) Hot, 3) Succulent 

b) Chaparral 1) Pinon, 2) Juniper, 3) Aspen 

c) Grassland 

d) Pine forest 

e) Stream edge 

o/ 0 major type a) 0-10%, b) 20-25%, c) 25-50%, d) 50-75%, e) 75%+ 
Density of ground cover a) 0-10%, b) 10-25%, c) 25-50%, 

d) 50-75%, e) 75% + 

51. Primary Obstacles in and around a site 

1) River a) on, b) near (distance) 

2) Hill a) on, b) near (distance) 

3) Mountain a) on, b) near (distance) 

4) Rugged Terrain a) on, b) near (distance) 

5) Arroyo a) on, b) near (distance) 

52. Number of visible sites from any given site 


A 
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TABULATION Ot' PRELIMINARY ANALYSIS 
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Additional Data Sheet for Variable 6 

Sites 

1) xxx 

2) Irregular mound indented in the middle of a large mound. Mesquite 
and cactus on it. Used by farmers for shade. 

3) In cornfield with tree — depression on second small mound. 

4) xxx 

5) Very small low mound with undifferentiated sides and top. Increase in 
the number of lithics and ceramics in the middle of the field near an 
oak tree. 

6) xxx 

7) xxx 

8) Getting into a cluster of sites. Surrounding field richer as work moves 
toward the river. 

9) xxx 

10) xxx 

11) xxx 

12) xxx 

13) Partially destroyed by path and wall across it. 

14) xxx 

15) Discrepancy as to whether there are one or two sites due to apparent 
wash area in between top ones. 

16) There is a terrace not included or part of the site, but which is drawn. 
Altar 1 near the river. 

17) Very near dry river ravine. Cornfield with contour and high density 
(material) on NE side. 

18) Cornfield with no noticeable slope. 

19) Earth has washed away a higher point exposing rock. 

20) Top of hill. The top has illegal excavation and the material from 
inside is scattered on top. 

20 Road) This road runs down the hill towards the site (from site 20) 

850' long and 5'6" wide. 

20 Plaza) Plaza is west of mound. 
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21 Ballcourt) Possible Ballcourt 60 1 x 12 T . 

21) xxx 

21 Road) 549' long, downhill from site 21 same width as 20 road. 

22) Holes dug in top. 

23) xxx 

24) xxx 

25) See mound cut. 

26) Secondary mounding. 

27) xxx 

28) Side 1-5 are questionable since they fall into the 27/28 border. 

29) xxx 
3 0) xxx 

31) Part which has been undisturbed by plowing has been covered with 
new rocks. 

32) Middle and S side will probably continue to be covered by new rock. 

33) In middle of fallow field. Very rocky with clay or lava soil. Very 
exposed, no mounding. Both sherd and ceramic concentration. 

34) In field across arroyo from 33. 

35) Adjacent to 34 and 33. A small field above field which appears to 
be workshop area. 

36) In middle of fallow field. Obvious mound with both sherds and lithics. 

37) Terraced cornfield. Cactus shrub to W. 

3 8) No mounds or abundance of ceramics or lithics — only 3 samples 
taken. 

3 9) In valley E of shoulder of full hill. 

40) xxx 

41) Stone walled well 60 1 deep in center. 

42) xxx 

43) xxx 

44) Large lithic concentration. No mounding, therefore no samples 
for the top and side. 

45) xxx 
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46) xxx 

47) xxx 

48) Circular. 

49) Only sherd scatter. 

50) Small site, almost part of 51. 

51) Seems to be along terrace edge along with several other sites. 

52) Mound evident from field tq south and no mounding on plane plaza. 

53) On terrace w/ 51, 52. 

53p) xxx 

54) Modern rock cover make sampling impossible. 

55) This site is right next to excavation site 43 to W. 

56) 3 00 yds W of site 43. 

57) xxx 

57 Plaza) 5 samples taken from plaza which lies between 57 and 58. 

58) xxx 

59) One large pothole, most samples taken on border. 

60) xxx 

61) xxx 

62) Sherd scatter - 5 samples. 

63) xxx 

64) Treated as sherd scatter due to modern rocks. 

65) xxx 

66) xxx 

67) xxx 

68) xxx 

69) Large sherd scatter to S not included in claculations or collections. 
71 ) - 73B) xxx 
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Additinnal Data Sheet for Variable 28_ 

Sites 

1) xxx 

2) xxx 

3) xxx 

4) More lithics than pots on North side, more pots on S. 

5) xxx 

6) xxx 

7) Debitage and pottery piled in one spot. Definite disturbance. 

8) Middle of site is being excavated. 

9) xxx 

10) xxx 

11) xxx 

12) xxx 

13) xxx 

15) xxx 

16 ) xxx 
Altar 1) xxx 

17 ) - 20 Road) xxx 

20 Plaza) fields are on lower slope of hill. 

21 B ale our t) xxx 

22 Road) xxx 
23) -27) xxx 

28) Probable human long bone in side 5. 

29) -31) xxx 

32) Extremely rocky field 

33) Site undifferentiated from fields where it is located. 

34) In fallow field sandy soil w/humus. 

35) Field surrounding mound is very rocky, hasn't been cleared recen y 

36) Fields fallow and used for grazing. 
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37) Site is terraced cornfield. 

38) -42) xxx 

) Previous cultivation in plaza, now used for grazing. 

44) Isolated by heavy mesquite, cactus, and brush. Borders on cornfield. 

45) -48) xxx 

49) Sherd scatter in cornfield 

50) -51) xxx 

52) Top of mound is part of adjacent plaza - wall of terrace falls 


toward the river. 





33) Field used for grazing, remains of irrigation ditch. 

53 Plaza) Previously cultivated and irrigated, now follow and grazed. 
54)-57 Plaza) xxx 
58)-73B) xxx 
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A-dditional Data Sheet for Variable 37 


Sites 

1) - 13) xxx 

15) 8 mesquite, 1 napel cactus, 2 canoe plus corn 

16) Agricultural field. 

17) -28) xxx 

29) River 50 yards. 

30) River 100 yards. 

31) -45) xxx 

46) Water 200'. 

47) -49) xxx 

50) River less than 100 yds. 

51) River 100 yds. 

52) On river. 

53) River yields gravel 500 yds. 

53 Plaza) River for irrigation close. 


54)-73B) xxx 
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Appendix B 

Calculations of Tabulations 

There are three methods used for calculating the tabulations: 

1) a graph which depicted individual site characteristics, 2) a simple 
chart totaling the characteristics for all sites or for a given site group,^ 
3) in three cases, variables 34, 51, and 52, a chart was made. A high 
and a low point were ascertained. For the calculation of the middle point 
and a description of the high and low .points for each variable, see below. 

Variable 34 asks for the compositions of rock types for each site. 

In each case, the recorder marked all of the rock types found at each 
site and of those rocks marked, the percentage that each contributes 
to the total composition of.each site. A chart was made to determine 
the total number of appearances each rock made per site group. Then 
these totals, for each rock, were combined to determine the total for 
all sites. Next charted was 1) the highest percentage of each rock type 
for all 73 sites, 2) the lowest percentage of each rock type for all 73 
sites, and 3) the average of all percentages recorded for each rock type. 

Variable 54 asks for the distances of each site from several ob¬ 
stacles (e.g. river, hill, mountain). For each site, it was recorded 
whether the obstacle was "on" or "near" a site. A chart was made to 
determine the number of "on" and "near" sites, per site group, per 
obstacle. Next, a chart was made which showed per site group for 
each obstacle, the lowest recorded distance, the highest recorded dis¬ 
tance, and the distance which occurred most of the points in between. 

This results in three trends: 1) the highest recorded distance for each 
obstacle (this is the distance from a given obstacle to the furthest site 
within a given site group), 2) the lowest recorded distance for each 
obstacle, and 3) all the points in between occurring 4 or more times. 

If, of all the middle points, a given distance didn't occur 4 or more 
times and the recordings were scattered, no middle trend was given. 

It was decided that a trend could only be valid if it contained 4 or more 
of the same recorded distance. The assumption of validity is that if a 
number occurs in more than fifty percent of the cases examined, then 
a trend exists. If each site group contains 15 sites and four to five of 
them are included in determining the high and low trends, then 4 or 
more is roughly 50 percent of the remaining sites. 


There are five site groups, none of which contain more than 18 
sites. The grouping of sites was arbitrary. Group 1) included sites 

•t?; 4, 2) sites 15 " 22 ’ 3 ) sites 23-35, 4) sites 36-57p, 5) sites 58-73B. 
The site-composition of each site group varies for certain variables, 
see calculation charts. 
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Variable 1 
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Variable 1 continued 


Distance of nearest trail in yards 
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Variable 2 


Distance of nearest building in yards or miles 
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Variable 2 continued 


Distance of nearest building in yards or miles 
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Variable 3 
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Variable 3 continued 


No. of site 
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nearest site 
is recorded 
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Calculation Charts 


Variable 4: 100% C 


Variable 5: 
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55-73B 
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18 

1 

0 





Total (all sites) 
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Variable 7: 

A 

B 

C 

D 

E 

F 

G 

J 

Site groups: 1-14 
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Variable 8: 

A 
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Variable 10: 

A_ 

B 

C 






Site groups: 1-14 
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Calculation charts continued 


Variable 13: 


A 

B 

C 

D 

A 

B 

C 

A 

B 

C 

D 

Site groups: 

1-14 


4 

7 

5 

1 
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5 
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15-21 
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3 



22-35 
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6 



36-54 
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55-73B 
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12 
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Totals 
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11 

14 
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Variable 14: 


A 
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D 

E 



A 

B_ 

C 



Site groups: 

1-14 
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Variable 15: 100% D (no perishables) 


Variable 16: wall on sites: 2, 
drain on site: 15 
road on site: 20 

8, 
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29, 

30, 36, 43 
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Variable 18: A 

B 
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2 
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Variable 19: 1 piece of negative painted pottery 

in site group 1 




















- 186 - 


Calculation charts continued 


no. of appearances amount at each appear. 


Variable 20: 

Site groups: 1-14 
15-21 
22-35 
36-54 
55-73B 

Totals 


A B C D E F G 

9 1 11 

2 2 
12 11 1 
6 1 113 

5 1 1 

23 4 1 3 1 4 3 


A B, _C JD _E F G_ 
17 1 11 

3 2 
2 2 11 1 

11 1 1 1 12 

5 1 1 

38 4 1 3 1 14 3 


Variable 21: 

A 
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Variable 22: 

A_ 
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Variable 23; 

A_ 

B_ 

C_ 

D_ 
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Variable 24: 
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Calculation charts continued 


Variable 25: 

A 

B_ 

C_ 

D 

E_ 
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Variable 26: Only four trails 

= ir 

rigation and grazing, 2 unknown 

Variable 27: 
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Variable 29: 








Soil types 

A 
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Variable 31: 








Erosion 
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Variable 32: 

A 
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water drainage 
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Calculation charts continued 



* The ecotype for all sites is succulent desert 
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Calculation charts continued 

Variable 51b: 


Site groups: 



i - 

Lb - 
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0 
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20 

25 0 


60 


mid 

arroyo 

60 





high 


Variable 52: 
Site groups: 

1-14 

15-22 

23-35 

36-57p 

58-73B 


low 

0 

0 

0 

0 

0 


mid 

4 

2 

3 

0 

1 


high 

5 

3 

9 

9 

5 
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Theoretical Background 

In this section, the theoretical basis for the model proposed by 
this paper is explained and explored. Although the description of 
Schumpeter's model is incomplete, the fundamental points are pre¬ 
sented. ^ 

It should be kept in mind that several different models are dis¬ 
cussed in this paper. First, Schumpeter's model, then an adaption of 
it, and finally and explicit simulation are discussed. 

Schnmgeter's Model of Unstable Economic Growth 

Schumpeter's model of economic development is a model which 
views the process of development as driven by innovation, with a re¬ 
sultant growth pattern which is not uniform. The innovation-expansion- 
(growth)-equilibrium cycle operates as follows: 

The model starts with many firms in perfect competition, with the 
economy in equilibrium. Equilibrium, the stable condition, here means 
primarily that profits are zero. Furthermore, it means that the mar¬ 
ginal utility of additional investment, with the current technology, is 
also zero. 

Now, for reasons presented later, let us assume that a techno¬ 
logical innovation is introduced. Acceptance of the innovation will re¬ 
quire construction, since the innovation requires new investment for 
implementation. New construction will be undertaken: first, because 
of expectation of monopoly profits in the utilization of the innovation, 
and because of increased marginal efficiency of capital; second, be¬ 
cause old production centers will become obsolete; and third, because 
the innovation may cause an increased propensity to consume, so con¬ 
sumers will dip into their savings, and demand will increase. 

Schumpeter claims that after the first firms utilize the innovation, 
the remainder of the firms will follow in clusters. The industry will 
gradually retool, with a resultant higher investment level. As the 
relevant industry and related industries near total acceptance of the 
innovation, the rate of new investment will decrease to a maintenance 
level. As the supply increases, prices fall, debts are paid off and a 
new equilibrium is reached at a new technology/investment level. 

Schumpeter further argues that innovations will occur in stable 
conditions for two reasons. First, the pressure to innovate is largest 
at equilibrium, (perhaps because people's desire to engage in additional 
activity to produce growth is greatest when no activity to this end is 
being executed or the role of science at equilibrium differs from its 

^Information used in this paper concerning Schumpeter's model 
was obtained mainly from Higgins, 1968. 


A 









role in a time of crisis or instability). Second, the risk of innovating 
is at a minimum (as the marginal efficiency of capital is nil). 

fdi.ggi.ns (1968) presents a set of pseudoformal relations concern¬ 
ing variables in the model. These relations are summarized in Fig. 1. 
In this diagram, an arrow indicates an unspecified functional depen¬ 
dency. (a) (b) indicates that "b" is dependent on "a". 

Reformulation of Schumpeter's Model 

One thesis of this paper is that Schumpeter's unstable growth model 
can be reformulated to apply to subsistence, non-market economies in 
a meaningful way (Zubrow, 1973a). 

First, as our unit of analysis, we will use the subsistence unit, 
and its environmental resource base. As our analog to the "perfect 
competition" assumption, we must assume that the system is closed. 

In a market economy, money is a convenient measure of such ent¬ 
ities as output, investment, and debt. In a non-market economy a dif¬ 
ferent unit of measurement must be found to evaluate these entities. 

We have chosen the concept of energy to accomplish this (Clark & 
Haswell, 1970). 

Output, or production, is the quantity of consumeable energy pro¬ 
duced or gathered in some time period. If 1000kg of wheat were pro¬ 
duced in a year, this would represent an output of 3440 Kcal/year, 
as wheat contains roughly 3.440 Kcal per kg (Flannery, 1971). Growth, 
then, is a simple increase in output energy of the system under consi- 
deration. 

For the purposes of this paper, investment is a quantity of human 
energy which has been expended to influence the capacity or produc¬ 
tivity of an output unit. Therefore, cumulative labor going into the 
construction of an irrigation system, is considered to be investment, 
while energy going into annual production is not. 

Savings for Schumpeter is different than investment. Savings will 
be considered to be a quantity of stored consumeable energy, that is, 
energy not tied up in investment. Schumpeter considers an innovation 
to be technological progress or resource discovery. Technology is 
here defined as the ratio of output energy to input energy of some pro¬ 
ductive system (e.g. if X Kcals are expended in hunting and Y Kcals 
of energy are obtained in meat, then the technology evidenced is Y/X, 
which is simply an efficiency). Hence, a technological innovation 
increases the output/input ratio. 

The resource base is viewed as having an absolute limit of poten¬ 
tial energy producing capacity. Of that, some fraction is considered 
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to be utilizable in the sense that it is known, and can be used at will. 
Therefore a resource discovery has the effect of increasing the frac¬ 
tion of energy potential which is usable. 

To start the model, we have a subsistence unit in equilibrium 
with its environment. All production is consumed, yet no one starves. 
There is no advantage at this time to increased investment, as the 
current technology is fully realized for the labor force. 

Now, we propose, an innovation occurs. We assume that the inno¬ 
vation will require an investment in order to be realized. The subsist¬ 
ence group will dip into existing savings to finance the investment m 
expectation of higher production. 

As the fraction of the labor force using the innovation produces a 
surplus (relative to the old technology), this surplus may be invested 
in the technology until the technology is fully implemented. 

If a group strives to maintain a certain level of savings for use in 
times of adversity or in readiness for an innovation, the lack of this 
level of savings may be considered a debt. The magnitude of the debt, 
then, is the difference between the level of savings desired by the group 
and the current level of savings. 

Production in excess of that required for the subsistence of the 
group and the maintenance of the desired savings may be considered 
surplus. ("Surplus" is used in a slightly different sense m the simu¬ 
lation. ) 

As a technology nears full implementation, there will be production 
which is not directly consumed, and not invested. This will go into re¬ 
payment of the debt or into additional savings. 

If the population increases to consume this surplus, or if the 
surplus is otherwise consumed, the unit will eventually reach an equi¬ 
librium state,like that initially specified. 

This process description is summarized in Fig. 2. 
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Figure 2. Reformulated Growth Model 
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Heuristic Example 

To enable the reader to get a better feeling for this model, a 
simplistic example of its operation in a hypothetical subsistence unit 
is presented below* 

Let us suppose that we have an island populated by one subsist¬ 
ence unit in one ecological zone. This zone is being entirely utilized 
in the production of .one food crop at the highest level of technology 
available. This situation has resulted in a stable population. 

An inventive native hypothesizes that irrigating the fields would 
increase crop production (i.e. resource utilization). He digs his ditches 
and produces crops from his irrigated fields. To the amazement of 
his neighbors, his production is 25 percent better per unit area. 

Gradually others decide to implement this innovation. For some 
this requires going into debt, since in order to use this innovation, 
one must invest significant time in digging irrigation ditches before 
any production is realized. On the individual level, the debt may be 
in the form of an obligation, to those who helped him in his construc¬ 
tion. After a given individual finishes his initial investment, and crops 
are produced, he must repay his debts with his excess production. His 
repayment enables others to invest in the innovation. 

On a less direct level, other segments of the social structure must 
adapt to irrigation. A selective set of craftsmen, specializing in sluice 
gates, may appear, or roads may require rerouting, around irrigation 
ditches. 

Eventually most producers will adopt irrigation, and their debts 
will be paid off. At this point we are at a new equilibrium, with a 
higher output, and a higher level of investment. 


The Simulation Model 

In this section a simulation of the reformulation of Schumpeter’s 
model is described. First there is a general discussion of the nature 
of the process of the transformation of Schumpeter’s model, through an 
intermediate step, to a final simulation and test. Then the operation 
of the model is discussed, along with the assumptions of the model. 
Finally, the behavior of the simulation model is described. 

Thus far two "models’ 1 have been discussed, Schumpeter's and its 
reformulation. Schumpeter’s model is highly intuitive but utterly un- 
testable. One can conceive of testable implications of Schumpter’s 
model (Zubrow, 1973a: 11-29 ff). However, even to test these 


E. g. innovations will occur in times of economic stability. 
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implications one must make a large number of assumptions. In addi¬ 
tion it is clear that a test of these implications cannot validate Schum¬ 
peter, since the operation of the model is not adequately specified. ^ 

The second model presented above was an attempt to partially 
remedy the under specification and to introduce the reader to concepts 
of the simulation. This second model was elaborated to form the final 
simulation model. The simulation model is formally described in 
Appendix A in an algorithmic language, ALGOL W. The ALGOL 
de5criptioa is, in dlect, ct dcociiptiun of the model to an ignorant 

computer. This language has the advantages of logical precision and 

reasonable readability. 

In addition to being a formal description of the model it is a com¬ 
puter program for the simulation model. This has two immediate con¬ 
sequences. First, if the model has inadequate assumptions built in, 
the simulation will not run. Put in another way, we know the output 
produced is exactly a product of the model including its assumptions 
(Dutton & Starbuck, 1971: 4). Secondly, we can run the simulation and 
see how the model behaves with different inputs. 

Operation 

This description of the operation of the model will, as much as 
possible, follow the structure of the formal description given in Appendix 
A, and use similar terminology. The description presented below dem¬ 
onstrates the central features and assumptions of the model. However, 
it neglects details. The ALGOL description was written with read¬ 
ability in mind and should be consulted for detailed information. 

(Names in capital letters refer either to variable names or pro¬ 
cedure (subroutine) names in the computer model. ) 

INITIALIZE. Assume the initial state of the economy to be stable 
and initialize variables accordingly. This means that as long as tech¬ 
nology remains the same, production, population, savings and invest¬ 
ment should all stay constant. Several initial variable conditions are 
input by the user (e. g. RESOURCE^JPOTENTIAL — absolute energy 
potential of the resource base, HUMANJSUBSISTENCE ENERGY - 
energy required for individual human subsistence). Others are 


That is, Schumpeter deals with a large number of variables, let 
us say xp . . , x n . To say that xj = f (x^-x^ . . . , x m ); x^ is a func¬ 
tion of x^-, . . . , x m is not adequate. To deal with a model we must know 
the operation of f, the function. Even to say that x^ is an increasing 
function of xj is not sufficient. If we know x^ = f(x^); x^ = f (xq); X£ incr => 
X 1 incr; xq incr X 3 deer, we may want to know what happens to x^ + X 3 
if X£ and xq are increased by one unit. 
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computed (e. g. TECHNOLOGY — output/input ratio, RESOURCE 
CAPACITY - energy of resource base available at the current technol¬ 
ogy), The result of this is that an equilibrium need not be calculated 
by the user of the model, but results from the input and calculated 
variable values. 

PRODUCE. Using LABOR ENERGY — the energy expended by 
the population for production, the INVESTMENT^ LEVEL in each 
technology level, the RESOURCEjCAPACITY, and the TECHNOLOGY 
available, compute the OUTPUT for this time period. The INVEST¬ 
MENT __LEVEL (TECH__LEVEL) is the fraction of the RESOURCE 
CAPACITY which has been improved up to the technology level indicated 
by TECH__LEVEL. That is, an innovation will at times be incompletely 
implemented, because the necessary investment has not been made to 
utilize the innovation. An innovation produces a new maximum 
TECH__LEVEL. 

Computing the OUTPUT happens as follows: at the top level of 
technology available (CUR_TECH_LEVEL) the production at this tech¬ 
nological level equals LABOR (TECH_LEVEL)*TECHNOLOGY. LABOR 
represents the amount of labor energy available which can be used at 
this TECH_LEVEL to utilize all of the RESOURCE_CAPACITY which 
has been improved to this TECH LEVEL. 

At the next lower TECH_LEVEL the same procedure is executed, 
utilizing only the remaining RESOURCE_CAPACITY and LABOR_ENERGY. 

OUTPUT is then the sum of production at each TECH__LEVEL going 
down until all labor energy is utilized, or all the RESOURCE CAPACITY 
used. 

INNOVATE. Production in the last two time periods is examined. 

If there has been less than a set percentage growth (STABILITY) then 
with a set probability (PROPENSITY_TO_INNOVATE) an innovation 
will occur. If the subsistence unit is completely utilizing the current 
RESOURCEJCAPACITY, then the innovation will be in the form of a 
RESOURCE_EXPANSION, otherwise it will be a TECHNOLOGICAL 
INNOVATION (see the section on the reformulation of Schumpeter ! s 
model for explanation). This last assumption is made since, in each 
case, the innovation the unit receives is the only one it can use to in¬ 
crease production.^ Either process results in an increase of the CUR_ 
TECHJLEVEL. 

^ If the subsistence unit is operating at the resource capacity, then 
all of the utilizable energy is being drawn from the environment, and 
although an increase of efficiency would allow a decreased labor force 
it could not increase production. Likewise, an increase in capacity 
alone cannot increase production unless production already utilizes the 
entire resource capacity. 







- 200 - 


TECHNOLOGICALJLNNOVATION. A technological innovation is 
an increase of TECHNOLOGY, the output/input ratio. This increase 
is for lack of anything better a RANDOM fraction of a maximum in¬ 
crease in technology (MAX_TECH_INNOV) . 

RESOURCE EXPANSION. A resource expansion is an increase in 
RESOURCE UTLLIZ, the usable fraction of the RESOURCE_POTEN- 
TIAL. This”increase is a RANDOM fraction of a maximum increase 
(MAX RES EXPANSION). For obvious reasons RESOURCE_UTILIZ 

oannot gkgq©c1 1_ 

As an important assumption, both of the above innovation processes 

will require an increased level of investment per output unit (INVEST_ 

PER RES). To explain this, a digression to further explain investment 
is required. 

There is an initial investment per output unit required for produc¬ 
tion. For example if the base investment is 1, this implies that for each 
unit of output the resource base supplies, there must have been one unit 
of labor invested at some time. An increase in technology requires an 
added investment to the base investment. For example, if technology 
level 1 has a required investment of 1, level 2 may require an addi¬ 
tional *1 and level 3 an additional .15. Hence, to make unimproved 
land useful at the third technology level, 1 + . 1 + . 15 = 1. 25 units per 
output unit investment are required. 

For lack of a better procedure, the size of the increment to the 
investment is a function of the percent increase in TECHNOLOGY or 
in RESOURCE UTILIZ and a constant factor, CAPITAL_INCR_FACTOR. 
This factor relates the strength of the innovation to incremental invest¬ 
ment. (We are using CAPITAL_INCR_FACTOR = 1.5). 

ADJUST_POPULATION controls BIRTHS, DEATHS, and STAR¬ 
VATIONS. BIRTHRATE is determined differently depending on whether 
GROWTH exceeds STAGNANT. If the economy is growing, the BIRTH¬ 
RATE approaches a limit. The change in BIRTHRATE is based on the 
amount of SURPLUS in this time period. (SURPLUS = OUTPUT - 
NECESSARY_CONSUMPTION) (NECESSARY_C'.ONSUMPTION is the 
energy required to feed the POPULATION). 

If the economy is STAGNANT (with growth less than a fixed per¬ 
centage) then the BIRTHRATE depends in the same way upon SAVINGS 
rather than SURPLUS. In this case POPULATION tends to a size that 
allows a minimum SAVINGS_LIMIT to be maintained, (see reformula¬ 
tion of Schumpeter's model) 

DEATHRATE in this model is constant. STARVATIONS are in addi¬ 
tion to DEATHS produced by the DEATHRATE. If in a time period, 
NECESSARY_CONSUMPTION cannot be met from OUTPUT + SAVINGS, 
then the number of people who cannot be fed, starve. The model assumes 
that children starve first. 


i 
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This routine also adjusts the LABOR POOL. Individuals enter 
the LABORPOOL at a fixed MATURE AGE and leave it at OLD AGE. 

CONSUMEJENVEST. OUTPUT in this model goes 4 places. First 
OUTPUT is used to feed the POPULATION. Next the remainder of 
OUTPUT and accumulated SAVINGS are invested, so far as that is 
useful. The remainder, up to a limit (MAXjSAVINGS)^ is saved and 
the rest is somehow consumed or wasted. 

The amount going into INVESTMENT is that amount which, when 
made, allows the maximum proportion of the LABOR__POOL to use 
the TECHNOLOGY and RESOURCEJJTILIZ of the CUR_TECH_LEVEL. 
As explained above, in order to do this, INVESTMENT must be made 
at all lower levels, to be used at the top level. Hence INVESTMENT 
is made in order from the bottom level up to CUR__TECH__LEVEL. 

Behavior 

The simulation described in the section on operation (and Appendix 
A) has been run a number of times. In the remainder of this subsec¬ 
tion, the general behavior of the model will be discussed, and the reac¬ 
tion of the model to changes in different input parameters examined. 

The simulation has a large number of input parameters (20). Tests 
to determine behavior were done by taking a standard set of initializa- 
tions, and varying one parameter and comparing outputs. 

In preparation for a test of the model, research was done to deter¬ 
mine the input parameters for the test case at hand (Hay Hollow Valley, 
Arizona 100-13 00 A.D. ). The values used as standard input parameters 
came from the data discovered in preparing the test. (A discussion of 
this data is presented in the Data ). These standard input parameters 
are given along with a short description on pages 6-8 of Appendix A. 
Note that the time interval chosen represents one year. 

Hence the tests run were for a range of variables subjectively 
considered to be reasonable. 

Six variables were considered central in the analysis of this simu¬ 
lation. They are output, population, labor pool, starvations, invest¬ 
ment, and savings. A typical set of results is given in Appendix B. 
(Beware of the scale on the graphed output, similar graphs may differ 
in scale). In discussing variables affecting other variables, one must 
be aware of the fact that very many of them are interdependent. 


c 

This is done on the assumption that energy stored as food, decays 
over time, and can only be saved in limited quantities. 
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Output is characterized by cycles of stability followed by approx 
imately exponential growth leveling off at a higher stable output level. 
The growth is a response to an innovation. 

The time lag required for the innovation to significantly afiect the 
output, is affected by the amount of investment required for the innova¬ 
tion to be used, (small investment rapid response), the savings on 
hand which can go into investment (large saving rapid response) and 
the magnitude of the innovation (large innovation => rapid response) 
(however 3 large innovation => investment). 


The rate of growth is affected by a number of factors. An increase 
in maximum birthrate will increase growth rate, however, this can 
cause later instability in the output; a decrease in the maximum sav¬ 
ings which can be maintained, seems to increase the growth rate 
(probably by way of interaction with the population adjustor). Increasing 
the length of time workers stay in the work force also increases the 
growth rate. The factor by which investment is increased with a tech¬ 
nology increase has very little effect on the growth rate. Changing the 
initial investment level or minimum savings level has a similarly small 
effect. As would be expected, the greater the innovation, the larger the 

growth rate. 


Growth levels off to a new level of stability either when the output 
reaches the resource capacity, or when all available labor used at the 
top level of technology can just supply the population with subsisten 
energy (whichever occurs first). 


An elaboration of this point may be in order here. In the initial 
state, a fixed proportion of the population, the labor force, is exactly 
able to produce for the entire population. What a technological inno¬ 
vation does, is provide temporary overhead necessary to expand pro¬ 
duction. That is, in the initial state, any birth, over and deaths, 

would starve. Since new births do not enter the labor pool for perhaps 
10 years, they do not increase the production, only the consumption. 
However, by hypothesis, in the initial.stable state all production 
consumed by the population and there is no surplus. 

Population, in the first analysis, grows as output grows, with a 
significant time lag. (In these experiments 12-43 years), bu * 
of leveling off, settles into a cyclic unstable equilibrium.. Like outpu , 
population exhibits an exponential growth. The cyclic variation is a 
product of delay in entering the labor pool of new births and the co 
stant interaction between birthrate and output. This eye ic p enomen 
was not forseen. With variable combinations that make for unstable 
output, the variation in population may be so great that a stable popu¬ 
lation level is found well below a peak population carrying capacity. 

Direct population variables, e.g. birthrate, deathrate, have the 
expected effects on population. By the nature of the population adjust¬ 
ment function any increased output growth rate will result in increase 
population growth. 
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The labor pool grows much as population does, with the time lag 
introduced by the fact that individuals cannot enter the labor pool im¬ 
mediately after birth. What is perhaps surprising, is that even though 
population may cyclically vary 15 - 20 percent, the labor pool may 
remain constant, or may vary even more while maintaining a constant 
output. (The latter implies surpuls labor). 

Starvations seem mainly affected by the rapidity of change in 
population as discussed above. Because of this, increasing the sav¬ 
ings limit may increase starvations in the long run, since savings 
influence birthrate, and savings will cycle with excess production. 

Investment varies much as expected. It does not tend to cycle 
and always reaches 0 sufficiently long after an innovation, since there 
is a limited investment to be made at each technology increase. As 
mentioned above, increases in the factor of capital investment for re¬ 
source expansion and technological innovation only serve to retard 
production somewhat, and raise the level of investment at each time 
period. 

Savings, in general, grows smoothly along with the last leg of out¬ 
put growth. As population cycles, savings cycles, offset by 180 degrees. 
That is, low populations imply high savings. As population increases,^ 
it uses the savings. The savings drop to zero and some of the population 
starves. The birthrate is meanwhile decreasing and savings again grow. 

Most of the surprising action and instability occur by way of the 
constant delayed feedback between population and output. These fluc¬ 
tuations represent serious problems in the proposed model, as the 
degree of instability produced seems intuitively quite unlikely. 

Validation of the Simulation Model 

The simulation model presented, claims to explain and predict 
long term economic growth, mainly by reference to technology. The 
purpose of a validation is to determine how good a predictive tool the 
model is. In this section we do not attempt to present conclusive evi¬ 
dence. Rather, this section attempts to provide a view of what a vali¬ 
dation of the model would look like, and secondly to present what data 
has been collected which bears on the question of validity. 


Validation Procedure 

One strategy for validating a simulation model is to use it as a 
black box to predict known data, given other known, related data. Thus, 
if we can tell the simulation, the situation at time t, and it can reliably 
predict the state of its universe in later time periods without further 
input, then it has succeeded. In a given empirical situation the simu¬ 
lation needs certain information as input. If the model is true, it shou 
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gene rate further information concerning the situation at hand which 
fits the data. It should be noted that there may be any number of dif¬ 
ferent black boxes which will produce the same output from the input 
given 0 If however, one uses the black box on a variety of empirical^ 
situations and receives good fit, one gains confidence in the model. 

A better strategy, which is really only an elaboration of the first, 
would be to dissect tlie model and test its components, then reassemble 
them and test the unit. This strategy then puts the blackbox problem at 
the component level, rather than at the complete simulation level. For 
the simulation model of this paper, one might then wnnt to take the pop* 
ulation adjustment function and deal with it separately. It could be tested 
and improved. Then the other segments could be tested and the improved 
model (the combination of the improved segments) could be tested. If 
the fit with data was not good, and one had confidence in the subpro¬ 
cesses, one might then know to look for problems in the relations be¬ 
tween the subprocess, not in the processes themselves. In the first 
strategy one may have no idea where to look to solve problems that 
arise. 

Nevertheless, presented below is an instance of the first strategy. 

Data 

The simulation model, which is the subject of this paper, takes as 
input several initial conditions concerning physical/biological constants 
and social structural variables, generates random innovations and pro¬ 
duces, among other things, a prediction of population and output. 

In this section, we wish to take an empirical example, find the 
inputs and compare the predictions with the actual data. However, in 
addition to the required inputs we wish to provide the simulation with 
information concerning actual innovations. 6 * 8 


6 There exist statistical techniques to quantify some of these 
judgments, but they will not be discussed here. 

^ One can, of course, separate the model into any number of sub- 
assemblies and components, and test at each level. However, with 
limited time, knowledge, and resources this repeated subdivision rap¬ 
idly becomes infeasible. 

8 

In this simulation we do not claim to have any model of the genera¬ 
tion of innovation. It is characterized as a random process, for lack of 
any better alternative. Even if innovation is in some sense random, there 
is no reason to believe that our random number generator is the same as 
nature ! s (God l s). 

0 
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As the reader may have guessed, the data required to test this 
model over a long time period, say 1000 years, is not readily avail¬ 
able. The best job possible has been done within time constraints to 
collect as good data as possible for this validation. 

The empirical situation chosen is the Hay Hollow Valley, Arizona, 
from A.D. 100 to A.D. 1300. This area has been extensively studied 
by archaeologists directed by Dr. Paul S. Martin of the Field Museum 
of Natural History. “ In the time period under consideration, it appears 
that the population grew from about 100 to about 2500. In this perio 
the economy went from primary dependence on hunting and gathering to 
an irrigated agricultural system. The area considered is roughly 18.5 
square miles. 

Below are presented the data collected to initiate this validation. 


1. Environmental Variables 

a Resource potential. This may be thought of as the maximum 
energy production possible for the area. The resource potential is the 
sum of the resource potential in each of six ecological zones. The 
resource potential of an ecological zone is defined to be the product of 
the annual productivity per unit area, the energy per unit biomass, an 
the area. Productivity and energy per unit biomass come from Zubrow 
(1971b: 200-201). Area figures are also from Zubrow (1971a: 133). 


Resource 

Zone 

1 

2 

3 

4 

5 
7 

Totals 


Area 

mi 2 

Productivity 
gm biomass/mo 

Kcal/ 

gm biomass 

Months 

Potential 
10 6 cal 

1.17 

10. 08 

4 

4 

5.25*10 5 

£ 

2.15 

12.4 

4 

4 

1. 11 *10 

. 92 

14. 66 

4 

4 

5.5 0 *10 

c 

1.43 

7. 18 

4 

4 

4.28*10 

c 

2. 84 

2. 14 

4 

4 

2.49*10 

6 

9. 92 

18.43 

22. 70 

4 

4 

9. 93*10 

6 

12. 20*10 


Figure 3. Resource Productivity of the Hay Hollow Valley 


2. Population Data 

a Population . (There are a number of problems with this data 
discussed by Schiffer (1968:4-21) and Zubrow (1971b: 150-151.) Population 
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at a given time is estimated by the total number of habitation rooms 
used in that period times an estimate of the number of people/habitation 
room. The valley has been surveyed several times and an estimate of 
habitation rooms per time period is given by Zubrow (1973b). (The sites 
are mainly pottery dated by Dr. Paul Martin. ) The number of inhabit¬ 
ants per room is a matter of some controversy but we will use the 
estimate given in Zubrow (196 9:40). 


Time 

Number of Rooms 

Occupants/ room 

Population 

100-200 

22 1/2 

5.13 

115 

200-300 

35 

5.13 

179 

300-400 

37 1/2 

5.13 

193 

400-500 

70 

5.13 

359 

500-600 

85 

5.13 

436 

600-700 

110 

5.13 

565 

700-800 

80 

5.13 

410 

800-900 

157 1/2 

5.13 

805 

900-1000 

357 1/2 

5.13 

183 0 

1000-1100 

412 1/2 

5.13 

2115 

1100-1200 

430 1/2 

5.13 

2209 

1200-1300 

497 1/2 

5.13 

2558 


Figure 4. Population of the Hay Hollow Valley 

b. The fraction of mature individuals in the labor pool is esti¬ 
mated at 1/2. This assumes there is a strict male/female division of 
labor for work on energy production. This assumption is based on 
Pueblo ethnographic literature (e.g. Dozier, 1966). 

c. Dozier (1966) implies that Hopi boys start to go to the fields 
with their fathers at age 7 or 8. Based on this we will set 10 as the 
age for productively entering the work force. 

d. The average age at which a worker leaves the work force is 
estimated to be 42. 

e. The maximum birthrate allowable will be set at 3 1/3 percent 
which is within the bounds of birthrates of developing countries. 

f. We estimate initial birth and deathrate to be 2 percent. This 
assumption is noncrucial, in that birthrate easily adjusts. 
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3. Human Energy Variables 

activity by fn^ndi^dull^I’ 24 "ho^T^ USed in seden tary 

(Clarkfl970M 7 r hour day, is set at 2000 Kcal/day 

(above s b e'deltart e e ° UtpU ' ener ^> the ene ^y expended in active work 
(above sedentary energy) is set at 170 Kcal/hr (Clark, 1970: 17). 

tiannf l m The hours worked per day is estimated at 4, in considera¬ 
tion of the arguments of Clark et ah (1970:15) and the short agricul¬ 
tural season of perhaps 4 months. snort agricul- 

4 " Technol ogy, Investment and Savings 
arbitrarily * ‘ ^ WhiCh> “ «“ »• -v.d ia set 


output. 


b* Tiie minimum savings is arbitrarily 5 


percent of the last 


c. 


The level of growth deemed stagnant is .1 percent. 

. d ' . We assume that investment per output unit in hunting and 

loo m TA 1S ■ i>e - L C1 “ k (1970:39) indicates it may require 

man hours/hectare to clear Danish farmland with primitive loots 
I a square meter of land yields 91. 29 Kcal (Zubrow 1971b) and if 

Ht unit of W ° rk r6preSentS 170 Kcals - g«t an investment’per out- 


170 Kcal/hr * 200 hr/hectare 

91.29 Kcal/m a *10 4 m 1 hectare * 


05 


34000 

45000 


= . 75 
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deviations in one or another may not seriously effect the operation of 
a i it -mav well be, however, that the combined error 
sufficient to destroy the predictive power of the model. There was^ 
insufficient time and money available to determine a range o va 
tion of parameters which yield roughly the same results. 

Results 

The simulation was run for 1300 time periods (57 second .^ 
Stanford's IBM 360/67). Output and population onto,, were generated 
from the input parameters specified in the last section. 

Basically, at A.D. 0, the model was started at population 100 with 
^ ' t first stable period after A.D. 

a hunting and gathering economy. JT** per cent increase in technology, 
mo agriculture was introduced, with a 1U percent _ 

In the first stable period after A.D. 1000. irrigation was introduced 
a resource expansion allowing 20 percent better utilisation of the poten 

tial. 

The results of this simulation are given in Appendix B. M 
prediction which can be checked directly is population a 100 year time 
intervals. The real data and the prediction are plotted m Figure 

I think overanalysis of these results would not be worthwhile. 
AsstLtag the data are correct, there are definite problems with the 
model. The population function is so important for outpu 
Turi n population almost certainly Implies failure in prediction of 

output. 


Conclusion 

Dealing with a theory in terms of a simulation forced two things. 
Firsthand most important, it forced the construction of a reasonably 
complete model. Having to link the model, together sufficren y - well 
t w 1 operated in a remotely sensible fashion required directly deal 

ing with many assumptions about the operation of economic and social 
mg witn m y * of the mode l to the data influenced the con¬ 

struction of the model in a way that made it more testable than it might 
otherwise have been. 

The model presented in this paper was certainly not a spectacular 
success. From an optimist's standpoint, its major weakness seems., 
to be its population regulating procedure. This procedure 
problems in the debugging and use of the simulation. yie 
which were intuitively problematic and which seeme oversen * 

With a more successful model of population growth it may be worth¬ 
while to reexamine Schumpeter's ideas m the framewor presen 

here. 

















Population 



Figure 5. Predictions of the Model, and Empirical Data 










- 210 - 


On the other hand, the model may be just plain incorrect from a 
theoretical standpoint. Along this line, George Cowgill has provided 
some thought-provoking criticism (Cowgill; personal communication). 

He suggests first that such "developments" as the replacement of 
hunting and collecting by rainfall agriculture, and its replacement by 
irrigation may in face worsen the ratio of labor input to food output. 

Thus, these changes would not be accepted on a random, or automatic 
basis; rather the advantages or pressures which mitigate for the ac¬ 
ceptance of such "developments" should be studied. 

The foregoing discussion brings up more general points which con¬ 
cern some basic assumption's for this model. The model assumes 
first, an unrelenting tendency for economic growth, and second, a 
similar tendency for the population to increase to the point that the 
incremental population consumes the incremental production. The 
former is an element of Schumpeter's model, the latter, an element 
of this particular formalization of his model. 

Finally, while a model of growth, not surprisingly, predicts growth, 
perhaps the cultural problem should be removed one level to ask the 
answer the "Why" as well as the "How" question of growth. 
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APPENDIX A 

SIMULATION PROGRAM LISTING 
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